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 321 – NEMO

Changeset 321


Ignore:
Timestamp:
2005-10-03T11:20:35+02:00 (19 years ago)
Author:
opalod
Message:

nemo_v1_bugfix_014:RB: update control print for 1D model

Location:
trunk/NEMO/C1D_SRC
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/C1D_SRC/dyncor1d.F90

    r257 r321  
    1515   USE oce            ! ocean dynamics and tracers 
    1616   USE dom_oce        ! ocean space and time domain 
    17    USE phycst          ! physical constants 
     17   USE phycst         ! physical constants 
    1818   USE in_out_manager ! I/O manager 
     19   USE prtctl         ! Print control 
    1920 
    2021   IMPLICIT NONE 
    … …  
    138139         END DO 
    139140      END DO    
    140       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    141          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    142          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    143          WRITE(numout,*) ' cor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    144          u_ctl = zua   ;   v_ctl = zva 
     141 
     142      IF(ln_ctl)   THEN 
     143         CALL prt_ctl(tab3d_1=ua, clinfo1=' cor  - Ua: ', mask1=umask, & 
     144            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask) 
    145145      ENDIF 
     146 
     147!     IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     148!        zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     149!        zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
     150!        WRITE(numout,*) ' cor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
     151!        u_ctl = zua   ;   v_ctl = zva 
     152!     ENDIF 
     153 
    146154   END SUBROUTINE dyn_cor_1d 
    147155 

  • trunk/NEMO/C1D_SRC/dynnxt1d.F90

    r257 r321  
    1616   USE in_out_manager  ! I/O manager 
    1717   USE lbclnk          ! lateral boundary condition (or mpp link) 
     18   USE prtctl          ! Print control 
    1819 
    1920   IMPLICIT NONE 
    … …  
    128129      !                                                ! =============== 
    129130 
    130       IF(l_ctl)   WRITE(numout,*) ' nxt  - Un: ', SUM(un(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)), & 
    131       &                                  ' Vn: ', SUM(vn(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
     131      IF(ln_ctl)   THEN 
     132         CALL prt_ctl(tab3d_1=un, clinfo1=' nxt_1d  - Un: ', mask1=umask, & 
     133            &         tab3d_2=vn, clinfo2=' Vn: ', mask2=vmask) 
     134      ENDIF 
     135 
     136!     IF(l_ctl)   WRITE(numout,*) ' nxt  - Un: ', SUM(un(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)), & 
     137!     &                                  ' Vn: ', SUM(vn(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
    132138 
    133139   END SUBROUTINE dyn_nxt_1d 

  • trunk/NEMO/C1D_SRC/icestp1d.F90

    r296 r321  
    4040   USE dtasal 
    4141   USE ocfzpt 
    42  
     42   USE prtctl          ! Print control 
    4343 
    4444 
    … …  
    138138          
    139139          
    140          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    141             WRITE(numout,*) 'Ice Forcings ' 
    142             WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (2:nictl,2:njctl) ), ' qsr_ice  : ', SUM( qsr_ice (2:nictl,2:njctl) ) 
    143             WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(2:nictl,2:njctl) ), ' qnsr_ice : ', SUM( qnsr_ice(2:nictl,2:njctl) ) 
    144             WRITE(numout,*) ' evap     : ', SUM( evap    (2:nictl,2:njctl) ) 
    145             WRITE(numout,*) ' precip   : ', SUM( tprecip(2:nictl,2:njctl)  ), ' Snow     : ', SUM( sprecip (2:nictl,2:njctl) ) 
    146             WRITE(numout,*) ' u-stress : ', SUM( gtaux  (2:nictl,2:njctl)  ), ' v-stress : ', SUM( gtauy  (2:nictl,2:njctl) ) 
    147             WRITE(numout,*) ' sst      : ', SUM( sst_io (2:nictl,2:njctl)  ), ' sss      : ', SUM( sss_io (2:nictl,2:njctl) ) 
    148             WRITE(numout,*) ' u_io     : ', SUM( u_io   (2:nictl,2:njctl)  ), ' v_io     : ', SUM( v_io   (2:nictl,2:njctl) ) 
    149             WRITE(numout,*) ' tio_u  1 : ', SUM( tio_u  (2:nictl,2:njctl)  ), ' tio_v    : ', SUM( tio_v  (2:nictl,2:njctl) ) 
    150             WRITE(numout,*) ' hsnif  1 : ', SUM( hsnif  (2:nictl,2:njctl)  ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    151             WRITE(numout,*) ' frld   1 : ', SUM( frld   (2:nictl,2:njctl)  ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
     140         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     141            CALL prt_ctl_info('Ice Forcings ') 
     142            CALL prt_ctl(tab2d_1=qsr_oce ,clinfo1=' qsr_oce  : ', tab2d_2=qsr_ice , clinfo2=' qsr_ice   : ') 
     143            CALL prt_ctl(tab2d_1=qnsr_oce,clinfo1=' qnsr_oce : ', tab2d_2=qnsr_ice, clinfo2=' qnsr_ice  : ') 
     144            CALL prt_ctl(tab2d_1=evap    ,clinfo1=' evap     : ') 
     145            CALL prt_ctl(tab2d_1=tprecip ,clinfo1=' precip   : ', tab2d_2=sprecip , clinfo2=' Snow      : ') 
     146            CALL prt_ctl(tab2d_1=gtaux   ,clinfo1=' u-stress : ', tab2d_2=gtauy   , clinfo2=' v-stress  : ') 
     147            CALL prt_ctl(tab2d_1=sst_io  ,clinfo1=' sst      : ', tab2d_2=sss_io  , clinfo2=' sss       : ') 
     148            CALL prt_ctl(tab2d_1=u_io    ,clinfo1=' u_io     : ', tab2d_2=v_io    , clinfo2=' v_io      : ') 
     149            CALL prt_ctl(tab2d_1=hsnif   ,clinfo1=' hsnif  1 : ', tab2d_2=hicif   , clinfo2=' hicif     : ') 
     150            CALL prt_ctl(tab2d_1=frld    ,clinfo1=' frld   1 : ', tab2d_2=sist    , clinfo2=' sist      : ') 
    152151         ENDIF 
    153152 
    … …  
    167166         CALL lim_thd                                                ! Ice thermodynamics ! 
    168167         !                                                           !--------------------! 
    169          IF(l_ctl) THEN 
    170             WRITE(numout,*) ' hsnif  4 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    171             WRITE(numout,*) ' frld   4 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
    172             WRITE(numout,*) ' u_io   4 : ', SUM( u_io   (2:nictl,2:njctl) ), ' v_io     : ', SUM( v_io   (2:nictl,2:njctl) ) 
    173             WRITE(numout,*) ' tio_u  4 : ', SUM( tio_u  (2:nictl,2:njctl) ), ' tio_v    : ', SUM( tio_v  (2:nictl,2:njctl) ) 
     168         IF(ln_ctl) THEN 
     169            CALL prt_ctl(tab2d_1=hsnif ,clinfo1=' hsnif  2 : ', tab2d_2=hicif , clinfo2=' hicif     : ') 
     170            CALL prt_ctl(tab2d_1=frld  ,clinfo1=' frld   2 : ', tab2d_2=sist  , clinfo2=' sist      : ') 
     171            CALL prt_ctl(tab2d_1=u_io  ,clinfo1=' u_io   4 : ', tab2d_2=v_io  , clinfo2=' v_io      : ') 
     172            CALL prt_ctl(tab2d_1=tio_u  ,clinfo1=' tio_u  4 : ', tab2d_2=tio_v  , clinfo2=' tio_v     : ') 
    174173         ENDIF 
     174 
    175175 
    176176 
    … …  
    180180         !                                                           !------------------------------! 
    181181 
    182          IF(l_ctl) THEN 
    183             WRITE(numout,*) ' hsnif  7 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    184             WRITE(numout,*) ' frld   7 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
    185             WRITE(numout,*) ' tio_u  7 : ', SUM( tio_u  (2:nictl,2:njctl) ), ' tio_v    : ', SUM( tio_v  (2:nictl,2:njctl) ) 
     182         IF(ln_ctl) THEN 
     183            CALL prt_ctl(tab2d_1=hsnif ,clinfo1=' hsnif  7 : ', tab2d_2=hicif , clinfo2=' hicif   : ') 
     184            CALL prt_ctl(tab2d_1=frld  ,clinfo1=' frld   7 : ', tab2d_2=sist  , clinfo2=' sist      : ') 
     185            CALL prt_ctl(tab2d_1=tio_u  ,clinfo1=' tio_u  7 : ', tab2d_2=tio_v  , clinfo2=' tio_v     : ') 
    186186         ENDIF 
    187  
    188187         !                                                           !-------------! 
    189188         CALL lim_wri                                                ! Ice outputs ! 
    … …  
    386385      END DO 
    387386 
    388       IF(l_ctl)   WRITE(numout,*) ' oce_sbc_dmp : a+ = ', aplus, ' a- = ', aminus 
     387      IF(ln_ctl)   WRITE(numout,*) ' oce_sbc_dmp : a+ = ', aplus, ' a- = ', aminus 
    389388#else 
    390389      ! Rigid-lid (emp=emps=E-P-R+Erp) 

  • trunk/NEMO/C1D_SRC/step1d.F90

    r257 r321  
    6969   USE stpctl          ! time stepping control            (stp_ctl routine) 
    7070   USE restart         ! ocean restart                    (rst_wri routine) 
    71   
     71   USE prtctl          ! Print control                    (prt_ctl routine) 
    7272   IMPLICIT NONE 
    7373   PRIVATE 
    … …  
    160160      ENDIF 
    161161 
    162       IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    163          WRITE(numout,*) ' emp  -   : ', SUM( emp   (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
    164          WRITE(numout,*) ' emps -   : ', SUM( emps  (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
    165          WRITE(numout,*) ' qt   -   : ', SUM( qt    (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
    166          WRITE(numout,*) ' qsr  -   : ', SUM( qsr   (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
    167          WRITE(numout,*) ' runoff   : ', SUM( runoff(1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
    168          WRITE(numout,*) ' tmask    : ', SUM( tmask (1:nictl+1,1:njctl+1,:) ) 
    169          WRITE(numout,*) ' sst  -   : ', SUM( tn    (1:nictl+1,1:njctl+1,1) * tmask(1:nictl+1,1:njctl+1,1) ) 
    170          WRITE(numout,*) ' sss  -   : ', SUM( sn    (1:nictl+1,1:njctl+1,1) * tmask(1:nictl+1,1:njctl+1,1) ) 
    171          WRITE(numout,*) ' tau  - x : ', SUM( taux  (1:nictl+1,1:njctl+1) ), ' - y : ', SUM( tauy(1:nictl+1,1:njctl+1) ) 
    172       ENDIF 
     162      IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     163         CALL prt_ctl(tab2d_1=emp    , clinfo1=' emp  -   : ', mask1=tmask, ovlap=1) 
     164         CALL prt_ctl(tab2d_1=emps   , clinfo1=' emps -   : ', mask1=tmask, ovlap=1) 
     165         CALL prt_ctl(tab2d_1=qt     , clinfo1=' qt   -   : ', mask1=tmask, ovlap=1) 
     166         CALL prt_ctl(tab2d_1=qsr    , clinfo1=' qsr  -   : ', mask1=tmask, ovlap=1) 
     167         CALL prt_ctl(tab2d_1=runoff , clinfo1=' runoff   : ', mask1=tmask, ovlap=1) 
     168         CALL prt_ctl(tab3d_1=tmask  , clinfo1=' tmask    : ', mask1=tmask, ovlap=1, kdim=jpk) 
     169         CALL prt_ctl(tab3d_1=tn     , clinfo1=' sst  -   : ', mask1=tmask, ovlap=1, kdim=1) 
     170         CALL prt_ctl(tab3d_1=sn     , clinfo1=' sss  -   : ', mask1=tmask, ovlap=1, kdim=1) 
     171         CALL prt_ctl(tab2d_1=taux   , clinfo1=' tau  - x : ', tab2d_2=tauy, clinfo2='      - y : ', ovlap=1) 
     172      ENDIF 
     173 
    173174 
    174175 
Note: See TracChangeset for help on using the changeset viewer.