Changeset 9335


Ignore:
Timestamp:
2018-02-19T16:29:35+01:00 (2 years ago)
Author:
jpalmier
Message:

Import 3D carb-chem calling conditions for tests in the coupled model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_CarbFix_DMSMin_3DMOCSY/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r8649 r9335  
    8080                                            gdept_0, gdept_n,               & 
    8181                                            gdepw_0, gdepw_n,               & 
    82                                             nday_year, nsec_day, nyear,     & 
    83                                             rdt, tmask, mig, mjg 
     82                                            nday_year, nsec_day,            & 
     83                                            nyear, nyear_len, ndastp,       & 
     84                                            nsec_month,                     & 
     85                                            rdt, tmask, mig, mjg, nimpp,    & 
     86                                            njmpp  
    8487      USE in_out_manager,             ONLY: lwp, numout, nn_date0 
    8588# if defined key_iomput 
     
    353356      !!          x * 30d + 1*rdt(i.e: mod = rdt)    
    354357      !!          ++ need to pass carb-chem output var through restarts 
    355       If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                        & 
    356            ((86400*mod(nn_date0,100) + mod(kt*rdt,2592000.)) == rdt) ) THEN 
     358      !!If ( (kt == nitt8rc000 .AND. .NOT.ln_rsttr) .OR.          & 
     359      !!     ( (mod(kt*rdt,2592000.)) == rdt) THEN 
     360      !!============================= 
     361      !! (Jpalm -- updated for restartability issues) 
     362      !! We want this to be start of month or if starting afresh from   
     363      !! climatology - marc 20/6/17  
     364      !!If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                         &  
     365      !!     ((86400*mod(nn_date0,100) + mod(kt*rdt,2592000.)) == rdt) ) THEN  
     366      !!============================= 
     367      !! Jpalm -- 15-02-2018 -- need to change 3D carb-chem call freq again. 
     368      !! previous call did not work, probably the (86400*mod(nn_date0,100) part 
     369      !! should not be in... 
     370      !! now use the NEMO calendar tool : nsec_month to be sure to call  
     371      !! at the beginning of a new month . 
     372      IF ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                        & 
     373           ( nsec_month .LE. INT(rdt) ) )  THEN 
     374           IF ( lwp )  WRITE(numout,*)                                       &          
     375                              ' *** 3D carb chem call *** -- kt:', kt,       & 
     376                              ';  current date:', ndastp  
    357377         !!--------------------------------------------------------------- 
    358378         !! Calculate the carbonate chemistry for the whole ocean on the first 
Note: See TracChangeset for help on using the changeset viewer.