Changeset 3897 for branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90

Timestamp:

2013-05-02T16:29:49+02:00 (11 years ago)

Author:

cetlod

Message:

2013/dev_r3411_CNRS4_IOCRS: 2nd step create a new configuration with age tracer ; starting point for TOP coarsening

File:

• branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90 (modified) (1 diff)

branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90

@@ -44 +44 @@
       IF(lwp) WRITE(numout,*) ' trc_ini_my_trc: initialisation of MY_TRC model'
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
+
+!! BUG in s-coordinate this does not work!
+      ! deepest/shallowest W level Above/Below age_depth
+      nlb_age = MINLOC( gdepw_0, mask = gdepw_0 > age_depth, dim = 1 )   ! shallowest W level Below age_depth
+                                                                         !  = shallowest T level wholly below age_depth
+      nl_age = nlb_age - 1                                               ! deepest    W level Above age_depth
+                                                                         !  = T level surrounding age_depth
+      nla_age = nl_age - 1                                               ! deepest    T level wholly above age_depth
+      frac_kill_age = (age_depth - gdepw_0(nl_age))/e3t_0(nl_age)         ! fraction of level nl_age above age_depth
+      frac_add_age = 1._wp -  frac_kill_age                              ! fraction of level nl_age below age_depth
+!! end bug
+
+      rryear = 1._wp / ( nyear_len(1) * rday )
       
       IF( .NOT. ln_rsttr ) trn(:,:,:,jp_myt0:jp_myt1) = 0.