Changeset 1701 for palm/trunk/SOURCE/radiation_model.f90

Timestamp:

Nov 2, 2015 7:43:04 AM (8 years ago)

Author:

maronga

Message:

minor bugfixes for radiation model. bugfix in subjob

File:

• palm/trunk/SOURCE/radiation_model.f90 (modified) (7 diffs)

palm/trunk/SOURCE/radiation_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfixes: wrong index for output of timeseries, setting of nz_snd_end
 ! 
 ! Former revisions:
@@ -628 +628 @@
        dots_num = dots_num + 5
 
-       dots_label(dots_rad+1) = "rad_net"
-       dots_label(dots_rad+2) = "rad_lw_in"
-       dots_label(dots_rad+3) = "rad_lw_out"
-       dots_label(dots_rad+4) = "rad_sw_in"
-       dots_label(dots_rad+5) = "rad_sw_out"
+       dots_label(dots_rad) = "rad_net"
+       dots_label(dots_rad+1) = "rad_lw_in"
+       dots_label(dots_rad+2) = "rad_lw_out"
+       dots_label(dots_rad+3) = "rad_sw_in"
+       dots_label(dots_rad+4) = "rad_sw_out"
        dots_unit(dots_rad:dots_rad+4) = "W/m2"
 
@@ -1161 +1161 @@
 !
 ! !--    Allocate temporary array for storing pressure data
-       ALLOCATE( hyp_snd_tmp(nzb+1:nz_snd) )
+       ALLOCATE( hyp_snd_tmp(1:nz_snd) )
        hyp_snd_tmp = 0.0_wp
 
@@ -1173 +1173 @@
 !
 !--    Allocate temporary array for storing temperature data
-       ALLOCATE( t_snd_tmp(nzb+1:nz_snd) )
+       ALLOCATE( t_snd_tmp(1:nz_snd) )
        t_snd_tmp = 0.0_wp
 
@@ -1186 +1186 @@
 !--    Calculate start of sounding data
        nz_snd_start = nz_snd + 1
-       nz_snd_end   = nz_snd_end
+       nz_snd_end   = nz_snd + 1
 
 !
@@ -1199 +1199 @@
 
        IF ( nz_snd_start <= nz_snd )  THEN
-          nz_snd_end = nz_snd - 1
+          nz_snd_end = nz_snd
        END IF
 
@@ -1205 +1205 @@
 !
 !--    Calculate of total grid points for RRTMG calculations
-       nzt_rad = nzt + nz_snd_end - nz_snd_start + 2
+       nzt_rad = nzt + nz_snd_end - nz_snd_start + 1
 
 !