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

Add possibility to output surface data during spinup

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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.