- Timestamp:
- 2015-12-15T10:46:14+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_CMCC_merge_2015/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r5930 r6051 42 42 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 43 43 USE in_out_manager ! I/O manager 44 USE diadimg ! dimg direct access file format output45 44 USE iom 46 45 USE ioipsl … … 97 96 END FUNCTION dia_wri_alloc 98 97 99 #if defined key_dimgout100 !!----------------------------------------------------------------------101 !! 'key_dimgout' DIMG output file102 !!----------------------------------------------------------------------103 # include "diawri_dimg.h90"104 105 #else106 98 !!---------------------------------------------------------------------- 107 99 !! Default option NetCDF output file 108 100 !!---------------------------------------------------------------------- 109 # 101 #if defined key_iomput 110 102 !!---------------------------------------------------------------------- 111 103 !! 'key_iomput' use IOM library … … 410 402 INTEGER :: iimi, iima, ipk, it, itmod, ijmi, ijma ! local integers 411 403 INTEGER :: jn, ierror ! local integers 412 REAL(wp) :: zsto, zout, zmax, zjulian , zdt! local scalars404 REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars 413 405 ! 414 406 REAL(wp), POINTER, DIMENSION(:,:) :: zw2d ! 2D workspace … … 435 427 436 428 ! Define frequency of output and means 437 zdt = rdt438 IF( nacc == 1 ) zdt = rdtmin439 429 clop = "x" ! no use of the mask value (require less cpu time and otherwise the model crashes) 440 430 #if defined key_diainstant 441 zsto = nwrite * zdt431 zsto = nwrite * rdt 442 432 clop = "inst("//TRIM(clop)//")" 443 433 #else 444 zsto= zdt434 zsto=rdt 445 435 clop = "ave("//TRIM(clop)//")" 446 436 #endif 447 zout = nwrite * zdt448 zmax = ( nitend - nit000 + 1 ) * zdt437 zout = nwrite * rdt 438 zmax = ( nitend - nit000 + 1 ) * rdt 449 439 450 440 ! Define indices of the horizontal output zoom and vertical limit storage … … 488 478 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 489 479 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 490 & nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set )480 & nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 491 481 CALL histvert( nid_T, "deptht", "Vertical T levels", & ! Vertical grid: gdept 492 482 & "m", ipk, gdept_1d, nz_T, "down" ) … … 524 514 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu 525 515 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 526 & nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set )516 & nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 527 517 CALL histvert( nid_U, "depthu", "Vertical U levels", & ! Vertical grid: gdept 528 518 & "m", ipk, gdept_1d, nz_U, "down" ) … … 537 527 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv 538 528 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 539 & nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set )529 & nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 540 530 CALL histvert( nid_V, "depthv", "Vertical V levels", & ! Vertical grid : gdept 541 531 & "m", ipk, gdept_1d, nz_V, "down" ) … … 550 540 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 551 541 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & 552 & nit000-1, zjulian, zdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set )542 & nit000-1, zjulian, rdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 553 543 CALL histvert( nid_W, "depthw", "Vertical W levels", & ! Vertical grid: gdepw 554 544 & "m", ipk, gdepw_1d, nz_W, "down" ) … … 843 833 ! 844 834 END SUBROUTINE dia_wri 845 # endif846 847 835 #endif 848 836 … … 867 855 INTEGER :: id_i , nz_i, nh_i 868 856 INTEGER, DIMENSION(1) :: idex ! local workspace 869 REAL(wp) :: zsto, zout, zmax, zjulian , zdt857 REAL(wp) :: zsto, zout, zmax, zjulian 870 858 !!---------------------------------------------------------------------- 871 859 ! … … 876 864 clname = cdfile_name 877 865 IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 878 zdt = rdt879 866 zsto = rdt 880 867 clop = "inst(x)" ! no use of the mask value (require less cpu time) 881 868 zout = rdt 882 zmax = ( nitend - nit000 + 1 ) * zdt869 zmax = ( nitend - nit000 + 1 ) * rdt 883 870 884 871 IF(lwp) WRITE(numout,*) … … 895 882 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 896 883 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 gphit884 1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 898 885 CALL histvert( id_i, "deptht", "Vertical T levels", & ! Vertical grid : gdept 899 886 "m", jpk, gdept_1d, nz_i, "down") … … 962 949 ! ----------------- 963 950 CALL histclo( id_i ) 964 #if ! defined key_iomput && ! defined key_dimgout951 #if ! defined key_iomput 965 952 IF( ninist /= 1 ) THEN 966 953 CALL histclo( nid_T )
Note: See TracChangeset
for help on using the changeset viewer.