- Timestamp:
- 2015-05-18T17:19:50+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90
r4990 r5282 21 21 USE timing ! preformance summary 22 22 USE wrk_nemo ! working arrays 23 USE fldread ! type FLD_N 24 USE phycst ! physical constant 25 USE in_out_manager ! I/O manager 23 26 24 27 IMPLICIT NONE … … 197 200 REAL(wp) :: zztmp 198 201 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: zsaldta ! Jan/Dec levitus salinity 202 ! reading initial file 203 LOGICAL :: ln_tsd_init !: T & S data flag 204 LOGICAL :: ln_tsd_tradmp !: internal damping toward input data flag 205 CHARACTER(len=100) :: cn_dir 206 TYPE(FLD_N) :: sn_tem,sn_sal 207 INTEGER :: ios=0 208 209 NAMELIST/namtsd/ ln_tsd_init,ln_tsd_tradmp,cn_dir,sn_tem,sn_sal 210 ! 211 212 REWIND( numnam_ref ) ! Namelist namtsd in reference namelist : 213 READ ( numnam_ref, namtsd, IOSTAT = ios, ERR = 901) 214 901 IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtsd in reference namelist for dia_ar5', lwp ) 215 REWIND( numnam_cfg ) ! Namelist namtsd in configuration namelist : Parameters of the run 216 READ ( numnam_cfg, namtsd, IOSTAT = ios, ERR = 902 ) 217 902 IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtsd in configuration namelist for dia_ar5', lwp ) 218 IF(lwm) WRITE ( numond, namtsd ) 219 ! 199 220 !!---------------------------------------------------------------------- 200 221 ! … … 216 237 END DO 217 238 IF( lk_mpp ) CALL mpp_sum( vol0 ) 218 219 CALL iom_open ( 'data_1m_salinity_nomask', inum )220 CALL iom_get ( inum, jpdom_data, 'vosaline', zsaldta(:,:,:,1), 1 )221 CALL iom_get ( inum, jpdom_data, 'vosaline', zsaldta(:,:,:,2), 12 )239 240 CALL iom_open ( TRIM( cn_dir )//TRIM(sn_sal%clname), inum ) 241 CALL iom_get ( inum, jpdom_data, TRIM(sn_sal%clvar), zsaldta(:,:,:,1), 1 ) 242 CALL iom_get ( inum, jpdom_data, TRIM(sn_sal%clvar), zsaldta(:,:,:,2), 12 ) 222 243 CALL iom_close( inum ) 223 244 sn0(:,:,:) = 0.5_wp * ( zsaldta(:,:,:,1) + zsaldta(:,:,:,2) )
Note: See TracChangeset
for help on using the changeset viewer.