Changeset 8213 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2017-06-25T11:27:43+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- split trcbio - mergeable-ish MEDUSA branch

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (10 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -24 +24 @@
    USE trcini_pisces   ! PISCES   initialisation
    USE trcini_c14b     ! C14 bomb initialisation
+   USE trcini_age      ! AGE      initialisation
    USE trcini_my_trc   ! MY_TRC   initialisation
+   USE trcini_idtra    ! idealize tracer initialisation
    USE trcini_medusa   ! MEDUSA   initialisation
-   USE trcini_idtra    ! idealize tracer initialisation
-   USE trcini_age      ! AGE      initialisation
    USE trcdta          ! initialisation from files
    USE daymod          ! calendar manager
@@ -81 +81 @@
          &   CALL ctl_warn(' Coupling with passive tracers and used of diurnal cycle. &
          & Computation of a daily mean shortwave for some biogeochemical models) ')
-          !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-          !!!!! CHECK For MEDUSA
-          !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
       IF( nn_cla == 1 )   &
          &  CALL ctl_stop( ' Cross Land Advection not yet implemented with passive tracer ; nn_cla must be 0' )
@@ -112 +110 @@
 
       CALL trc_ice_ini                                 ! Tracers in sea ice
-
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' in trc_init'
-         IF (lwp) write (numout,*) ' sms init OK'
-         IF (lwp) write (numout,*) ' next: open tracer.stat'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
 
       IF( ln_ctl ) THEN
@@ -135 +123 @@
       ENDIF
 
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' in trc_init'
-         IF (lwp) write (numout,*) 'open tracer.stat -- OK'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
-
-
       IF( ln_trcdta ) THEN
-#if defined key_medusa
-         IF(lwp) WRITE(numout,*) 'AXY: calling trc_dta_init'
-         IF(lwp) CALL flush(numout)
-#endif
          CALL trc_dta_init(jptra)
       ENDIF
@@ -155 +129 @@
       IF( ln_rsttr ) THEN
         !
-#if defined key_medusa
-        IF(lwp) WRITE(numout,*) 'AXY: calling trc_rst_read'
-        IF(lwp) CALL flush(numout)
-#endif
         CALL trc_rst_read              ! restart from a file
         !
       ELSE
-        !
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' Init from file -- will call trc_dta'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
         !
         IF( ln_trcdta .AND. nb_trcdta > 0 ) THEN  ! Initialisation of tracer from a file that may also be used for damping
@@ -193 +155 @@
         ENDIF
         !
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' in trc_init'
-         IF (lwp) write (numout,*) ' before trb = trn'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
-        !
         trb(:,:,:,:) = trn(:,:,:,:)
-        ! 
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' in trc_init'
-         IF (lwp) write (numout,*) ' trb = trn -- OK'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
         ! 
       ENDIF
@@ -220 +164 @@
       IF( ln_zps .AND. .NOT. lk_c1d .AND.       ln_isfcav )   &
         &    CALL zps_hde_isf( nit000, jptra, trn, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi )       ! tracers at the bottom ocean level
-      !
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' in trc_init'
-         IF (lwp) write (numout,*) ' partial step -- OK'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
+
+
       !
       IF( nn_dttrc /= 1 )        CALL trc_sub_ini      ! Initialize variables for substepping passive tracers
       !
-# if defined key_debug_medusa
-         IF (lwp) write (numout,*) '------------------------------'
-         IF (lwp) write (numout,*) 'Jpalm - debug'
-         IF (lwp) write (numout,*) ' in trc_init'
-         IF (lwp) write (numout,*) ' before initiate tracer contents'
-         IF (lwp) write (numout,*) ' '
-         CALL flush(numout)
-# endif
-      !
+
       trai(:) = 0._wp                                                   ! initial content of all tracers
       DO jn = 1, jptra
@@ -253 +182 @@
          WRITE(numout,*) '          *** Total inital content of all tracers '
          WRITE(numout,*)
-# if defined key_debug_medusa
-         CALL flush(numout)
-# endif
-         !
-# if defined key_debug_medusa
-         WRITE(numout,*) ' litle check :  ', ctrcnm(1)
-         CALL flush(numout)
-# endif
          DO jn = 1, jptra
             WRITE(numout,9000) jn, TRIM( ctrcnm(jn) ), trai(jn)
@@ -304 +225 @@
       USE trdmxl_trc    , ONLY:   trd_mxl_trc_alloc
 #endif
+# if defined key_medusa
+      USE bio_medusa_mod, ONLY:   bio_medusa_alloc
+# endif
+
       !
       INTEGER :: ierr
@@ -316 +241 @@
       ierr = ierr + trd_mxl_trc_alloc()
 #endif
+#if defined key_medusa
+      ierr = ierr + bio_medusa_alloc()
+#endif
       !
       IF( lk_mpp    )   CALL mpp_sum( ierr )