Ignore:
Timestamp:
Mar 19, 2007 8:20:46 AM (17 years ago)
Author:
raasch
Message:

preliminary changes for precipitation output

File:
1 edited

Legend:

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

    r63 r72  
    88! -----------------
    99! New initializing action "by_user" calls user_init_3d_model,
    10 ! ts_value is allocated, +module netcdf_control,
     10! precipitation_amount/rate, ts_value are allocated, +module netcdf_control,
    1111! initial velocities at nzb+1 are regarded for volume
    1212! flow control in case they have been set zero before (to avoid small timesteps)
     
    4040    USE arrays_3d
    4141    USE averaging
     42    USE cloud_parameters
    4243    USE constants
    4344    USE control_parameters
     
    157158!--          Liquid water content
    158159             ALLOCATE ( ql_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     160!
     161!--          Precipitation amount and rate (only needed if output is switched)
     162             ALLOCATE( precipitation_amount(nys-1:nyn+1,nxl-1:nxr+1), &
     163                       precipitation_rate(nys-1:nyn+1,nxl-1:nxr+1) )
    159164          ENDIF
    160165
     
    642647
    643648!
    644 !--    Treating cloud physics, liquid water content is zero at beginning of
    645 !--    the simulation
    646        IF ( cloud_physics )  ql = 0.0
     649!--    Treating cloud physics, liquid water content and precipitation amount
     650!--    are zero at beginning of the simulation
     651       IF ( cloud_physics )  THEN
     652          ql = 0.0
     653          IF ( precipitation )  precipitation_amount = 0.0
     654       ENDIF
    647655
    648656!
Note: See TracChangeset for help on using the changeset viewer.