Ignore:
Timestamp:
Oct 24, 2018 6:39:32 PM (3 years ago)
Author:
gronemeier
Message:

new surface-data output; renamed output variables (pt to theta, rho_air to rho, rho_ocean to rho_sea_water)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/time_integration.f90

    r3418 r3421  
    2525! -----------------
    2626! $Id$
     27! Surface data output
     28!
     29! 3418 2018-10-24 16:07:39Z kanani
    2730! call to material_heat_model now with check if spinup runs (rvtils)
    2831!
     
    432435               skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,        &
    433436               skip_time_do3d, skip_time_domask, skip_time_dopr,               &
    434                skip_time_data_output_av, sloping_surface,                      &
    435                stop_dt, terminate_coupled, terminate_run, timestep_scheme,     &
     437               skip_time_data_output_av, sloping_surface, stop_dt,             &
     438               surface_data_output, terminate_coupled, terminate_run,          &
     439               timestep_scheme,                                                &
    436440               time_coupling, time_do2d_xy, time_do2d_xz, time_do2d_yz,        &
    437441               time_do3d, time_domask, time_dopr, time_dopr_av,                &
     
    519523    USE surface_mod,                                                           &
    520524        ONLY:  surf_def_h, surf_lsm_h, surf_usm_h
     525       
     526    USE surface_output_mod,                                                    &
     527        ONLY:  average_count_surf, averaging_interval_surf, dt_dosurf,         &
     528               dt_dosurf_av, surface_output, surface_output_averaging,         &
     529               skip_time_dosurf, skip_time_dosurf_av, time_dosurf,             &
     530               time_dosurf_av
    521531
    522532    USE turbulence_closure_mod,                                                &
     
    12351245       time_run_control   = time_run_control + dt_3d
    12361246!
     1247!--    Increment time-counter for surface output
     1248       IF ( surface_data_output )  THEN
     1249          IF ( simulated_time >= skip_time_dosurf )  THEN
     1250             time_dosurf    = time_dosurf + dt_3d
     1251          ENDIF
     1252          IF ( simulated_time >= skip_time_dosurf_av )  THEN
     1253             time_dosurf_av = time_dosurf_av + dt_3d
     1254          ENDIF
     1255       ENDIF
     1256!
    12371257!--    In case of synthetic turbulence generation and parametrized turbulence
    12381258!--    information, update the time counter and if required, adjust the
     
    13241344             average_count_3d = average_count_3d + 1
    13251345             time_do_sla = MOD( time_do_sla, MAX( dt_averaging_input, dt_3d ) )
     1346          ENDIF
     1347       ENDIF
     1348!
     1349!--    Average surface data
     1350       IF ( surface_data_output )  THEN
     1351          IF ( averaging_interval_surf /= 0.0_wp  .AND.                        &
     1352               ( dt_dosurf_av - time_dosurf_av ) <= averaging_interval_surf    &
     1353          .AND.  simulated_time >= skip_time_dosurf_av )  THEN
     1354             IF ( time_dosurf_av >= dt_averaging_input )  THEN       
     1355                CALL surface_output_averaging
     1356                average_count_surf = average_count_surf + 1
     1357             ENDIF
    13261358          ENDIF
    13271359       ENDIF
     
    14241456          time_do_av = MOD( time_do_av, MAX( dt_data_output_av, dt_3d ) )
    14251457       ENDIF
     1458!
     1459!--    Output of surface data, instantaneous and averaged data
     1460       IF ( surface_data_output )  THEN
     1461          IF ( time_dosurf >= dt_dosurf )  THEN
     1462             CALL surface_output( 0 )
     1463             time_dosurf = MOD( time_dosurf, MAX( dt_dosurf, dt_3d ) )
     1464          ENDIF
     1465          IF ( time_dosurf_av >= dt_dosurf_av )  THEN
     1466             CALL surface_output( 1 )
     1467             time_dosurf_av = MOD( time_dosurf_av, MAX( dt_dosurf_av, dt_3d ) )
     1468          ENDIF
     1469       ENDIF
    14261470
    14271471!
Note: See TracChangeset for help on using the changeset viewer.