Ignore:
Timestamp:
2016-11-04T08:17:27+01:00 (4 years ago)
Author:
gm
Message:

#1692 - branch SIMPLIF_2_usrdef: e3.=dk[dep.] TOOLS update + domzgr style

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r7188 r7189  
    216216      IF( z_cav == 0._wp ) THEN   ;   ld_isfcav = .false.   ;   ELSE   ;   ld_isfcav = .true.   ;   ENDIF 
    217217      ! 
    218       !                                   ! reference 1D-coordinate 
     218      !                                   ! 1D vertical scale factors (reference coordinate) 
    219219      CALL iom_get( inum, jpdom_unknown, 'e3t_1d'  , pe3t_1d  ) 
    220220      CALL iom_get( inum, jpdom_unknown, 'e3w_1d'  , pe3w_1d  ) 
    221       ! 
    222       IF(  iom_varid( inum, 'gdept_1d', ldstop = .FALSE. ) > 0  .AND.  & 
    223          & iom_varid( inum, 'gdepw_1d', ldstop = .FALSE. ) > 0    ) THEN 
    224          CALL ctl_warn( 'zgr_read : old definition of 1d depths and scale factors used ', &  
    225             &                      'depths at t- and w-points read in domain_cfg file') 
    226          CALL iom_get( inum, jpdom_unknown, 'gdept_1d', pdept_1d )    
    227          CALL iom_get( inum, jpdom_unknown, 'gdepw_1d', pdepw_1d ) 
    228       ELSE  
    229          CALL e3_to_depth( pe3t_1d, pe3w_1d, pdept_1d, pdepw_1d ) 
    230          IF(lwp) THEN                        ! control print 
    231             WRITE(numout,*) 
    232             WRITE(numout,*) '              Reference 1D z-coordinate depth and scale factors:' 
    233             WRITE(numout, "(9x,' level  gdept_1d  gdepw_1d  e3t_1d   e3w_1d  ')" ) 
    234             WRITE(numout, "(10x, i4, 4f9.2)" ) ( jk, pdept_1d(jk), pdepw_1d(jk), pe3t_1d(jk), pe3w_1d(jk), jk = 1, jpk ) 
    235          ENDIF 
    236       ENDIF 
    237221      ! 
    238222      !                                   ! 3D vertical scale factors 
     
    245229      CALL iom_get( inum, jpdom_data, 'e3vw_0' , pe3vw , lrowattr=ln_use_jattr ) 
    246230      ! 
    247       !                                   ! 3D depths 
    248       IF(  iom_varid( inum, 'gdept_0', ldstop = .FALSE. ) > 0  .AND.  & 
    249          & iom_varid( inum, 'gdepw_0', ldstop = .FALSE. ) > 0    ) THEN 
     231      !                                   ! 1D & 3D depths 
     232      ! 
     233      !                                         ! old depth definition (obsolescent feature) 
     234      IF(  iom_varid( inum, 'gdept_1d', ldstop = .FALSE. ) > 0  .AND.  & 
     235         & iom_varid( inum, 'gdepw_1d', ldstop = .FALSE. ) > 0  .AND.  & 
     236         & iom_varid( inum, 'gdept_0' , ldstop = .FALSE. ) > 0  .AND.  & 
     237         & iom_varid( inum, 'gdepw_0' , ldstop = .FALSE. ) > 0    ) THEN 
    250238         CALL ctl_warn( 'zgr_read : old definition of depths and scale factors used ', &  
    251239            &                      'depths at t- and w-points read in domain_cfg file') 
    252          CALL iom_get( inum, jpdom_data, 'gdept_0' , pdept , lrowattr=ln_use_jattr ) 
    253          CALL iom_get( inum, jpdom_data, 'gdepw_0' , pdepw , lrowattr=ln_use_jattr ) 
    254       ELSE 
    255          CALL e3_to_depth( pe3t, pe3w, pdept, pdepw ) 
     240         CALL iom_get( inum, jpdom_unknown, 'gdept_1d', pdept_1d )    
     241         CALL iom_get( inum, jpdom_unknown, 'gdepw_1d', pdepw_1d ) 
     242         CALL iom_get( inum, jpdom_data   , 'gdept_0' , pdept , lrowattr=ln_use_jattr ) 
     243         CALL iom_get( inum, jpdom_data   , 'gdepw_0' , pdepw , lrowattr=ln_use_jattr ) 
     244         ! 
     245      ELSE                                      ! depths computed from e3. scale factors 
     246         ! 
     247         CALL e3_to_depth( pe3t_1d, pe3w_1d, pdept_1d, pdepw_1d )    ! 1D reference depth 
     248         IF(lwp) THEN 
     249            WRITE(numout,*) 
     250            WRITE(numout,*) '              Reference 1D z-coordinate depth and scale factors:' 
     251            WRITE(numout, "(9x,' level  gdept_1d  gdepw_1d  e3t_1d   e3w_1d  ')" ) 
     252            WRITE(numout, "(10x, i4, 4f9.2)" ) ( jk, pdept_1d(jk), pdepw_1d(jk), pe3t_1d(jk), pe3w_1d(jk), jk = 1, jpk ) 
     253         ENDIF 
     254         !                                       
     255         CALL e3_to_depth( pe3t   , pe3w   , pdept   , pdepw    )    ! 3D depths 
     256         ! 
    256257      ENDIF 
    257258      ! 
Note: See TracChangeset for help on using the changeset viewer.