Changeset 4654 for palm


Ignore:
Timestamp:
Aug 28, 2020 1:47:23 PM (3 months ago)
Author:
pavelkrc
Message:

Add possibility to output surface data during spinup

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4653 r4654  
    2828! -----------------
    2929! $Id$
     30! Add possibility to output surface data during spinup (author: J. Resler)
     31!
     32! 4653 2020-08-27 08:54:43Z pavelkrc
    3033! RTM version 4.0:
    3134! - Support full-3D geometry:
     
    3538! - Remove separate direction indices for urban and land faces
    3639! - Bugfix: skyvft output
     40! Authors: P. Krc, J. Resler (Institute of Computer Science, Prague)
    3741!
    3842! 4644 2020-08-22 15:05:46Z pavelkrc
  • palm/trunk/SOURCE/time_integration_spinup.f90

    r4635 r4654  
    116116               spinup_pt_mean,                                                                     &
    117117               spinup_time,                                                                        &
     118               surface_output,                                                                     &
    118119               timestep_count,                                                                     &
    119120               time_dopr,                                                                          &
     
    174175    USE statistics,                                                                                &
    175176        ONLY:  flow_statistics_called
     177
     178    USE surface_data_output_mod,                                                                   &
     179        ONLY:  dt_dosurf,                                                                          &
     180               skip_time_dosurf,                                                                   &
     181               surface_data_output,                                                                &
     182               time_dosurf
    176183
    177184    USE surface_layer_fluxes_mod,                                                                  &
     
    490497             time_dopr         = time_dopr        + dt_3d
    491498          ENDIF
     499          IF ( surface_output )  THEN
     500             IF ( simulated_time >= skip_time_dosurf )  THEN
     501                time_dosurf    = time_dosurf + dt_3d
     502             ENDIF
     503          ENDIF
    492504          time_run_control     = time_run_control + dt_3d
    493505
     
    545557             time_do3d = MOD( time_do3d, MAX( dt_do3d, dt_3d ) )
    546558          ENDIF
    547 
     559!
     560!--       Output of surface data
     561          IF ( surface_output )  THEN
     562             IF ( time_dosurf >= dt_dosurf )  THEN
     563                CALL surface_data_output( 0 )
     564                time_dosurf = MOD( time_dosurf, MAX( dt_dosurf, dt_3d ) )
     565             ENDIF
     566          ENDIF
    548567
    549568       ENDIF
Note: See TracChangeset for help on using the changeset viewer.