Changeset 7041 for branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2016-10-19T12:55:07+02:00 (8 years ago)

Author:

cetlod

Message:

ROBUST5_CNRS : implementation of part I of new TOP interface - 1st step -, see ticket #1782

File:

• branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (3 diffs)

branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -59 +59 @@
       IF(lwp) WRITE(numout,*) '~~~~~~~'
 
+      IF( .NOT.lk_offline )   CALL trc_nam_run     ! Parameters of the run 
       !
       CALL top_alloc()   ! allocate TOP arrays
       !
       CALL trc_ini_ctl   ! control 
-      !
       CALL trc_nam       ! read passive tracers namelists
       !
+      !
       IF(lwp) WRITE(numout,*)
-      !
       IF( ln_rsttr .AND. .NOT. lk_offline ) CALL trc_rst_cal( nit000, 'READ' )   ! calendar
-      !
       IF(lwp) WRITE(numout,*)
       !
       CALL trc_ini_sms   ! SMS
-      !
       CALL trc_ini_trp   ! passive tracers transport
-      !
       CALL trc_ice_ini   ! Tracers in sea ice
       !
-      IF( lwp )  &
-         &  CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea )
+      IF(lwp) CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea )
       !
       CALL trc_ini_state  !  passive tracers initialisation : from a restart or from clim
-      !
-      IF( nn_dttrc /= 1 )        CALL trc_sub_ini      ! Initialize variables for substepping passive tracers
-      !
+      IF( nn_dttrc /= 1 ) &
+      CALL trc_sub_ini      ! Initialize variables for substepping passive tracers
       CALL trc_ini_inv   ! Inventories
       !
@@ -120 +115 @@
          cvol(:,:,jk) = e1e2t(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk)
       END DO
-      IF( lk_degrad )   cvol(:,:,:) = cvol(:,:,:) * facvol(:,:,:)    ! degrad option: reduction by facvol
       !                                                              ! total volume of the ocean 
       areatot = glob_sum( cvol(:,:,:) )
@@ -158 +152 @@
       !! ** Purpose :   SMS initialisation
       !!----------------------------------------------------------------------
-      USE trcini_cfc      ! CFC      initialisation
-      USE trcini_pisces   ! PISCES   initialisation
-      USE trcini_c14b     ! C14 bomb initialisation
-      USE trcini_my_trc   ! MY_TRC   initialisation
-      !!----------------------------------------------------------------------
-      IF( lk_pisces  )       CALL trc_ini_pisces       ! PISCES  bio-model
-      IF( lk_cfc     )       CALL trc_ini_cfc          ! CFC     tracers
-      IF( lk_c14b    )       CALL trc_ini_c14b         ! C14 bomb  tracer
-      IF( lk_my_trc  )       CALL trc_ini_my_trc       ! MY_TRC  tracers
+      USE par_trc
+      USE trcini_pisces  ! PISCES   initialisation
+      USE trcini_cfc     ! CFC      initialisation
+      USE trcini_c14     ! C14  initialisation
+      USE trcini_age     ! age initialisation
+      USE trcini_my_trc  ! MY_TRC   initialisation
+      !!----------------------------------------------------------------------
+      jp_lc = 0
+      IF( lk_pisces  )  THEN       ! PISCES bio-model
+         CALL trc_ini_pisces    
+         jp_lc       = jp_pisces 
+      ELSE
+         jp_pisces   = 0
+         jp_pcs0     = 0              
+         jp_pcs1     = 0         
+      ENDIF
+      IF( lk_my_trc )  THEN      ! MY_TRC
+         CALL trc_ini_my_trc
+         jp_lc       = jp_lc + jp_my_trc  
+      ELSE
+         jp_my_trc   = 0
+         jp_myt0     = 0       
+         jp_myt1     = 0      
+      ENDIF
+      IF( ll_cfc    )  THEN          ! CFC
+         CALL trc_ini_cfc    
+         jp_lc       = jp_lc + jp_cfc
+      ELSE
+         jp_cfc      = 0
+         jp_cfc0     = 0       
+         jp_cfc1     = 0      
+      ENDIF
+      IF( ln_c14    )  THEN      ! C14
+         CALL trc_ini_c14
+         jp_lc       = jp_lc + jp_c14
+      ELSE
+         jp_c14      = 0
+         jp_c140     = 0       
+      ENDIF
+      IF( ln_age    )  THEN      ! AGE
+         CALL trc_ini_age
+         jp_lc       = jp_lc + jp_age
+      ELSE
+         jp_age      = 0
+         jp_age0     = 0       
+      ENDIF
+      !
+      IF( jp_lc /= jptra ) THEN
+        IF( lwp ) WRITE(numout,*) ' The total number of tracers in namelist_top jptra = ', jptra
+        IF( lwp ) WRITE(numout,*) ' The total number of tracers to be used by the model is jp_lc = ', jp_lc
+        CALL ctl_stop( " Change the value of total number of tracers in namelist " )
+      ENDIF
       !
    END SUBROUTINE trc_ini_sms