- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/SAS_SRC/diawri.F90
r5930 r7351 36 36 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 37 37 USE in_out_manager ! I/O manager 38 USE diadimg ! dimg direct access file format output39 38 USE diaar5, ONLY : lk_diaar5 40 39 USE iom … … 61 60 62 61 !! * Substitutions 63 # include "zdfddm_substitute.h90"64 # include "domzgr_substitute.h90"65 62 # include "vectopt_loop_substitute.h90" 66 63 !!---------------------------------------------------------------------- … … 81 78 END FUNCTION dia_wri_alloc 82 79 83 #if defined key_dimgout84 !!----------------------------------------------------------------------85 !! 'key_dimgout' DIMG output file86 !!----------------------------------------------------------------------87 # include "diawri_dimg.h90"88 89 #else90 80 !!---------------------------------------------------------------------- 91 81 !! Default option NetCDF output file 92 82 !!---------------------------------------------------------------------- 93 # 83 #if defined key_iomput 94 84 !!---------------------------------------------------------------------- 95 85 !! 'key_iomput' use IOM library … … 106 96 !! ** Method : use iom_put 107 97 !!---------------------------------------------------------------------- 108 !! 109 INTEGER, INTENT( in ) :: kt ! ocean time-step index 98 INTEGER, INTENT(in) :: kt ! ocean time-step index 110 99 !!---------------------------------------------------------------------- 111 100 ! … … 140 129 INTEGER :: ierr ! error code return from allocation 141 130 INTEGER :: iimi, iima, ipk, it, itmod, ijmi, ijma ! local integers 142 REAL(wp) :: zsto, zout, zmax, zjulian , zdt! local scalars131 REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars 143 132 !!---------------------------------------------------------------------- 144 133 ! … … 159 148 160 149 ! Define frequency of output and means 161 zdt = rdt162 IF( nacc == 1 ) zdt = rdtmin163 150 IF( ln_mskland ) THEN ; clop = "only(x)" ! put 1.e+20 on land (very expensive!!) 164 151 ELSE ; clop = "x" ! no use of the mask value (require less cpu time) 165 152 ENDIF 166 153 #if defined key_diainstant 167 zsto = nwrite * zdt154 zsto = nwrite * rdt 168 155 clop = "inst("//TRIM(clop)//")" 169 156 #else 170 zsto= zdt157 zsto=rdt 171 158 clop = "ave("//TRIM(clop)//")" 172 159 #endif 173 zout = nwrite * zdt174 zmax = ( nitend - nit000 + 1 ) * zdt160 zout = nwrite * rdt 161 zmax = ( nitend - nit000 + 1 ) * rdt 175 162 176 163 ! Define indices of the horizontal output zoom and vertical limit storage … … 214 201 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 215 202 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 216 & nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set )203 & nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 217 204 CALL histvert( nid_T, "deptht", "Vertical T levels", & ! Vertical grid: gdept 218 205 & "m", ipk, gdept_1d, nz_T, "down" ) … … 226 213 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu 227 214 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 228 & nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set )215 & nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 229 216 CALL histvert( nid_U, "depthu", "Vertical U levels", & ! Vertical grid: gdept 230 217 & "m", ipk, gdept_1d, nz_U, "down" ) … … 238 225 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv 239 226 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 240 & nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set )227 & nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 241 228 CALL histvert( nid_V, "depthv", "Vertical V levels", & ! Vertical grid : gdept 242 229 & "m", ipk, gdept_1d, nz_V, "down" ) … … 334 321 ! 335 322 END SUBROUTINE dia_wri 336 # endif337 338 323 #endif 339 324 … … 358 343 INTEGER :: id_i , nz_i, nh_i 359 344 INTEGER, DIMENSION(1) :: idex ! local workspace 360 REAL(wp) :: zsto, zout, zmax, zjulian , zdt345 REAL(wp) :: zsto, zout, zmax, zjulian 361 346 !!---------------------------------------------------------------------- 362 347 ! … … 369 354 clname = cdfile_name 370 355 IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 371 zdt = rdt372 356 zsto = rdt 373 357 clop = "inst(x)" ! no use of the mask value (require less cpu time) 374 358 zout = rdt 375 zmax = ( nitend - nit000 + 1 ) * zdt359 zmax = ( nitend - nit000 + 1 ) * rdt 376 360 377 361 IF(lwp) WRITE(numout,*) … … 388 372 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 389 373 CALL histbeg( clname, jpi, glamt, jpj, gphit, & 390 1, jpi, 1, jpj, nit000-1, zjulian, zdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit374 1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 391 375 CALL histvert( id_i, "deptht", "Vertical T levels", & ! Vertical grid : gdept 392 376 "m", jpk, gdept_1d, nz_i, "down") … … 431 415 ! ----------------- 432 416 CALL histclo( id_i ) 433 #if ! defined key_iomput && ! defined key_dimgout417 #if ! defined key_iomput 434 418 IF( ninist /= 1 ) THEN 435 419 CALL histclo( nid_T )
Note: See TracChangeset
for help on using the changeset viewer.
