Changeset 7017 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

Timestamp:

2016-10-11T16:20:23+02:00 (8 years ago)

Author:

jpalmier

Message:

Jpalm -- 11-10-2016 -- Improve MEDUSA reproducibility - change carb-chem calling/restart

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (1 diff)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -1039 +1039 @@
       !!                 be calculated and done properly
       !! IF( kt == nit000 .or. mod(kt,1920) == 0 ) THEN
-      IF( kt == nittrc000 .or. mod(kt,960) == 0 ) THEN
+      !! IF( kt == nittrc000 .or. mod(kt,960) == 0 ) THEN 
+      !!=============================
+      !! Jpalm -- 07-10-2016 -- need to change carb-chem frequency call :
+      !!          we don't want to call on the first time-step of all run submission, 
+      !!          but only on the very first time-step, and then every month
+      !!          So we call on nittrc000 if not restarted run, 
+      !!          else if one month after last call.
+      !!          assume one month is 30d --> 3600*24*30 : 2592000s
+      !!          try to call carb-chem at 1st month's tm-stp : x * 30d + 1*rdt(i.e: mod = rdt)   
+      !!          ++ need to pass carb-chem output var through restarts
+      If ( ( kt == nittrc000 .AND. .NOT.ln_rsttr ) .OR. mod(kt*rdt,2592000.) == rdt ) THEN
          !!----------------------------------------------------------------------
          !! Calculate the carbonate chemistry for the whole ocean on the first