Changeset 5601 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90

Timestamp:

2015-07-16T11:04:29+02:00 (9 years ago)

Author:

cbricaud

Message:

commit changes/bugfix/... for crs ; ok with time-splitting/fixed volume

File:

• branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90 (modified) (2 diffs)

branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90

@@ -16 +16 @@
    USE trc
    USE trcsms_my_trc
+   USE dom_oce, ONLY : gdepw_1d,e3t_1d,nyear_len
 
    IMPLICIT NONE
@@ -44 +45 @@
       IF(lwp) WRITE(numout,*) ' trc_ini_my_trc: initialisation of MY_TRC model'
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
+
+      nlb_age = MINLOC( gdepw_1d, mask = gdepw_1d > age_depth, dim = 1  )
+      nl_age = nlb_age - 1
+      nla_age = nl_age - 1
+      frac_kill_age = (age_depth - gdepw_1d(nl_age))/e3t_1d(nl_age)
+      frac_add_age = 1._wp -  frac_kill_age
+      rryear = 1._wp / ( nyear_len(1) * rday )
+
       
+      IF( .NOT. ln_rsttr ) trb(:,:,:,jp_myt0:jp_myt1) = 0.
       IF( .NOT. ln_rsttr ) trn(:,:,:,jp_myt0:jp_myt1) = 0.
       !