Changeset 7693 for branches/UKMO/dev_r5518_MEDUSA_optim_RH/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2017-02-17T16:38:26+01:00 (7 years ago)

Author:

frrh

Message:

Merge in branches/NERC/dev_r5518_NOC_MEDUSA_Stable from revision range
r 5711:7611 (i.e. the entire branch)

File:

• branches/UKMO/dev_r5518_MEDUSA_optim_RH/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (9 diffs)

branches/UKMO/dev_r5518_MEDUSA_optim_RH/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
@@ -77 +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' )
@@ -98 +104 @@
 
       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 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( lwp ) THEN
@@ -110 +129 @@
       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
         !
@@ -141 +179 @@
         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
@@ -150 +206 @@
       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
@@ -168 +239 @@
          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)
@@ -180 +259 @@
          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)
       !