Changeset 1316 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2009-02-17T16:12:46+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diaptr.F90
r1312 r1316 419 419 CHARACTER (len=40) :: clhstnam, clop ! temporary names 420 420 INTEGER :: iline, it, ji ! 421 REAL(wp) :: zsto, zout, zdt, z max, zjulian ! temporary scalars421 REAL(wp) :: zsto, zout, zdt, zjulian ! temporary scalars 422 422 REAL(wp), DIMENSION(jpj) :: zphi, zfoo 423 423 !!---------------------------------------------------------------------- 424 424 425 ! Define frequency of output and means426 zdt = rdt427 IF( nacc == 1 ) zdt = rdtmin428 IF( ln_mskland ) THEN ; clop = "only(x)" ! put 1.e+20 on land (very expensive!!)429 ELSE ; clop = "x" ! no use of the mask value (require less cpu time)430 ENDIF431 #if defined key_diainstant432 zsto = nf_ptr * zdt433 clop = "inst("//TRIM(clop)//")"434 #else435 zsto=zdt436 clop = "ave("//TRIM(clop)//")"437 #endif438 zout = nf_ptr * zdt439 zmax = ( nitend - nit000 + 1 ) * zdt440 441 442 425 ! define time axis 443 426 it = kt - nit000 + 1 … … 446 429 ! -------------- 447 430 IF( kt == nit000 ) THEN 448 449 zdt = rdt450 IF( nacc == 1 ) zdt = rdtmin431 432 zdt = rdt 433 IF( nacc == 1 ) zdt = rdtmin 451 434 452 435 ! Reference latitude … … 482 465 ! Define frequency of output and means 483 466 zsto = nf_ptr * zdt 484 clop = "ave(x)" 467 IF( ln_mskland ) THEN ; clop = "ave(only(x))" ! put 1.e+20 on land (very expensive!!) 468 ELSE ; clop = "ave(x)" ! no use of the mask value (require less cpu time) 469 ENDIF 485 470 zout = nf_ptr * zdt 486 471 zfoo(:) = 0.e0
Note: See TracChangeset
for help on using the changeset viewer.