- Timestamp:
- 2016-01-08T10:35:19+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/SAS_SRC/diawri.F90
- Property svn:keywords set to Id
r4292 r6225 26 26 USE dom_oce ! ocean space and time domain 27 27 USE zdf_oce ! ocean vertical physics 28 USE ldftra_oce ! ocean active tracers: lateral physics29 USE ldfdyn_oce ! ocean dynamics: lateral physics30 USE traldf_iso_grif, ONLY : psix_eiv, psiy_eiv31 USE sol_oce ! solver variables32 28 USE sbc_oce ! Surface boundary condition: ocean fields 33 29 USE sbc_ice ! Surface boundary condition: ice fields … … 40 36 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 41 37 USE in_out_manager ! I/O manager 42 USE diadimg ! dimg direct access file format output43 38 USE diaar5, ONLY : lk_diaar5 44 39 USE iom … … 65 60 66 61 !! * Substitutions 67 # include "zdfddm_substitute.h90"68 # include "domzgr_substitute.h90"69 62 # include "vectopt_loop_substitute.h90" 70 63 !!---------------------------------------------------------------------- 71 64 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 72 !! $Id 65 !! $Id$ 73 66 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 74 67 !!---------------------------------------------------------------------- … … 85 78 END FUNCTION dia_wri_alloc 86 79 87 #if defined key_dimgout88 !!----------------------------------------------------------------------89 !! 'key_dimgout' DIMG output file90 !!----------------------------------------------------------------------91 # include "diawri_dimg.h90"92 93 #else94 80 !!---------------------------------------------------------------------- 95 81 !! Default option NetCDF output file 96 82 !!---------------------------------------------------------------------- 97 # 83 #if defined key_iomput 98 84 !!---------------------------------------------------------------------- 99 85 !! 'key_iomput' use IOM library … … 110 96 !! ** Method : use iom_put 111 97 !!---------------------------------------------------------------------- 112 !! 113 INTEGER, INTENT( in ) :: kt ! ocean time-step index 98 INTEGER, INTENT(in) :: kt ! ocean time-step index 114 99 !!---------------------------------------------------------------------- 115 100 ! … … 144 129 INTEGER :: ierr ! error code return from allocation 145 130 INTEGER :: iimi, iima, ipk, it, itmod, ijmi, ijma ! local integers 146 REAL(wp) :: zsto, zout, zmax, zjulian , zdt! local scalars131 REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars 147 132 !!---------------------------------------------------------------------- 148 133 ! … … 163 148 164 149 ! Define frequency of output and means 165 zdt = rdt166 IF( nacc == 1 ) zdt = rdtmin167 150 IF( ln_mskland ) THEN ; clop = "only(x)" ! put 1.e+20 on land (very expensive!!) 168 151 ELSE ; clop = "x" ! no use of the mask value (require less cpu time) 169 152 ENDIF 170 153 #if defined key_diainstant 171 zsto = nwrite * zdt154 zsto = nwrite * rdt 172 155 clop = "inst("//TRIM(clop)//")" 173 156 #else 174 zsto= zdt157 zsto=rdt 175 158 clop = "ave("//TRIM(clop)//")" 176 159 #endif 177 zout = nwrite * zdt178 zmax = ( nitend - nit000 + 1 ) * zdt160 zout = nwrite * rdt 161 zmax = ( nitend - nit000 + 1 ) * rdt 179 162 180 163 ! Define indices of the horizontal output zoom and vertical limit storage … … 218 201 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 219 202 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 220 & 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 ) 221 204 CALL histvert( nid_T, "deptht", "Vertical T levels", & ! Vertical grid: gdept 222 205 & "m", ipk, gdept_1d, nz_T, "down" ) … … 230 213 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu 231 214 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 232 & 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 ) 233 216 CALL histvert( nid_U, "depthu", "Vertical U levels", & ! Vertical grid: gdept 234 217 & "m", ipk, gdept_1d, nz_U, "down" ) … … 242 225 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv 243 226 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 244 & 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 ) 245 228 CALL histvert( nid_V, "depthv", "Vertical V levels", & ! Vertical grid : gdept 246 229 & "m", ipk, gdept_1d, nz_V, "down" ) … … 338 321 ! 339 322 END SUBROUTINE dia_wri 340 # endif341 342 323 #endif 343 324 … … 362 343 INTEGER :: id_i , nz_i, nh_i 363 344 INTEGER, DIMENSION(1) :: idex ! local workspace 364 REAL(wp) :: zsto, zout, zmax, zjulian , zdt345 REAL(wp) :: zsto, zout, zmax, zjulian 365 346 !!---------------------------------------------------------------------- 366 347 ! … … 373 354 clname = cdfile_name 374 355 IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 375 zdt = rdt376 356 zsto = rdt 377 357 clop = "inst(x)" ! no use of the mask value (require less cpu time) 378 358 zout = rdt 379 zmax = ( nitend - nit000 + 1 ) * zdt359 zmax = ( nitend - nit000 + 1 ) * rdt 380 360 381 361 IF(lwp) WRITE(numout,*) … … 392 372 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 393 373 CALL histbeg( clname, jpi, glamt, jpj, gphit, & 394 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 395 375 CALL histvert( id_i, "deptht", "Vertical T levels", & ! Vertical grid : gdept 396 376 "m", jpk, gdept_1d, nz_i, "down") … … 435 415 ! ----------------- 436 416 CALL histclo( id_i ) 437 #if ! defined key_iomput && ! defined key_dimgout417 #if ! defined key_iomput 438 418 IF( ninist /= 1 ) THEN 439 419 CALL histclo( nid_T )
Note: See TracChangeset
for help on using the changeset viewer.