Changeset 11358 for NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DIA/diawri.F90
- Timestamp:
- 2019-07-26T16:10:19+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DIA/diawri.F90
r11325 r11358 426 426 !! define all the NETCDF files and fields 427 427 !! At each time step call histdef to compute the mean if ncessary 428 !! Each n write time step, output the instantaneous or mean fields428 !! Each nn_write time step, output the instantaneous or mean fields 429 429 !!---------------------------------------------------------------------- 430 430 INTEGER, INTENT( in ) :: kt ! ocean time-step index … … 442 442 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d ! 3D workspace 443 443 !!---------------------------------------------------------------------- 444 !445 IF( ln_timing ) CALL timing_start('dia_wri')446 444 ! 447 445 IF( ninist == 1 ) THEN !== Output the initial state and forcings ==! … … 450 448 ENDIF 451 449 ! 450 IF( nn_write == -1 ) RETURN ! we will never do any output 451 ! 452 IF( ln_timing ) CALL timing_start('dia_wri') 453 ! 452 454 ! 0. Initialisation 453 455 ! ----------------- … … 459 461 clop = "x" ! no use of the mask value (require less cpu time and otherwise the model crashes) 460 462 #if defined key_diainstant 461 zsto = n write * rdt463 zsto = nn_write * rdt 462 464 clop = "inst("//TRIM(clop)//")" 463 465 #else … … 465 467 clop = "ave("//TRIM(clop)//")" 466 468 #endif 467 zout = n write * rdt469 zout = nn_write * rdt 468 470 zmax = ( nitend - nit000 + 1 ) * rdt 469 471 … … 496 498 ! WRITE root name in date.file for use by postpro 497 499 IF(lwp) THEN 498 CALL dia_nam( clhstnam, n write,' ' )500 CALL dia_nam( clhstnam, nn_write,' ' ) 499 501 CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 500 502 WRITE(inum,*) clhstnam … … 504 506 ! Define the T grid FILE ( nid_T ) 505 507 506 CALL dia_nam( clhstnam, n write, 'grid_T' )508 CALL dia_nam( clhstnam, nn_write, 'grid_T' ) 507 509 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 508 510 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit … … 540 542 ! Define the U grid FILE ( nid_U ) 541 543 542 CALL dia_nam( clhstnam, n write, 'grid_U' )544 CALL dia_nam( clhstnam, nn_write, 'grid_U' ) 543 545 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 544 546 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu … … 553 555 ! Define the V grid FILE ( nid_V ) 554 556 555 CALL dia_nam( clhstnam, n write, 'grid_V' ) ! filename557 CALL dia_nam( clhstnam, nn_write, 'grid_V' ) ! filename 556 558 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 557 559 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv … … 566 568 ! Define the W grid FILE ( nid_W ) 567 569 568 CALL dia_nam( clhstnam, n write, 'grid_W' ) ! filename570 CALL dia_nam( clhstnam, nn_write, 'grid_W' ) ! filename 569 571 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 570 572 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit … … 741 743 ! donne le nombre d'elements, et ndex la liste des indices a sortir 742 744 743 IF( lwp .AND. MOD( itmod, n write ) == 0 ) THEN745 IF( lwp .AND. MOD( itmod, nn_write ) == 0 ) THEN 744 746 WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step' 745 747 WRITE(numout,*) '~~~~~~ '
Note: See TracChangeset
for help on using the changeset viewer.