Changeset 341 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys
- Timestamp:
- Nov 3, 2011, 2:28:57 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
r324 r341 412 412 ENDIF 413 413 414 !!!! for 'subgrid' temporal diagnostics415 test2 = MODULO(elaps,history_interval*100.)414 !!!!! for 'subgrid' temporal diagnostics 415 !test2 = MODULO(elaps,history_interval*100.) 416 416 417 417 !!!******!! … … 478 478 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 479 479 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 480 IF (test2.EQ.0) THEN481 print *, 'compute stats'482 print *, 'RESET'483 uave = uave*0.484 vave = vave*0.485 tave = tave*0.486 wave = wave*0.487 ustd = ustd*0.488 vstd = vstd*0.489 tstd = tstd*0.490 wstd = wstd*0.491 ENDIF492 uave = uave + u * dt / (float(history_interval)*100.)493 vave = vave + v * dt / (float(history_interval)*100.)494 tave = tave + th * dt / (float(history_interval)*100.)495 wave = wave + w * dt / (float(history_interval)*100.)496 ustd = ustd + u * u * dt / (float(history_interval)*100.)497 vstd = vstd + v * v * dt / (float(history_interval)*100.)498 tstd = tstd + th * th * dt / (float(history_interval)*100.)499 wstd = wstd + w * w * dt / (float(history_interval)*100.)480 !IF (test2.EQ.0) THEN 481 ! print *, 'compute stats' 482 ! print *, 'RESET' 483 ! uave = uave*0. 484 ! vave = vave*0. 485 ! tave = tave*0. 486 ! wave = wave*0. 487 ! ustd = ustd*0. 488 ! vstd = vstd*0. 489 ! tstd = tstd*0. 490 ! wstd = wstd*0. 491 !ENDIF 492 ! uave = uave + u * dt / (float(history_interval)*100.) 493 ! vave = vave + v * dt / (float(history_interval)*100.) 494 ! tave = tave + th * dt / (float(history_interval)*100.) 495 ! wave = wave + w * dt / (float(history_interval)*100.) 496 ! ustd = ustd + u * u * dt / (float(history_interval)*100.) 497 ! vstd = vstd + v * v * dt / (float(history_interval)*100.) 498 ! tstd = tstd + th * th * dt / (float(history_interval)*100.) 499 ! wstd = wstd + w * w * dt / (float(history_interval)*100.) 500 500 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 501 501 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Note: See TracChangeset
for help on using the changeset viewer.