- Timestamp:
- 2018-10-29T15:20:26+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r10249 r10251 44 44 USE in_out_manager ! I/O manager 45 45 USE diadimg ! dimg direct access file format output 46 USE diatmb ! Top,middle,bottom output47 USE dia25h ! 25h Mean output48 USE diaopfoam ! Diaopfoam output49 46 USE iom 50 47 USE ioipsl 51 48 USE dynspg_oce, ONLY: un_adv, vn_adv ! barotropic velocities 52 USE eosbn2 ! equation of state (eos_bn2 routine)53 54 49 55 50 #if defined key_lim2 … … 137 132 REAL(wp), POINTER, DIMENSION(:,:) :: z2d ! 2D workspace 138 133 REAL(wp), POINTER, DIMENSION(:,:,:) :: z3d ! 3D workspace 139 REAL(wp), POINTER, DIMENSION(:,:,:) :: zrhd , zrhop ! 3D workspace140 134 !!---------------------------------------------------------------------- 141 135 ! … … 144 138 CALL wrk_alloc( jpi , jpj , z2d ) 145 139 CALL wrk_alloc( jpi , jpj, jpk , z3d ) 146 CALL wrk_alloc( jpi , jpj, jpk , zrhd , zrhop )147 140 ! 148 141 ! Output the initial state and forcings … … 383 376 CALL iom_put( "v_salttr", 0.5 * z2d ) ! heat transport in j-direction 384 377 ENDIF 385 386 IF( iom_use("rhop") ) THEN387 CALL eos( tsn, zrhd, zrhop, fsdept_n(:,:,:) ) ! now in situ and potential density388 zrhop(:,:,jpk) = 0._wp389 CALL iom_put( 'rhop', zrhop )390 ENDIF391 392 378 ! 393 379 CALL wrk_dealloc( jpi , jpj , z2d ) 394 380 CALL wrk_dealloc( jpi , jpj, jpk , z3d ) 395 CALL wrk_dealloc( jpi , jpj, jpk , zrhd , zrhop )396 !397 ! If we want tmb values398 399 IF (ln_diatmb) THEN400 CALL dia_tmb401 ENDIF402 IF (ln_dia25h) THEN403 CALL dia_25h( kt )404 ENDIF405 IF (ln_diaopfoam) THEN406 CALL dia_diaopfoam407 ENDIF408 381 ! 409 382 IF( nn_timing == 1 ) CALL timing_stop('dia_wri') … … 465 438 zdt = rdt 466 439 IF( nacc == 1 ) zdt = rdtmin 467 clop = "x" ! no use of the mask value (require less cpu time, and otherwise the model crashes) 440 IF( ln_mskland ) THEN ; clop = "only(x)" ! put 1.e+20 on land (very expensive!!) 441 ELSE ; clop = "x" ! no use of the mask value (require less cpu time) 442 ENDIF 468 443 #if defined key_diainstant 469 444 zsto = nwrite * zdt … … 1045 1020 CALL histdef( id_i, "vovvldep", "T point depth" , "m" , & ! t-point depth 1046 1021 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 1047 CALL histdef( id_i, "vovvle3t", "T point thickness" , "m" , & ! t-point depth1048 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout )1049 1022 END IF 1050 1023 … … 1077 1050 CALL histwrite( id_i, "sozotaux", kt, utau , jpi*jpj , idex ) ! i-wind stress 1078 1051 CALL histwrite( id_i, "sometauy", kt, vtau , jpi*jpj , idex ) ! j-wind stress 1079 IF( lk_vvl ) THEN1080 CALL histwrite( id_i, "vovvldep", kt, fsdept_n(:,:,:), jpi*jpj*jpk, idex )! T-cell depth1081 CALL histwrite( id_i, "vovvle3t", kt, fse3t_n (:,:,:), jpi*jpj*jpk, idex )! T-cell thickness1082 END IF1083 1052 1084 1053 ! 3. Close the file
Note: See TracChangeset
for help on using the changeset viewer.