New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 6201 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable – NEMO

Ignore:
Timestamp:
2016-01-04T14:06:11+01:00 (8 years ago)
Author:
jpalmier
Message:

JPALM -- 04-01-2016 -- add debugg prints under debugg_key

Location:
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r6164 r6201  
    200200             WRITE(NUMOUT,*) ' ' 
    201201         END DO   
     202# if defined key_debug_medusa 
     203         CALL flush(numout) 
     204# endif 
    202205      ENDIF 
    203206        ! 

  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcsms_idtra.F90

    r6164 r6201  
    8282      !!---------------------------------------------------------------------- 
    8383      IF (kt == nittrc000) THEN 
    84          IF(lwp) WRITE(numout,*) '   - JPALM - verif :' 
     84         IF(lwp) WRITE(numout,*) '   trcsms_idtra :' 
    8585         IF(lwp) WRITE(numout,*) '   ~~~~~~~~~~~~~~~~~' 
    8686         IF(lwp) WRITE(numout,*) '   - idtra decay factor : ', FDEC 
     87# if defined key_debug_medusa 
     88      CALL flush(numout) 
     89# endif 
    8790      ENDIF 
    8891 
    … …  
    160163      IF( nn_timing == 1 )  CALL timing_stop('trc_sms_idtra') 
    161164      ! 
     165# if defined key_debug_medusa 
     166         IF(lwp) WRITE(numout,*) '   IDTRA DONE trc_sms_idtra -- ' 
     167      CALL flush(numout) 
     168# endif 
     169      ! 
    162170   END SUBROUTINE trc_sms_idtra 
    163171 

  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcnam.F90

    r6164 r6201  
    270270      ELSE                    ;   IF(lwp) WRITE(numout,*) '          CFC not used' 
    271271      ENDIF 
    272  
     272      ! 
     273# if defined key_debug_medusa 
     274      CALL flush(numout) 
     275      IF (lwp) write (numout,*) '------------------------------' 
     276      IF (lwp) write (numout,*) 'Jpalm - debug' 
     277      IF (lwp) write (numout,*) 'CALL trc_nam_cfc -- OK' 
     278      IF (lwp) write (numout,*) 'in trc_nam - CALL trc_nam OK' 
     279      IF (lwp) write (numout,*) ' ' 
     280# endif 
     281      ! 
    273282      IF( lk_c14b     ) THEN   ;   CALL trc_nam_c14b         ! C14 bomb     tracers 
    274283      ELSE                    ;   IF(lwp) WRITE(numout,*) '          C14 not used' 

  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r6174 r6201  
    8686      IF( kt == nitrst - 2*nn_dttrc .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc .AND. .NOT. lrst_trc ) ) THEN 
    8787         IF ( ln_rstdate ) THEN 
    88             !! JPALM -- 22-12-2015 -- modif to get the good date on restart trc 
    89             !file name 
    90             !!                     -- the condition to open the rst file is not 
    91             !the same than for the dynamic rst. 
    92             !!                     -- here it - for an obscure reason - is open 
    93             !2 time-step before the restart writing process 
     88            !! JPALM -- 22-12-2015 -- modif to get the good date on restart trc file name 
     89            !!                     -- the condition to open the rst file is not the same than for the dynamic rst. 
     90            !!                     -- here it - for an obscure reason - is open 2 time-step before the restart writing process 
    9491            !!                     instead of 1. 
    95             !!                     -- i am not sure if someone forgot +1 in the 
    96             !if loop condition as 
    97             !!                     it is writen in all comments nitrst - 
    98             !2*nn_dttrc + 1 and the condition is nitrst - 2*nn_dttrc 
    99             !!                     -- nevertheless we didn't wanted to broke 
    100             !something already working  
     92            !!                     -- i am not sure if someone forgot +1 in the if loop condition as 
     93            !!                     it is writen in all comments nitrst -2*nn_dttrc + 1 and the condition is  
     94            !!                     nitrst - 2*nn_dttrc 
     95            !!                     -- nevertheless we didn't wanted to broke something already working  
    10196            !!                     and just adapted the part we added. 
    102             !!                     -- So instead of calling ju2ymds( fjulday + 
    103             !(rdttra(1)) .... 
    104             !!                     we call ju2ymds( fjulday + (2*rdttra(1)) 
    105             !..... 
     97            !!                     -- So instead of calling ju2ymds( fjulday + (rdttra(1))  
     98            !!                     we call ju2ymds( fjulday + (2*rdttra(1))  
    10699            !!--------------------------------------------------------------------       
    107100            CALL ju2ymds( fjulday + (2*rdttra(1)) / rday, iyear, imonth, iday, zsec ) 
    … …  
    219212      if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 
    220213         &        fq0, fq1, fq2 
     214      !! 
     215      !! AXY (07/07/15): read in temporally averaged fields for DMS 
     216      !!                 calculations 
     217      !! 
     218      IF( iom_varid( numrtr, 'B_DMS_CHN', ldstop = .FALSE. ) > 0 ) THEN 
     219         !! YES; in which case read them 
     220         !! 
     221         IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS present - reading in ...' 
     222         CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_CHN',  zb_dms_chn(:,:)  ) 
     223         CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_CHN',  zn_dms_chn(:,:)  ) 
     224         CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_CHD',  zb_dms_chd(:,:)  ) 
     225         CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_CHD',  zn_dms_chd(:,:)  ) 
     226         CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_MLD',  zb_dms_mld(:,:)  ) 
     227         CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_MLD',  zn_dms_mld(:,:)  ) 
     228         CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_QSR',  zb_dms_qsr(:,:)  ) 
     229         CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_QSR',  zn_dms_qsr(:,:)  ) 
     230         CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_DIN',  zb_dms_din(:,:)  ) 
     231         CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_DIN',  zn_dms_din(:,:)  ) 
     232      ELSE 
     233         !! NO; in which case set them to zero 
     234         !! 
     235         IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS absent - setting to zero ...' 
     236         zb_dms_chn(:,:)  = 0.0   !! CHN 
     237         zn_dms_chn(:,:)  = 0.0 
     238         zb_dms_chd(:,:)  = 0.0   !! CHD 
     239         zn_dms_chd(:,:)  = 0.0 
     240         zb_dms_mld(:,:)  = 0.0   !! MLD 
     241         zn_dms_mld(:,:)  = 0.0 
     242         zb_dms_qsr(:,:)  = 0.0   !! QSR 
     243         zn_dms_qsr(:,:)  = 0.0 
     244         zb_dms_din(:,:)  = 0.0   !! DIN 
     245         zn_dms_din(:,:)  = 0.0 
     246      ENDIF 
     247      !! 
     248      !! calculate stats on these fields 
     249      IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 
     250      fq0 = MINVAL(zn_dms_chn(:,:)) 
     251      fq1 = MAXVAL(zn_dms_chn(:,:)) 
     252      fq2 = SUM(zn_dms_chn(:,:)) 
     253      if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 
     254      fq0 = MINVAL(zn_dms_chd(:,:)) 
     255      fq1 = MAXVAL(zn_dms_chd(:,:)) 
     256      fq2 = SUM(zn_dms_chd(:,:)) 
     257      if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 
     258      fq0 = MINVAL(zn_dms_mld(:,:)) 
     259      fq1 = MAXVAL(zn_dms_mld(:,:)) 
     260      fq2 = SUM(zn_dms_mld(:,:)) 
     261      if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 
     262      fq0 = MINVAL(zn_dms_qsr(:,:)) 
     263      fq1 = MAXVAL(zn_dms_qsr(:,:)) 
     264      fq2 = SUM(zn_dms_qsr(:,:)) 
     265      if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 
     266      fq0 = MINVAL(zn_dms_din(:,:)) 
     267      fq1 = MAXVAL(zn_dms_din(:,:)) 
     268      fq2 = SUM(zn_dms_din(:,:)) 
     269      if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 
     270 
    221271#endif 
    222272  
    … …  
    300350      if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 
    301351         &        fq0, fq1, fq2 
     352      !! 
     353      !! AXY (07/07/15): write out temporally averaged fields for DMS 
     354      !!                 calculations 
     355      !! 
     356      IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS - writing out ...' 
     357      CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_CHN',  zb_dms_chn(:,:)  ) 
     358      CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_CHN',  zn_dms_chn(:,:)  ) 
     359      CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_CHD',  zb_dms_chd(:,:)  ) 
     360      CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_CHD',  zn_dms_chd(:,:)  ) 
     361      CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_MLD',  zb_dms_mld(:,:)  ) 
     362      CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_MLD',  zn_dms_mld(:,:)  ) 
     363      CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_QSR',  zb_dms_qsr(:,:)  ) 
     364      CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_QSR',  zn_dms_qsr(:,:)  ) 
     365      CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_DIN',  zb_dms_din(:,:)  ) 
     366      CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_DIN',  zn_dms_din(:,:)  ) 
     367      !! 
     368      !! calculate stats on these fields 
     369      IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 
     370      fq0 = MINVAL(zn_dms_chn(:,:)) 
     371      fq1 = MAXVAL(zn_dms_chn(:,:)) 
     372      fq2 = SUM(zn_dms_chn(:,:)) 
     373      if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 
     374      fq0 = MINVAL(zn_dms_chd(:,:)) 
     375      fq1 = MAXVAL(zn_dms_chd(:,:)) 
     376      fq2 = SUM(zn_dms_chd(:,:)) 
     377      if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 
     378      fq0 = MINVAL(zn_dms_mld(:,:)) 
     379      fq1 = MAXVAL(zn_dms_mld(:,:)) 
     380      fq2 = SUM(zn_dms_mld(:,:)) 
     381      if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 
     382      fq0 = MINVAL(zn_dms_qsr(:,:)) 
     383      fq1 = MAXVAL(zn_dms_qsr(:,:)) 
     384      fq2 = SUM(zn_dms_qsr(:,:)) 
     385      if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 
     386      fq0 = MINVAL(zn_dms_din(:,:)) 
     387      fq1 = MAXVAL(zn_dms_din(:,:)) 
     388      fq2 = SUM(zn_dms_din(:,:)) 
     389      if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 
     390      !!  
    302391#endif 
    303392 

  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcsms.F90

    r6164 r6201  
    5252      IF( lk_pisces  )   CALL trc_sms_pisces ( kt )    ! main program of PISCES  
    5353      IF( lk_medusa  )   CALL trc_sms_medusa ( kt )    ! MEDUSA  tracers 
     54# if defined key_debug_medusa 
     55         IF(lwp) WRITE(numout,*) '--trcsms : MEDUSA OK --  next IDTRA -- ' 
     56      CALL flush(numout) 
     57# endif 
    5458      IF( lk_idtra   )   CALL trc_sms_idtra  ( kt )    ! radioactive decay of Id. tracer 
     59# if defined key_debug_medusa 
     60         IF(lwp) WRITE(numout,*) '--trcsms : IDTRA OK --  next CFC -- ' 
     61      CALL flush(numout) 
     62# endif 
    5563      IF( lk_cfc     )   CALL trc_sms_cfc    ( kt )    ! surface fluxes of CFC 
     64# if defined key_debug_medusa 
     65         IF(lwp) WRITE(numout,*) '--trcsms : CFC OK --  continue -- ' 
     66      CALL flush(numout) 
     67# endif 
    5668      IF( lk_c14b    )   CALL trc_sms_c14b   ( kt )    ! surface fluxes of C14 
    5769      IF( lk_my_trc  )   CALL trc_sms_my_trc ( kt )    ! MY_TRC  tracers 
Note: See TracChangeset for help on using the changeset viewer.