Ignore:
Timestamp:
2019-10-29T18:14:49+01:00 (19 months ago)
Author:
laurent
Message:

Update the branch to r11830 of the trunk!

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/ISOMIP/MY_SRC/usrdef_nam.F90

    r10074 r11831  
    4242CONTAINS 
    4343 
    44    SUBROUTINE usr_def_nam( ldtxt, ldnam, cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 
     44   SUBROUTINE usr_def_nam( cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 
    4545      !!---------------------------------------------------------------------- 
    4646      !!                     ***  ROUTINE dom_nam  *** 
     
    5454      !! ** input   : - namusr_def namelist found in namelist_cfg 
    5555      !!---------------------------------------------------------------------- 
    56       CHARACTER(len=*), DIMENSION(:), INTENT(out) ::   ldtxt, ldnam    ! stored print information 
    5756      CHARACTER(len=*)              , INTENT(out) ::   cd_cfg          ! configuration name 
    5857      INTEGER                       , INTENT(out) ::   kk_cfg          ! configuration resolution 
     
    6059      INTEGER                       , INTENT(out) ::   kperio          ! lateral global domain b.c.  
    6160      ! 
    62       INTEGER ::   ios, ii   ! Local integer 
     61      INTEGER ::   ios   ! Local integer 
    6362      !! 
    6463      NAMELIST/namusr_def/ ln_zco, ln_zps, ln_sco, rn_e1deg, rn_e2deg, rn_e3 
    6564      !!---------------------------------------------------------------------- 
    6665      ! 
    67       ii = 1 
    68       ! 
    6966      REWIND( numnam_cfg )          ! Namelist namusr_def (exist in namelist_cfg only) 
    7067      READ  ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 ) 
    71 902   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namusr_def in configuration namelist', .TRUE. ) 
     68902   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namusr_def in configuration namelist' ) 
    7269      ! 
    73       WRITE( ldnam(:), namusr_def ) 
     70      IF(lwm)   WRITE( numond, namusr_def ) 
    7471      ! 
    7572      cd_cfg = 'ISOMIP'           ! name & resolution (not used) 
     
    8178      kpk = INT( rbathy / rn_e3    ) + 1     ! add 1 for t-point in the seafloor 
    8279      ! 
    83       !                             ! control print 
    84       WRITE(ldtxt(ii),*) '   '                                                                              ;   ii = ii + 1 
    85       WRITE(ldtxt(ii),*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg'       ;   ii = ii + 1 
    86       WRITE(ldtxt(ii),*) '~~~~~~~~~~~ '                                                                     ;   ii = ii + 1 
    87       WRITE(ldtxt(ii),*) '   Namelist namusr_def : ISOMIP test case'                                        ;   ii = ii + 1 
    88       WRITE(ldtxt(ii),*) '      type of vertical coordinate : '                                             ;   ii = ii + 1 
    89       WRITE(ldtxt(ii),*) '         z-coordinate flag                     ln_zco   = ', ln_zco               ;   ii = ii + 1 
    90       WRITE(ldtxt(ii),*) '         z-partial-step coordinate flag        ln_zps   = ', ln_zps               ;   ii = ii + 1 
    91       WRITE(ldtxt(ii),*) '         s-coordinate flag                     ln_sco   = ', ln_sco               ;   ii = ii + 1 
    92       WRITE(ldtxt(ii),*) '      resolution'                                                                 ;   ii = ii + 1 
    93       WRITE(ldtxt(ii),*) '         zonal      resolution                 rn_e1deg = ', rn_e1deg, ' degrees' ;   ii = ii + 1 
    94       WRITE(ldtxt(ii),*) '         meridional resolution                 rn_e1deg = ', rn_e1deg, ' degrees' ;   ii = ii + 1 
    95       WRITE(ldtxt(ii),*) '         vertical   resolution                 rn_e3    = ', rn_e3   , ' meters'  ;   ii = ii + 1 
    96       WRITE(ldtxt(ii),*) '      ISOMIP domain = 15° x 10° x 900 m'                                          ;   ii = ii + 1 
    97       WRITE(ldtxt(ii),*) '         resulting global domain size :        jpiglo   = ', kpi                  ;   ii = ii + 1 
    98       WRITE(ldtxt(ii),*) '                                               jpjglo   = ', kpj                  ;   ii = ii + 1 
    99       WRITE(ldtxt(ii),*) '                                               jpkglo   = ', kpk                  ;   ii = ii + 1 
    100       ! 
    101       ! 
    10280      !                             ! Set the lateral boundary condition of the global domain 
    10381      kperio = 0                    ! ISOMIP configuration : close basin 
    10482      ! 
    105       WRITE(ldtxt(ii),*) '   '                                                                          ;   ii = ii + 1 
    106       WRITE(ldtxt(ii),*) '   Lateral boundary condition of the global domain'                           ;   ii = ii + 1 
    107       WRITE(ldtxt(ii),*) '      ISOMIP : closed basin                    jperio   = ', kperio             ;   ii = ii + 1 
     83      !                             ! control print 
     84      IF(lwp) THEN 
     85         WRITE(numout,*) '   ' 
     86         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     87         WRITE(numout,*) '~~~~~~~~~~~ ' 
     88         WRITE(numout,*) '   Namelist namusr_def : ISOMIP test case' 
     89         WRITE(numout,*) '      type of vertical coordinate : ' 
     90         WRITE(numout,*) '         z-coordinate flag                     ln_zco   = ', ln_zco 
     91         WRITE(numout,*) '         z-partial-step coordinate flag        ln_zps   = ', ln_zps 
     92         WRITE(numout,*) '         s-coordinate flag                     ln_sco   = ', ln_sco 
     93         WRITE(numout,*) '      resolution' 
     94         WRITE(numout,*) '         zonal      resolution                 rn_e1deg = ', rn_e1deg, ' degrees' 
     95         WRITE(numout,*) '         meridional resolution                 rn_e1deg = ', rn_e1deg, ' degrees' 
     96         WRITE(numout,*) '         vertical   resolution                 rn_e3    = ', rn_e3   , ' meters' 
     97         WRITE(numout,*) '      ISOMIP domain = 15° x 10° x 900 m' 
     98         WRITE(numout,*) '         resulting global domain size :        jpiglo   = ', kpi 
     99         WRITE(numout,*) '                                               jpjglo   = ', kpj 
     100         WRITE(numout,*) '                                               jpkglo   = ', kpk 
     101         WRITE(numout,*) '   ' 
     102         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     103         WRITE(numout,*) '      ISOMIP : closed basin                    jperio   = ', kperio 
     104      ENDIF 
    108105      ! 
    109106   END SUBROUTINE usr_def_nam 
Note: See TracChangeset for help on using the changeset viewer.