- 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/OPA_SRC/DIA/diawri.F90
r5930 r7351 30 30 USE zdf_oce ! ocean vertical physics 31 31 USE ldftra ! lateral physics: eddy diffusivity coef. 32 USE ldfdyn ! lateral physics: eddy viscosity coef. 32 33 USE sbc_oce ! Surface boundary condition: ocean fields 33 34 USE sbc_ice ! Surface boundary condition: ice fields … … 40 41 USE zdfddm ! vertical physics: double diffusion 41 42 USE diahth ! thermocline diagnostics 43 ! 42 44 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 43 45 USE in_out_manager ! I/O manager 44 USE diadimg ! dimg direct access file format output 46 USE diatmb ! Top,middle,bottom output 47 USE dia25h ! 25h Mean output 45 48 USE iom 46 49 USE ioipsl … … 53 56 USE lib_mpp ! MPP library 54 57 USE timing ! preformance summary 58 USE diurnal_bulk ! diurnal warm layer 59 USE cool_skin ! Cool skin 55 60 USE wrk_nemo ! working array 56 61 … … 74 79 !! * Substitutions 75 80 # include "zdfddm_substitute.h90" 76 # include "domzgr_substitute.h90"77 81 # include "vectopt_loop_substitute.h90" 78 82 !!---------------------------------------------------------------------- … … 97 101 END FUNCTION dia_wri_alloc 98 102 99 #if defined key_dimgout100 !!----------------------------------------------------------------------101 !! 'key_dimgout' DIMG output file102 !!----------------------------------------------------------------------103 # include "diawri_dimg.h90"104 105 #else106 103 !!---------------------------------------------------------------------- 107 104 !! Default option NetCDF output file 108 105 !!---------------------------------------------------------------------- 109 # 106 #if defined key_iomput 110 107 !!---------------------------------------------------------------------- 111 108 !! 'key_iomput' use IOM library … … 143 140 ENDIF 144 141 145 IF( .NOT.lk_vvl ) THEN 146 CALL iom_put( "e3t" , fse3t_n(:,:,:) ) 147 CALL iom_put( "e3u" , fse3u_n(:,:,:) ) 148 CALL iom_put( "e3v" , fse3v_n(:,:,:) ) 149 CALL iom_put( "e3w" , fse3w_n(:,:,:) ) 150 ENDIF 142 ! Output of initial vertical scale factor 143 CALL iom_put("e3t_0", e3t_0(:,:,:) ) 144 CALL iom_put("e3u_0", e3t_0(:,:,:) ) 145 CALL iom_put("e3v_0", e3t_0(:,:,:) ) 146 ! 147 CALL iom_put( "e3t" , e3t_n(:,:,:) ) 148 CALL iom_put( "e3u" , e3u_n(:,:,:) ) 149 CALL iom_put( "e3v" , e3v_n(:,:,:) ) 150 CALL iom_put( "e3w" , e3w_n(:,:,:) ) 151 IF( iom_use("e3tdef") ) & 152 CALL iom_put( "e3tdef" , ( ( e3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 ) 151 153 152 154 CALL iom_put( "ssh" , sshn ) ! sea surface height 153 if( iom_use('ssh2') ) CALL iom_put( "ssh2", sshn(:,:) * sshn(:,:) ) ! square of sea surface height154 155 155 156 CALL iom_put( "toce", tsn(:,:,:,jp_tem) ) ! 3D temperature … … 204 205 CALL iom_put( "sbu", z2d ) ! bottom i-current 205 206 ENDIF 206 207 IF ( ln_dynspg_ts ) THEN208 CALL iom_put( "ubar", un_adv(:,:) ) ! barotropic i-current209 ELSE210 CALL iom_put( "ubar", un_b(:,:) ) ! barotropic i-current211 ENDIF212 207 213 208 CALL iom_put( "voce", vn(:,:,:) ) ! 3D j-current … … 223 218 ENDIF 224 219 225 IF ( ln_dynspg_ts ) THEN226 CALL iom_put( "vbar", vn_adv(:,:) ) ! barotropic j-current227 ELSE228 CALL iom_put( "vbar", vn_b(:,:) ) ! barotropic j-current229 ENDIF230 231 220 CALL iom_put( "woce", wn ) ! vertical velocity 232 221 IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN ! vertical mass transport & its square value … … 243 232 CALL iom_put( "avm" , avmu ) ! T vert. eddy visc. coef. 244 233 CALL iom_put( "avs" , fsavs(:,:,:) ) ! S vert. eddy diff. coef. (useful only with key_zdfddm) 234 235 IF( iom_use('logavt') ) CALL iom_put( "logavt", LOG( MAX( 1.e-20_wp, avt (:,:,:) ) ) ) 236 IF( iom_use('logavs') ) CALL iom_put( "logavs", LOG( MAX( 1.e-20_wp, fsavs(:,:,:) ) ) ) 245 237 246 238 IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN … … 266 258 DO jj = 1, jpj 267 259 DO ji = 1, jpi 268 z2d(ji,jj) = z2d(ji,jj) + fse3t(ji,jj,jk) * tsn(ji,jj,jk,jp_tem) * tmask(ji,jj,jk)260 z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_tem) * tmask(ji,jj,jk) 269 261 END DO 270 262 END DO … … 278 270 DO jj = 1, jpj 279 271 DO ji = 1, jpi 280 z2d(ji,jj) = z2d(ji,jj) + fse3t(ji,jj,jk) * tsn(ji,jj,jk,jp_sal) * tmask(ji,jj,jk)272 z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_sal) * tmask(ji,jj,jk) 281 273 END DO 282 274 END DO … … 290 282 DO jj = 2, jpjm1 291 283 DO ji = fs_2, fs_jpim1 ! vector opt. 292 zztmp = 1._wp / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) )293 zztmpx = 0.5 * ( un(ji-1,jj,jk) * un(ji-1,jj,jk) * e2u(ji-1,jj) * fse3u(ji-1,jj,jk) &294 & + un(ji ,jj,jk) * un(ji ,jj,jk) * e2u(ji ,jj) * fse3u(ji ,jj,jk) ) &284 zztmp = 1._wp / ( e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) 285 zztmpx = 0.5 * ( un(ji-1,jj,jk) * un(ji-1,jj,jk) * e2u(ji-1,jj) * e3u_n(ji-1,jj,jk) & 286 & + un(ji ,jj,jk) * un(ji ,jj,jk) * e2u(ji ,jj) * e3u_n(ji ,jj,jk) ) & 295 287 & * zztmp 296 288 ! 297 zztmpy = 0.5 * ( vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1v(ji,jj-1) * fse3v(ji,jj-1,jk) &298 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1v(ji,jj ) * fse3v(ji,jj ,jk) ) &289 zztmpy = 0.5 * ( vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1v(ji,jj-1) * e3v_n(ji,jj-1,jk) & 290 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1v(ji,jj ) * e3v_n(ji,jj ,jk) ) & 299 291 & * zztmp 300 292 ! … … 307 299 CALL iom_put( "eken", rke ) 308 300 ENDIF 309 301 ! 302 CALL iom_put( "hdiv", hdivn ) ! Horizontal divergence 303 ! 310 304 IF( iom_use("u_masstr") .OR. iom_use("u_heattr") .OR. iom_use("u_salttr") ) THEN 311 305 z3d(:,:,jpk) = 0.e0 312 306 DO jk = 1, jpkm1 313 z3d(:,:,jk) = rau0 * un(:,:,jk) * e2u(:,:) * fse3u(:,:,jk) * umask(:,:,jk)307 z3d(:,:,jk) = rau0 * un(:,:,jk) * e2u(:,:) * e3u_n(:,:,jk) * umask(:,:,jk) 314 308 END DO 315 309 CALL iom_put( "u_masstr", z3d ) ! mass transport in i-direction … … 346 340 z3d(:,:,jpk) = 0.e0 347 341 DO jk = 1, jpkm1 348 z3d(:,:,jk) = rau0 * vn(:,:,jk) * e1v(:,:) * fse3v(:,:,jk) * vmask(:,:,jk)342 z3d(:,:,jk) = rau0 * vn(:,:,jk) * e1v(:,:) * e3v_n(:,:,jk) * vmask(:,:,jk) 349 343 END DO 350 344 CALL iom_put( "v_masstr", z3d ) ! mass transport in j-direction … … 380 374 CALL wrk_dealloc( jpi , jpj, jpk , z3d ) 381 375 ! 376 ! If we want tmb values 377 378 IF (ln_diatmb) THEN 379 CALL dia_tmb 380 ENDIF 381 IF (ln_dia25h) THEN 382 CALL dia_25h( kt ) 383 ENDIF 384 382 385 IF( nn_timing == 1 ) CALL timing_stop('dia_wri') 383 386 ! … … 410 413 INTEGER :: iimi, iima, ipk, it, itmod, ijmi, ijma ! local integers 411 414 INTEGER :: jn, ierror ! local integers 412 REAL(wp) :: zsto, zout, zmax, zjulian , zdt! local scalars415 REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars 413 416 ! 414 417 REAL(wp), POINTER, DIMENSION(:,:) :: zw2d ! 2D workspace … … 418 421 IF( nn_timing == 1 ) CALL timing_start('dia_wri') 419 422 ! 420 CALL wrk_alloc( jpi,jpj , zw2d )421 IF( lk_vvl) CALL wrk_alloc( jpi,jpj,jpk , zw3d )423 CALL wrk_alloc( jpi,jpj , zw2d ) 424 IF( .NOT.ln_linssh ) CALL wrk_alloc( jpi,jpj,jpk , zw3d ) 422 425 ! 423 426 ! Output the initial state and forcings … … 435 438 436 439 ! Define frequency of output and means 437 zdt = rdt438 IF( nacc == 1 ) zdt = rdtmin439 440 clop = "x" ! no use of the mask value (require less cpu time and otherwise the model crashes) 440 441 #if defined key_diainstant 441 zsto = nwrite * zdt442 zsto = nwrite * rdt 442 443 clop = "inst("//TRIM(clop)//")" 443 444 #else 444 zsto= zdt445 zsto=rdt 445 446 clop = "ave("//TRIM(clop)//")" 446 447 #endif 447 zout = nwrite * zdt448 zmax = ( nitend - nit000 + 1 ) * zdt448 zout = nwrite * rdt 449 zmax = ( nitend - nit000 + 1 ) * rdt 449 450 450 451 ! Define indices of the horizontal output zoom and vertical limit storage … … 488 489 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 489 490 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 490 & nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set )491 & nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 491 492 CALL histvert( nid_T, "deptht", "Vertical T levels", & ! Vertical grid: gdept 492 493 & "m", ipk, gdept_1d, nz_T, "down" ) … … 524 525 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu 525 526 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 526 & nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set )527 & nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 527 528 CALL histvert( nid_U, "depthu", "Vertical U levels", & ! Vertical grid: gdept 528 529 & "m", ipk, gdept_1d, nz_U, "down" ) … … 537 538 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv 538 539 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 539 & nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set )540 & nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 540 541 CALL histvert( nid_V, "depthv", "Vertical V levels", & ! Vertical grid : gdept 541 542 & "m", ipk, gdept_1d, nz_V, "down" ) … … 550 551 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 551 552 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 552 & nit000-1, zjulian, zdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set )553 & nit000-1, zjulian, rdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 553 554 CALL histvert( nid_W, "depthw", "Vertical W levels", & ! Vertical grid: gdepw 554 555 & "m", ipk, gdepw_1d, nz_W, "down" ) … … 562 563 CALL histdef( nid_T, "vosaline", "Salinity" , "PSU" , & ! sn 563 564 & jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) 564 IF( lk_vvl) THEN565 IF( .NOT.ln_linssh ) THEN 565 566 CALL histdef( nid_T, "vovvle3t", "Level thickness" , "m" ,& ! e3t_n 566 567 & jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) … … 583 584 CALL histdef( nid_T, "sosfldow", "downward salt flux" , "PSU/m2/s", & ! sfx 584 585 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 585 IF( .NOT. lk_vvl) THEN586 IF( ln_linssh ) THEN 586 587 CALL histdef( nid_T, "sosst_cd", "Concentration/Dilution term on temperature" & ! emp * tsn(:,:,1,jp_tem) 587 588 & , "KgC/m2/s", & ! sosst_cd … … 729 730 ENDIF 730 731 731 IF( lk_vvl) THEN732 CALL histwrite( nid_T, "votemper", it, tsn(:,:,:,jp_tem) * fse3t_n(:,:,:) , ndim_T , ndex_T ) ! heat content733 CALL histwrite( nid_T, "vosaline", it, tsn(:,:,:,jp_sal) * fse3t_n(:,:,:) , ndim_T , ndex_T ) ! salt content734 CALL histwrite( nid_T, "sosstsst", it, tsn(:,:,1,jp_tem) * fse3t_n(:,:,1) , ndim_hT, ndex_hT ) ! sea surface heat content735 CALL histwrite( nid_T, "sosaline", it, tsn(:,:,1,jp_sal) * fse3t_n(:,:,1) , ndim_hT, ndex_hT ) ! sea surface salinity content732 IF( .NOT.ln_linssh ) THEN 733 CALL histwrite( nid_T, "votemper", it, tsn(:,:,:,jp_tem) * e3t_n(:,:,:) , ndim_T , ndex_T ) ! heat content 734 CALL histwrite( nid_T, "vosaline", it, tsn(:,:,:,jp_sal) * e3t_n(:,:,:) , ndim_T , ndex_T ) ! salt content 735 CALL histwrite( nid_T, "sosstsst", it, tsn(:,:,1,jp_tem) * e3t_n(:,:,1) , ndim_hT, ndex_hT ) ! sea surface heat content 736 CALL histwrite( nid_T, "sosaline", it, tsn(:,:,1,jp_sal) * e3t_n(:,:,1) , ndim_hT, ndex_hT ) ! sea surface salinity content 736 737 ELSE 737 738 CALL histwrite( nid_T, "votemper", it, tsn(:,:,:,jp_tem) , ndim_T , ndex_T ) ! temperature … … 740 741 CALL histwrite( nid_T, "sosaline", it, tsn(:,:,1,jp_sal) , ndim_hT, ndex_hT ) ! sea surface salinity 741 742 ENDIF 742 IF( lk_vvl) THEN743 zw3d(:,:,:) = ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2744 CALL histwrite( nid_T, "vovvle3t", it, fse3t_n (:,:,:) , ndim_T , ndex_T ) ! level thickness745 CALL histwrite( nid_T, "vovvldep", it, fsdept_n(:,:,:) , ndim_T , ndex_T ) ! t-point depth743 IF( .NOT.ln_linssh ) THEN 744 zw3d(:,:,:) = ( ( e3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 745 CALL histwrite( nid_T, "vovvle3t", it, e3t_n (:,:,:) , ndim_T , ndex_T ) ! level thickness 746 CALL histwrite( nid_T, "vovvldep", it, gdept_n(:,:,:) , ndim_T , ndex_T ) ! t-point depth 746 747 CALL histwrite( nid_T, "vovvldef", it, zw3d , ndim_T , ndex_T ) ! level thickness deformation 747 748 ENDIF … … 752 753 ! (includes virtual salt flux beneath ice 753 754 ! in linear free surface case) 754 IF( .NOT. lk_vvl) THEN755 IF( ln_linssh ) THEN 755 756 zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_tem) 756 757 CALL histwrite( nid_T, "sosst_cd", it, zw2d, ndim_hT, ndex_hT ) ! c/d term on sst … … 837 838 ENDIF 838 839 ! 839 CALL wrk_dealloc( jpi , jpj , zw2d )840 IF( lk_vvl) CALL wrk_dealloc( jpi , jpj , jpk , zw3d )840 CALL wrk_dealloc( jpi , jpj , zw2d ) 841 IF( .NOT.ln_linssh ) CALL wrk_dealloc( jpi , jpj , jpk , zw3d ) 841 842 ! 842 843 IF( nn_timing == 1 ) CALL timing_stop('dia_wri') 843 844 ! 844 845 END SUBROUTINE dia_wri 845 # endif846 847 846 #endif 848 847 … … 867 866 INTEGER :: id_i , nz_i, nh_i 868 867 INTEGER, DIMENSION(1) :: idex ! local workspace 869 REAL(wp) :: zsto, zout, zmax, zjulian , zdt868 REAL(wp) :: zsto, zout, zmax, zjulian 870 869 !!---------------------------------------------------------------------- 871 870 ! … … 876 875 clname = cdfile_name 877 876 IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 878 zdt = rdt879 877 zsto = rdt 880 878 clop = "inst(x)" ! no use of the mask value (require less cpu time) 881 879 zout = rdt 882 zmax = ( nitend - nit000 + 1 ) * zdt880 zmax = ( nitend - nit000 + 1 ) * rdt 883 881 884 882 IF(lwp) WRITE(numout,*) … … 895 893 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 896 894 CALL histbeg( clname, jpi, glamt, jpj, gphit, & 897 1, jpi, 1, jpj, nit000-1, zjulian, zdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit895 1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 898 896 CALL histvert( id_i, "deptht", "Vertical T levels", & ! Vertical grid : gdept 899 897 "m", jpk, gdept_1d, nz_i, "down") … … 913 911 CALL histdef( id_i, "vovecrtz", "Vertical Velocity" , "m/s" , & ! vertical current 914 912 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 913 ! 914 CALL histdef( id_i, "ahtu" , "u-eddy diffusivity" , "m2/s" , & ! zonal current 915 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 916 CALL histdef( id_i, "ahtv" , "v-eddy diffusivity" , "m2/s" , & ! meridonal current 917 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 918 CALL histdef( id_i, "ahmt" , "t-eddy viscosity" , "m2/s" , & ! zonal current 919 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 920 CALL histdef( id_i, "ahmf" , "f-eddy viscosity" , "m2/s" , & ! meridonal current 921 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 922 ! 915 923 CALL histdef( id_i, "sowaflup", "Net Upward Water Flux" , "Kg/m2/S", & ! net freshwater 916 924 & jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) … … 925 933 CALL histdef( id_i, "sometauy", "Meridional Wind Stress", "N/m2" , & ! j-wind stress 926 934 & jpi, jpj, nh_i, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 927 IF( lk_vvl ) THEN 928 CALL histdef( id_i, "vovvldep", "T point depth" , "m" , & ! t-point depth 935 IF( .NOT.ln_linssh ) THEN 936 CALL histdef( id_i, "vovvldep", "T point depth" , "m" , & ! t-point depth 937 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 938 CALL histdef( id_i, "vovvle3t", "T point thickness" , "m" , & ! t-point depth 929 939 & jpi, jpj, nh_i, jpk, 1, jpk, nz_i, 32, clop, zsto, zout ) 930 940 ENDIF … … 952 962 CALL histwrite( id_i, "vomecrty", kt, vn , jpi*jpj*jpk, idex ) ! now j-velocity 953 963 CALL histwrite( id_i, "vovecrtz", kt, wn , jpi*jpj*jpk, idex ) ! now k-velocity 964 ! 965 CALL histwrite( id_i, "ahtu" , kt, ahtu , jpi*jpj*jpk, idex ) ! aht at u-point 966 CALL histwrite( id_i, "ahtv" , kt, ahtv , jpi*jpj*jpk, idex ) ! - at v-point 967 CALL histwrite( id_i, "ahmt" , kt, ahmt , jpi*jpj*jpk, idex ) ! ahm at t-point 968 CALL histwrite( id_i, "ahmf" , kt, ahmf , jpi*jpj*jpk, idex ) ! - at f-point 969 ! 954 970 CALL histwrite( id_i, "sowaflup", kt, emp-rnf , jpi*jpj , idex ) ! freshwater budget 955 971 CALL histwrite( id_i, "sohefldo", kt, qsr + qns , jpi*jpj , idex ) ! total heat flux … … 959 975 CALL histwrite( id_i, "sometauy", kt, vtau , jpi*jpj , idex ) ! j-wind stress 960 976 977 IF( .NOT.ln_linssh ) THEN 978 CALL histwrite( id_i, "vovvldep", kt, gdept_n(:,:,:), jpi*jpj*jpk, idex )! T-cell depth 979 CALL histwrite( id_i, "vovvle3t", kt, e3t_n (:,:,:) , jpi*jpj*jpk, idex )! T-cell thickness 980 END IF 961 981 ! 3. Close the file 962 982 ! ----------------- 963 983 CALL histclo( id_i ) 964 #if ! defined key_iomput && ! defined key_dimgout984 #if ! defined key_iomput 965 985 IF( ninist /= 1 ) THEN 966 986 CALL histclo( nid_T ) … … 972 992 ! 973 993 END SUBROUTINE dia_wri_state 994 974 995 !!====================================================================== 975 996 END MODULE diawri
Note: See TracChangeset
for help on using the changeset viewer.