Changeset 8649


Ignore:
Timestamp:
2017-10-20T16:46:23+02:00 (3 years ago)
Author:
jpalmier
Message:

JPALM — 20-10-17 — ctl_warn text fix

Location:
branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

    r8521 r8649  
    3333      !!---------------------------------------------------------------------- 
    3434      USE bio_medusa_mod 
    35       USE dom_oce,           ONLY: atfp, atfp1, neuler, rdt, e3t_n, tmask 
     35      USE dom_oce,           ONLY: atfp, atfp1, neuler, rdt,            & 
     36                                   e3t_0, e3t_n, tmask 
    3637      USE in_out_manager,    ONLY: lwp, numout 
    3738# if defined key_iomput 
     
    5556      USE trcnam_trp,        ONLY: ln_trcadv_cen2, ln_trcadv_tvd 
    5657  
     58   !!* Substitution 
     59#  include "domzgr_substitute.h90" 
     60 
    5761      !! time (integer timestep) 
    5862      INTEGER, INTENT( in ) ::    kt 
     
    6165      INTEGER :: jn 
    6266 
    63       REAL(wp) :: fq0,fq1,fq2,fq3 
     67      REAL(wp) :: fq0,fq1,fq2,fq3,fq4 
    6468 
    6569# if defined key_roam                      
     
    148152                  fq2 = fq0 + fq1 
    149153                  fq3 = f_benout_n(ji,jj) 
    150                   if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     154                  if (lwp) write (numout,'(a,2i3,a,4f15.5)')                     & 
    151155                     'AXY N   cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',      & 
    152156                     fq0,fq1,fq2,fq3 
     
    162166                  fq2 = fq0 + fq1 
    163167                  fq3 = f_benout_si(ji,jj) 
    164                   if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     168                  if (lwp) write (numout,'(a,2i3,a,4f15.5)')                    & 
    165169                     'AXY Si  cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',     & 
    166170                     fq0,fq1,fq2,fq3 
     
    177181                  fq3 = fq0 + fq1 
    178182                  fq4 = f_benout_c(ji,jj) + f_benout_ca(ji,jj) 
    179                   if (lwp) write (numout,'a,2i3,a,5f15,5)')                   & 
     183                  if (lwp) write (numout,'(a,2i3,a,5f15.5)')                    & 
    180184                     'AXY C   cons: (i,j)=',ji,jj,', (flx,ben,asf,err,out)=', & 
    181185                     fq0,fq1,fq2,fq3,fq4 
     
    191195                   fq2 = fq0 + fq1 
    192196                   fq3 = 2.0 * f_benout_ca(ji,jj) 
    193                    if (lwp) write (numout,'a,2i3,a,4f15,5)')                   & 
     197                   if (lwp) write (numout,'(a,2i3,a,4f15.5)')                    & 
    194198                      'AXY alk cons: (i,j)=',ji,jj,', (flx,ben,err,out)=',     & 
    195199                      fq0,fq1,fq2,fq3 
  • branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r8643 r8649  
    181181      REAL(wp) ::    sumtsn, tsnavg 
    182182      INTEGER  ::    summask 
     183      CHARACTER(25) :: charout, charout2, charout3 
    183184      !! 
    184185      !!------------------------------------------------------------------ 
     
    477478             !! AXY (28/02/14): check input fields 
    478479                  if (ztmp(ji,jj) .lt. -3.0 .or. ztmp(ji,jj) .gt. 40.0 ) then 
     480                     write(charout,*)  tsn(ji,jj,jk,jp_tem) 
     481                     write(charout2,*) mig(ji), mjg(jj), jk, kt  
    479482                     Call ctl_warn(' trc_bio_medusa: T WARNING 3D : ',       & 
    480                                    tsn(ji,jj,jk,jp_tem), 'at I, J, K, kt :', &  
    481                                    mig(ji), mjg(jj), jk, kt ) 
     483                                   TRIM(charout), 'at I, J, K, kt :',        &  
     484                                   TRIM(charout2) ) 
    482485# if defined key_debug_medusa 
    483486                     !! temperature 
     
    543546                          ( tsnavg .GT. 40.0 ) ) THEN     
    544547                        IF (ztmp(ji,jj) .LT. -3.0 ) THEN 
    545                            CALL ctl_warn(                                    & 
    546                            ' trc_bio_medusa: T SWITCHING 3D, ',              & 
    547                            tsn(ji,jj,jk,jp_tem), ' -> -3.0 ' ) 
     548                           write(charout,*)  tsn(ji,jj,jk,jp_tem) 
     549                           Call ctl_warn(' trc_bio_medusa: T SWITCHING 3D : ',     & 
     550                                        TRIM(charout), ' -> -3.0 ') 
    548551                           ztmp(ji,jj) = -3.0 
    549552                        ENDIF 
    550553                        IF (ztmp(ji,jj) .GT. 40.0 ) THEN 
    551                            CALL ctl_warn(                                    & 
    552                            ' trc_bio_medusa: T SWITCHING 3D, ',              & 
    553                            tsn(ji,jj,jk,jp_tem), ' -> 40.0 ' ) 
     554                           write(charout,*)  tsn(ji,jj,jk,jp_tem) 
     555                           Call ctl_warn(' trc_bio_medusa: T SWITCHING 3D : ',     & 
     556                                        TRIM(charout), ' -> 40.0 ') 
    554557                           ztmp(ji,jj) = 40.0 
    555558                        ENDIF 
    556559                     ELSE  
    557                         CALL ctl_warn(                                       & 
    558                         ' trc_bio_medusa: T SWITCHING 3D, ',                 & 
    559                         tsn(ji,jj,jk,jp_tem), ' -> surounding avg : ', tsnavg ) 
     560                        write(charout,*)  tsn(ji,jj,jk,jp_tem) 
     561                        write(charout2,*) tsnavg 
     562                        Call ctl_warn(' trc_bio_medusa: T SWITCHING 3D : ',     & 
     563                                     TRIM(charout), ' -> surounding avg : ',    & 
     564                                     TRIM(charout2) ) 
    560565                        ztmp(ji,jj) = tsnavg 
    561566                     ENDIF 
     
    563568                  !! end T chack 
    564569                  if (zsal(ji,jj) .lt. 1.0 .or. zsal(ji,jj) .gt. 47.0 ) then 
    565                         CALL ctl_warn(                                       & 
    566                         ' trc_bio_medusa: S WARNING 2D, ',                   & 
    567                         tsn(ji,jj,jk,jp_sal),                                & 
    568                         ' at I, J, K, kt : ', mig(ji), mjg(jj), jk, kt) 
     570                     write(charout,*)  tsn(ji,jj,jk,jp_sal) 
     571                     write(charout2,*) mig(ji), mjg(jj), jk, kt  
     572                     Call ctl_warn(' trc_bio_medusa: S WARNING 3D : ',       & 
     573                                   TRIM(charout), 'at I, J, K, kt :',        &  
     574                                   TRIM(charout2) ) 
    569575                     !! Correct out of range values 
    570576                     IF (zsal(ji,jj) .LT. 1.0 ) THEN 
    571                         CALL ctl_warn(                                       & 
    572                         ' trc_bio_medusa: S SWITCHING 3D, ',                 & 
    573                         tsn(ji,jj,jk,jp_sal), ' -> 1.0 ') 
     577                        write(charout,*)  tsn(ji,jj,jk,jp_sal) 
     578                        Call ctl_warn(' trc_bio_medusa: S SWITCHING 3D : ',     & 
     579                                     TRIM(charout), ' -> 1.0 ') 
    574580                        zsal(ji,jj) = 1.0 
    575581                     ENDIF 
    576582                     IF (zsal(ji,jj) .GT. 47.0 ) THEN 
    577                         CALL ctl_warn(                                       & 
    578                         ' trc_bio_medusa: T SWITCHING 3D, ',                 & 
    579                         tsn(ji,jj,jk,jp_sal), ' -> 47.0 ') 
     583                        write(charout,*)  tsn(ji,jj,jk,jp_sal) 
     584                        Call ctl_warn(' trc_bio_medusa: S SWITCHING 3D : ',     & 
     585                                     TRIM(charout), ' -> 47.0 ') 
    580586                        zsal(ji,jj) = 47.0 
    581587                     ENDIF 
Note: See TracChangeset for help on using the changeset viewer.