Changeset 8280 for branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2017-07-05T10:28:51+02:00 (7 years ago)

Author:

timgraham

Message:

331: Merge of MEDUSA stable branch and HadGEM3 coupling branches into GO6 package branch.

File:

• branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (8 diffs)

branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -8 +8 @@
    !!            2.0  ! 2005-10 (C. Ethe, G. Madec) revised architecture
    !!            4.0  ! 2011-01 (A. R. Porter, STFC Daresbury) dynamical allocation
+   !!             -   ! 2014-06 (A. Yool, J. Palmieri) adding MEDUSA-2
    !!----------------------------------------------------------------------
 #if defined key_top
@@ -24 +25 @@
    USE trcini_c14b     ! C14 bomb initialisation
    USE trcini_my_trc   ! MY_TRC   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
@@ -76 +80 @@
          &   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' )
@@ -97 +103 @@
 
       IF( lk_pisces  )       CALL trc_ini_pisces       ! PISCES  bio-model
+      IF( lk_medusa  )       CALL trc_ini_medusa       ! MEDUSA  tracers
+      IF( lk_idtra   )       CALL trc_ini_idtra        ! Idealize tracers
       IF( lk_cfc     )       CALL trc_ini_cfc          ! CFC     tracers
       IF( lk_c14b    )       CALL trc_ini_c14b         ! C14 bomb  tracer
+      IF( lk_age     )       CALL trc_ini_age          ! AGE       tracer
       IF( lk_my_trc  )       CALL trc_ini_my_trc       ! MY_TRC  tracers
 
       CALL trc_ice_ini                                 ! Tracers in sea ice
 
-      IF( lwp ) THEN
+# 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
          !
-         CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea )
+         IF (narea == 1) THEN  
+            ! The tracer.stat file only contains global tracer sum values, if 
+            ! it contains anything at all. Hence it only needs to be opened 
+            ! and written to on the master PE, not on all PEs.  
+            CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE','FORMATTED',  & 
+                          'SEQUENTIAL', -1, numout, lwp , narea ) 
+         ENDIF  
          !
       ENDIF
 
-      IF( ln_trcdta )      CALL trc_dta_init(jptra)
-
+# 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
 
       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
@@ -137 +189 @@
         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
@@ -146 +216 @@
       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
@@ -164 +249 @@
          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)
@@ -176 +269 @@
          CALL prt_ctl_trc( tab4d=trn, mask=tmask, clinfo=ctrcnm )
       ENDIF
+
+      IF(lwp) WRITE(numout,*)
+      IF(lwp) WRITE(numout,*) 'trc_init : passive tracer set up completed'
+      IF(lwp) WRITE(numout,*) '~~~~~~~'
+      IF(lwp) CALL flush(numout)
+# if defined key_debug_medusa
+         CALL trc_rst_stat
+         CALL flush(numout)
+# endif
+
 9000  FORMAT(' tracer nb : ',i2,'      name :',a10,'      initial content :',e18.10)
       !