Changeset 181 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Jul 30, 2008 7:07:47 AM (16 years ago)

Author:

raasch

Message:

bugfixes + adjustments for SGI ICE system

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -7 +7 @@
 ! Actual revisions:
 ! -----------------
+! Bugfix: zero assignments to tendency arrays in case of restarts
 ! Further amendments to and modifications in the initialisation of the plant 
 ! canopy model 
 ! Allocation of hom_sum moved to parin, initialization of spectrum_x|y directly
 ! after allocating theses arrays,
-! read data for recycling added as new initialization option
+! read data for recycling added as new initialization option,
+! dummy allocation for diss
 !
 ! Former revisions:
@@ -235 +237 @@
     IF ( use_sgs_for_particles )  THEN
        ALLOCATE ( diss(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
+    ELSE
+       ALLOCATE ( diss(2,2,2) )  ! required because diss is used as a
+                                 ! formal parameter
     ENDIF
 
@@ -1046 +1051 @@
        IF ( humidity  .OR.  passive_scalar )  q_p = q
        IF ( ocean )  sa_p = sa
+
+!
+!--    Allthough tendency arrays are set in prognostic_equations, they have
+!--    have to be predefined here because they are used (but multiplied with 0)
+!--    there before they are set. 
+       IF ( timestep_scheme(1:5) == 'runge' )  THEN
+          te_m = 0.0; tpt_m = 0.0; tu_m = 0.0; tv_m = 0.0; tw_m = 0.0
+          IF ( humidity  .OR.  passive_scalar )  tq_m = 0.0
+          IF ( ocean )  tsa_m = 0.0
+       ENDIF
 
     ELSE