Changes between Version 6 and Version 7 of ticket/0846

Timestamp:

2011-09-12T15:12:36+02:00 (13 years ago)

Author:

cetlod

Comment:

--

ticket/0846

@@ -276 +276 @@
 }}}
 
+Christian's comments
+
+The NEMO architecture ( with the use of FCM to build dependancies ), allow to rewrite a routine without needing to change anything. In that case, we can use the MY_TRC framework to manage every type of biogeochemical model.
+
+In fact, when compiling with makenemo ( for example : ./makenemo -n ORCA2_LIM_MY_TRC -r ORCA2_LIM_PISCES del_key "key_iomput key_pisces" add_key "key_my_trc"), a sub-directory where we can copy and modify all the routines we need is created : NEMOGCM/CONFIG/ORCA2_LIM_MY_TRC/MY_SRC 
+
+ * trcini.F90
+ * trcstp.F90
+ * trcdia.F90
+ * trctrp.F90
+ * MY_TRC/*
+
+To go a little further we can add a new namelits parameter ln_nemo_io which is true if we want to use NEMO I/O manager and false if not. So we can modify the 2 modules trcini.F90 and trcstp.F90 as follow :
+
+trcini.F90
+{{
+      IF( lk_lobster ) THEN   ;   CALL trc_ini_lobster      ! LOBSTER bio-model
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          LOBSTER not used'
+      ENDIF
+
+      IF( lk_pisces  ) THEN   ;   CALL trc_ini_pisces       ! PISCES  bio-model
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          PISCES not used'
+      ENDIF
+
+      IF( lk_cfc     ) THEN   ;   CALL trc_ini_cfc          ! CFC     tracers
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          CFC not used'
+      ENDIF
+
+      IF( lk_c14b    ) THEN   ;   CALL trc_ini_c14b         ! C14 bomb  tracer
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          C14 not used'
+      ENDIF
+
+      IF( lk_my_trc  ) THEN   ;   CALL trc_ini_my_trc       ! MY_TRC  tracers
+      ELSE                    ;   IF(lwp) WRITE(numout,*) '          MY_TRC not used'
+      ENDIF
+
+      IF( ln_nemo_io ) THEN               ! Uses NEMO I/O and restart manager 
+         !
+         IF( ln_rsttr ) THEN
+           !
+           IF( lk_offline )  neuler = 1   ! Set time-step indicator at nit000 (leap-frog)
+           CALL trc_rst_read              ! restart from a file
+           !
+         ELSE
+           IF( lk_offline )  THEN
+              neuler = 0                  ! Set time-step indicator at nit000 (euler)
+              CALL day_init               ! set calendar
+           ENDIF
+#if defined key_dtatrc
+              CALL trc_dta( nit000 )      ! Initialization of tracer from a file that may also be used for damping
+              DO jn = 1, jptra
+              IF( lutini(jn) )   trn(:,:,:,jn) = trdta(:,:,:,jn) * tmask(:,:,:)   ! initialisation from file if required
+           END DO
+#endif
+           trb(:,:,:,:) = trn(:,:,:,:)
+           ! 
+         ENDIF
+         !
+      ENDIF
+
+      tra(:,:,:,:) = 0._wp
+
+
+}}
+
+trcstp.F90
+
+{{
+      IF( MOD( kt - 1 , nn_dttrc ) == 0 ) THEN      ! only every nn_dttrc time step
+         !
+         IF(ln_ctl) THEN
+            WRITE(charout,FMT="('kt =', I4,'  d/m/y =',I2,I2,I4)") kt, nday, nmonth, nyear
+            CALL prt_ctl_trc_info(charout)
+         ENDIF
+         !
+         tra(:,:,:,:) = 0.e0
+         !
+         IF( ln_nemo_io ) THEN                                  ! Uses NEMO I/O and restart manager 
+            !
+            IF( kt == nit000 .AND. lk_trdmld_trc  )  &
+               &                    CALL trd_mld_trc_init        ! trends: Mixed-layer
+                                    CALL trc_rst_opn( kt )       ! Open tracer restart file 
+         ENDIF
+         IF( lk_iomput ) THEN  ;    CALL trc_wri( kt )           ! output of passive tracers
+         ELSE                  ;    CALL trc_dia( kt )
+         ENDIF
+                                    CALL trc_sms( kt )       ! tracers: sink and source
+                                    CALL trc_trp( kt )       ! transport of passive tracers 
+
+         IF( ln_nemo_io ) THEN                                  ! Uses NEMO I/O and restart manager 
+            IF( kt == nit000 )      CALL iom_close( numrtr )     ! close input  passive tracers restart file
+            IF( lrst_trc )          CALL trc_rst_wri( kt )       ! write tracer restart file
+            IF( lk_trdmld_trc  )    CALL trd_mld_trc( kt )       ! trends: Mixed-layer
+         ENDIF
+         !
+      ENDIF
+}}
 === Status ===
 On the dev branch we have now committed the method we currently use and that has to be modified. It is just a placeholder for the status quo at CMCC and very similar to NERC and MetOffice.