Commit f80df24f75318c90da276f58bd8ea1e5076d0461

Authored by dumoda01
1 parent c91f14d2
Exists in master and in 1 other branch snow

Certains fichiers vides effaces

BUGS deleted
ChangeLog deleted
HISTORY deleted
... ... @@ -1,2 +0,0 @@
1   -Please read information and documentation on www.gotm.net
2   -
compilers/compiler.GFORTRAN
... ... @@ -6,15 +6,15 @@ FC=gfortran
6 6 DEFINES += -DFORTRAN95
7 7 can_do_F90 = true
8 8 F90_to_f90 = $(FC) -E $(F90FLAGS) $(EXTRA_FFLAGS) $< > $@
9   -#MODULES = -M $(MODDIR)
10   -MODULES = -J $(MODDIR)
  9 +MODULES = -M $(MODDIR)
  10 +#MODULES = -J $(MODDIR)
11 11 EXTRAS = -ffree-form -Wall -std=f95 -pedantic
12 12 DEBUG_FLAGS = -g
13 13 PROF_FLAGS = -p
14 14 PROD_FLAGS = -O0
15 15 DEFINES += -DREAL_4B=real\(4\)
16 16  
17   -NETCDFINC = /opt/local/netcdf-3.6.1/include
18   -NETCDFLIBDIR = /opt/local/netcdf-3.6.1/lib
  17 +#NETCDFINC = /opt/local/netcdf-3.6.1/include
  18 +#NETCDFLIBDIR = /opt/local/netcdf-3.6.1/lib
19 19  
20 20 endif
... ...
compilers/compiler.IFORT
1 1 # Set options for the Intel Fortran 95 compiler - ver. 8.
  2 +
2 3 ifeq ($(FORTRAN_COMPILER),IFORT)
3 4 FC=ifort
4   -DEFINES += -DFORTRAN95
5   -can_do_F90=true
6   -F90_to_f90=
7   -MODULES=
8   -MODULES=-module $(MODDIR)
9   -EXTRAS = -w95 -e95
10   -EXTRAS = -w95
  5 +
  6 +DEFINES += -DFORTRAN95
  7 +can_do_F90 =true
  8 +F90_to_f90 =
  9 +MODULES =
  10 +MODULES =-module $(MODDIR)
  11 +#EXTRAS = -w95 -e95
  12 +EXTRAS = -w95
11 13 DEBUG_FLAGS = -g -C
12 14 PROF_FLAGS = -qp -p
13 15 PROD_FLAGS = -O3
14   -DEFINES += -DREAL_4B=real\(4\)
  16 +DEFINES += -DREAL_4B=real\(4\)
  17 +
15 18 endif
... ...
scripts/ferret/SCM_contribution_to_total_PP.jnl
... ... @@ -8,7 +8,9 @@ use amdgulf
8 8  
9 9 !Variables definition.
10 10  
11   -define symbol alpha=0.04
  11 +define symbol D=15-may-2008
  12 +
  13 +define symbol alpha=0.11
12 14 define symbol inib=0.001
13 15 define symbol vp=0.3
14 16  
... ... @@ -16,18 +18,20 @@ define symbol p0=0.0001
16 18 define symbol k1=0.2
17 19 define symbol k2=0.8
18 20  
19   -define symbol C_N=8.82
20   -define symbol C_g_mol=12
  21 +define symbol C_mg_mmol=12.001
  22 +define symbol N_mg_mmol=14.007
  23 +define symbol C_N=(8.82/(($C_mg_mmol)/($N_mg_mmol)))
  24 +define symbol N_mmol_Chla=(15.54/($N_mg_mmol))
21 25  
22   -define symbol area=200000*1000^2
  26 +define symbol area=450000*1000^2
23 27  
24 28 define symbol tera=10^12
25 29  
26 30  
27 31  
28   -!Find the maximum concentration of phy (excluding the accumulation in the bottom layer; k=6:80)
  32 +!Find the maximum concentration of phy (if excluding the accumulation in the bottom layer; k=6:80)
29 33  
30   -let max_phy=phy[k=6:80@max]
  34 +let max_phy=phy[k=1:80@max]
31 35  
32 36 !plot max_phy
33 37 !frame/file=phy_concentration_SCM.gif
... ... @@ -38,11 +42,15 @@ let max_phy=phy[k=6:80@max]
38 42  
39 43 let norm_phy=phy/max_phy
40 44  
  45 +!fill norm_phy[k=6:80]
  46 +!frame/file=norm_phy.gif
  47 +
41 48  
42 49  
43 50 !Determination of the depth of the SCM.
44 51  
45 52 let depth_SCM=norm_phy[k=1:80@loc:1]
  53 +
46 54 !plot depth_SCM
47 55 !frame/file=depth_SCM.gif
48 56  
... ... @@ -53,6 +61,15 @@ let depth_SCM=norm_phy[k=1:80@loc:1]
53 61 let d_nit=nit[k=1:80@ddc]
54 62 let nitracline=d_nit[k=40:80@loc:-0.02]
55 63  
  64 +fill nit
  65 +frame/file=Nitrate.gif
  66 +
  67 +
  68 +
  69 +!The centered derivative of phy (excluding the accumulation in the bottom layer; k=6:80@ddc).
  70 +
  71 +let d_phy=phy[k=6:80@ddc]
  72 +
56 73  
57 74  
58 75 !The centered derivative of normalized phy (excluding the accumulation in the bottom layer; k=6:80@ddc).
... ... @@ -61,24 +78,54 @@ let d_norm_phy=norm_phy[k=6:80@ddc]
61 78  
62 79  
63 80  
64   -!Absolute values of d_norm_phy. To avoid "s" profiles.
  81 +!Absolute values of d_phy and d_norm_phy. To avoid "s" profiles.
  82 +
  83 +let abs_d_phy=abs(d_phy)
65 84  
66 85 let abs_d_norm_phy=abs(d_norm_phy)
67 86  
  87 +
68 88 !plot /T=($D) abs_d_norm_phy
69 89 !frame/file=Derivative_of_normalized_phy_absolute.gif
70 90  
71 91  
72 92  
  93 +!Smooth of of abs_d_phy to include norm_phy EQ 1.
  94 +
  95 +let smooth_abs_d_phy=abs_d_phy[k=6:80@sbx:8]
  96 +
  97 +
  98 +
73 99 !Graph of phytoplankton, SCM and nitracline depths and definition of the SCM zone.
74 100  
75   -fill phy
  101 +!Phytoplankton mmol N m-3
  102 +
  103 +fill /key=continous /pal=no_red phy[k=40:80]
  104 +
  105 +contour/overplot norm_phy EQ 1
  106 +!contour/overplot d_nit[k=40:80] EQ -0.02
  107 +contour/overplot phy[k=6:80] GE 0.1 and smooth_abs_d_phy[k=6:80] GE 0.002
  108 +
  109 +frame/file=phytoplankton_mmol_N_m-3.gif
  110 +
  111 +
  112 +!Phytoplankton ug chla l-1.
  113 +
  114 +let Chla_phy=phy/($N_mmol_Chla)
  115 +
  116 +fill /key=continous /pal=no_red Chla_phy[k=40:80]
76 117  
77 118 contour/overplot norm_phy EQ 1
78   -contour/overplot d_nit[k=40:80] EQ -0.02
79   -contour/overplot phy[k=6:80] GE 0.003 and abs_d_norm_phy[k=6:80] GE 0.005 or norm_phy[k=6:80] EQ 1
  119 +contour/overplot phy[k=6:80] GE 0.1 and smooth_abs_d_phy[k=6:80] GE 0.002
  120 +
  121 +frame/file=phytoplankton_chla_ug_l.gif
80 122  
81   -frame/file=phytoplankton.gif
  123 +
  124 +
  125 +plot /t=($D) phy
  126 +plot/overlay /t=($D) phy[k=6:80] GE 0.1 and smooth_abs_d_phy[k=6:80] GE 0.01
  127 +
  128 +frame/file=phy_daily.gif
82 129  
83 130  
84 131  
... ... @@ -131,9 +178,10 @@ let f_ratio = new_prod/tot_prod
131 178  
132 179 !Definition of the primary production at the SCM.
133 180  
134   -let good_locations=phy[k=6:80] GE 0.003 and abs_d_norm_phy[k=6:80] GE 0.005 or norm_phy[k=6:80] EQ 1
  181 +let good_locations=phy[k=6:80] GE 0.1 and smooth_abs_d_phy[k=6:80] GE 0.002
135 182  
136   -fill /level=(0,1,0.5) /pal=grayscale good_locations
  183 +fill /level=(0,1,0.5) /pal=grayscale good_locations[k=40:80]
  184 +!contour/overplot norm_phy EQ 1
137 185 frame/file=Good_locations.gif
138 186  
139 187 let SCM_prod=tot_prod*good_locations
... ... @@ -161,28 +209,51 @@ let PP_new_ratio=i_SCM_new_prod/i_new_prod
161 209 let PP_reg_ratio=i_SCM_reg_prod/i_reg_prod
162 210  
163 211  
164   -!Integrated primary production and contribution of the SCM.
  212 +!Integrated primary production and contribution of the SCM with hourly mean radiative fluxes.
165 213  
166   -plot i_prod
167   -plot/overlay i_SCM_prod
168   -plot/overlay PP_ratio
  214 +plot i_prod[l=1:2920@sbx:8]
  215 +plot/overlay i_SCM_prod[l=1:2920@sbx:8]
  216 +plot/overlay PP_ratio[l=1:2920@sbx:8]
169 217  
170 218 frame/file=Integrated_PP.gif
171 219  
172   -plot i_new_prod
173   -plot/overlay i_SCM_new_prod
174   -plot/overlay PP_new_ratio
  220 +plot i_new_prod[l=1:2920@sbx:8]
  221 +plot/overlay i_SCM_new_prod[l=1:2920@sbx:8]
  222 +plot/overlay PP_new_ratio[l=1:2920@sbx:8]
175 223  
176 224 frame/file=Integrated_PP_new.gif
177 225  
178   -plot i_reg_prod
179   -plot/overlay i_SCM_reg_prod
180   -plot/overlay PP_reg_ratio
  226 +plot i_reg_prod[l=1:2920@sbx:8]
  227 +plot/overlay i_SCM_reg_prod[l=1:2920@sbx:8]
  228 +plot/overlay PP_reg_ratio[l=1:2920@sbx:8]
181 229  
182 230 frame/file=Integrated_PP_reg.gif
183 231  
184 232  
185 233  
  234 +!Integrated primary production and contribution of the SCM with daily mean radiative fluxes.
  235 +
  236 +!plot i_prod
  237 +!plot/overlay i_SCM_prod
  238 +!plot/overlay PP_ratio
  239 +
  240 +!frame/file=Integrated_PP.gif
  241 +
  242 +!plot i_new_prod
  243 +!plot/overlay i_SCM_new_prod
  244 +!plot/overlay PP_new_ratio
  245 +
  246 +!frame/file=Integrated_PP_new.gif
  247 +
  248 +!plot i_reg_prod
  249 +!plot/overlay i_SCM_reg_prod
  250 +!plot/overlay PP_reg_ratio
  251 +
  252 +!frame/file=Integrated_PP_reg.gif
  253 +
  254 +
  255 +
  256 +
186 257 !Annual primary production.
187 258  
188 259 let annual_PP=i_prod[l=1:2920@sum]
... ... @@ -241,13 +312,13 @@ let C_SCM_annual_PP_reg=SCM_annual_PP_reg*($C_N)
241 312  
242 313 !To mmol C m-2 y-1 to g C m-2 y-1
243 314  
244   -let C_annual_PP_g=C_annual_PP*($C_g_mol)/1000
245   -let C_annual_PP_new_g=C_annual_PP_new*($C_g_mol)/1000
246   -let C_annual_PP_reg_g=C_annual_PP_reg*($C_g_mol)/1000
  315 +let C_annual_PP_g=C_annual_PP*($C_mg_mmol)/1000
  316 +let C_annual_PP_new_g=C_annual_PP_new*($C_mg_mmol)/1000
  317 +let C_annual_PP_reg_g=C_annual_PP_reg*($C_mg_mmol)/1000
247 318  
248   -let C_SCM_annual_PP_g=C_SCM_annual_PP*($C_g_mol)/1000
249   -let C_SCM_annual_PP_new_g=C_SCM_annual_PP_new*($C_g_mol)/1000
250   -let C_SCM_annual_PP_reg_g=C_SCM_annual_PP_reg*($C_g_mol)/1000
  319 +let C_SCM_annual_PP_g=C_SCM_annual_PP*($C_mg_mmol)/1000
  320 +let C_SCM_annual_PP_new_g=C_SCM_annual_PP_new*($C_mg_mmol)/1000
  321 +let C_SCM_annual_PP_reg_g=C_SCM_annual_PP_reg*($C_mg_mmol)/1000
251 322  
252 323  
253 324 !To g C m-2 y-1 to Tg C y-1
... ...