update sat/rivers

config/config_MARIA.json → config/config_FUMSECK.json

@@ -1,31 +1,33 @@
 {
 "_____":"GENERAL PARAMETERS",
 
-"date_ini":"2021-12-11",
-"date_end":"2021-12-12",
-"outdir":"/path/to/OUTPUT_SMP",
-"diagdir":"/path/to/diags/dir",
-"river_ini":"2015-01-03",
-"river_end":"2015-05-01",
+"date_ini":"2019-05-01",
+"date_end":"2019-05-15",
+"outdir":"/DISK2TB/MODEL/FUMSECK",
+"diagdir":"/DISK2TB/MODEL/DIAGS",
 
 
 "_____":"RIVER PARAMETERS",
 
 "french_list":["Argens","Var"],
-"italian_list":["Arno","Magra","Ombrone","Serchio"],
+"italian_list":["Arno","Magra","Magra","Serchio","Ombrone"],
+"river_ini":"2019-01-01",
+"river_end":"2019-12-31",
 
 
 "_____":"METEO PARAMETERS",
-"meteo_ini":"2015-03-10",
-"meteo_end":"2015-03-11",
+
+"meteo_ini":"2015-03-31",
+"meteo_end":"2015-05-30",
 "day_cycle":2,
-"n_cycle":6,
+"n_cycle":15,
 "imin":"500",
 "imax":"600",
 "jmin":"350",
 "jmax":"450",
 
-
+"_____":"MODIS PARAMETERS",
+"modis_dir":"/DISK2TB/DATA/NWMED/MODIS",
 
 "_____":"DYFAMED PARAMETERS",
 
@@ -35,47 +37,58 @@
 
 "some_par":"some_par",
 
+"_____":"CMEMS",
+
+
+"_____":"ARGO",
 
-"_____":"2D PLOT PARAMETERS",
+"argo_ds":"phy",
+
+
+"_____":"PLOT PARAMETERS",
 
 "resol":"10m",                       
 "fig_proj":"ccrs.PlateCarree()",      
 "fig_sx":"8",                         
 "fig_sy":"8",                        
 "fig_fmt":"jpg",                     
 "fig_res":"300",                  
-"tight":"False",                  
+"tight":"True",                  
 "fig_tck_size":"10",
 "fig_tcklbl_size":"10",
 "fig_lbl_size":"10",
 "cb_fraction_2D":"0.02",
 "cb_pad_2D":"0.12",
+"with_currents":"False",
+"with_whales":"False",
 
 "_____":"SECTION PLOT PARAMETERS",
 
+
 "fig_secx":"8",                         
 "fig_secy":"4",                        
-"sec_lon1":"8.5",
+"sec_lon1":"9",
 "sec_lat1":"43.",
-"sec_lon2":"10.",
-"sec_lat2":"43.5",
-"sec_dep":"600",
+"sec_lon2":"9.7",
+"sec_lat2":"44",
+"sec_dep":"300",
 "sec_resH":"0.05",  
 "sec_resV":"1",
 "cb_fraction_sec":"0.025",
 "cb_pad_sec":"0.06",
 
 "_____":"PROFILE PLOT PARAMETERS",
 
-"fig_prox":"8",                         
-"fig_proy":"4",                        
-"col1":"r",
-"col2":"b",
+"fig_prox":"6",
+"fig_proy":"11",
+"col_bcg":["g","b","r"],
+"col_phy":["b","r"],
 
 
 "_____":"INTERPOLATION PARAMETERS",
 
 "dump":"0.01",
-"itp_meth":"nearest"
+"itp_meth":"linear",
+"cmap":"viridis"
 
 }

config/variables_mit_FUMSECK.dat

@@ -0,0 +1,18 @@
+# Name tag cmap log vlim vmin vmax label units
+# ---------------------------------------------
+nppv BIOL jet False set 0. 1. lab un
+o2 BIOL jet False set  0. 1.  lab un
+ph CARB jet False set 0. 1. lab un
+talk CARB jet False set 0. 1. lab un
+spco2 CO2F jet False set 0. 1. lab un
+nh4 NUTR jet False set 0. 1. NH4 $mmol.m^{-3}$
+no3 NUTR cmc.hawaii_r False set 0.5 4. NO3 $mmol.m^{-3}$
+po4 NUTR cmc.batlow False set 0. 0.3 PO4 $mmol.m^{-3}$
+si NUTR jet False set 0. 1. SI $mmol.m^{-3}$
+chl PFTC cmc.imola True set 0.05 0.3 Chlorophyll $mg.m^{-3}$
+phyc PFTC jet False set 0. 1. Phytoplankton $mmol.m^{-3}$
+zooc PFTC cmc.hawaii_r False set 1. 3. Zooplankton $mg.m^{-3}$
+so PSAL jet False set 0. 1. Salinity none 
+uo RFVL cmc.vik False set -.3 .3 Eastward_current $m.s^{-1}$
+vo RFVL jet False set 0. 1. Northward_current $m.s^{-1}$
+thetao TEMP cmo.thermal False set 13.5 14.8 Temperature °C

config/variables_mit_TEMPLATE.dat

@@ -9,10 +9,10 @@ nh4 NUTR jet False set 0. 1. NH4 $mmol.m^{-3}$
 no3 NUTR cmc.hawaii_r False set 0.5 4. NO3 $mmol.m^{-3}$
 po4 NUTR cmc.batlow False set 0. 0.3 PO4 $mmol.m^{-3}$
 si NUTR jet False set 0. 1. SI $mmol.m^{-3}$
-chl PFTC cmc.imola True set 0.05 0.25 Chlorophyll $mg.m^{-3}$
+chl PFTC cmc.imola True set 0.05 0.5 Chlorophyll $mg.m^{-3}$
 phyc PFTC jet False set 0. 1. Phytoplankton $mmol.m^{-3}$
-zooc PFTC jet False set 0. 1.2 Zooplankton $mmol.m^{-3}$
+zooc PFTC cmc.broc False set 0. 3. Zooplankton $mg.m^{-3}$
 so PSAL jet False set 0. 1. Salinity none 
-uo RFVL jet False set 0. 1. Eastward_current $m.s^{-1}$
+uo RFVL cmc.vik False set -.3 .3 Eastward_current $m.s^{-1}$
 vo RFVL jet False set 0. 1. Northward_current $m.s^{-1}$
-thetao TEMP cmo.thermal False set 13. 16.5 Temperature °C
+thetao TEMP cmo.thermal False set 13.5 14.8 Temperature °C

datasets/ARGO/compare_argo.py

@@ -117,8 +117,10 @@
         = load_variable(config,mod_var)
 
     # Convert date
-    dstr = str(time[idx][0]).replace('.000000000','')
-    date = dt.datetime.strptime(dstr,'%Y-%m-%dT%H:%M:%S')
+    dstr = str(time[idx][0])
+    dstr = dstr.split(':')[0]
+
+    date = dt.datetime.strptime(dstr,'%Y-%m-%dT%H')
     jd = date.toordinal()
 
     hours,mnt = date.strftime('%H'),date.strftime('%M')
@@ -138,6 +140,7 @@
                  data_lon,data_lat,data_pres)
     except Exception as e:
       print('Interpolation failed')
+      print(e)
       val = data_pres
       val[:] = np.nan
 

datasets/CMEMS/get_data.py

@@ -0,0 +1,6 @@
+import pandas as pd
+
+
+
+
+

datasets/NASA/get_modis.py

@@ -0,0 +1,38 @@
+#############################
+# Routine to get MODIS data #
+#############################
+
+
+import numpy as np
+import os
+
+
+# Get args
+# --------
+config=sys.argv[1]  # Configuration name
+
+
+# Load parameters
+# ---------------
+load_config(config)
+from fun_gen import *
+
+os.system('mkdir -p '+modis_dir)
+
+
+# Loop on days
+for jd in range(jdini,jdend+1):
+
+    # Get date
+    d = dt.datetime.fromordinal(jd)
+    year = d.strftime('%Y')
+    doy = d.strftime('%j').zfill(3)
+
+    # Download
+    # !!!! TO FINISH !!!!!
+
+
+
+
+
+