Changeset 1709 for palm/trunk/SOURCE/read_3d_binary.f90

Timestamp:

Nov 4, 2015 2:47:01 PM (8 years ago)

Author:

maronga

Message:

several bugfixes related to the new surface layer routine and land-surface-radiation interaction

File:

• palm/trunk/SOURCE/read_3d_binary.f90 (modified) (4 diffs)

palm/trunk/SOURCE/read_3d_binary.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added rad_lw_out_change_0, increased binary_version
 ! 
 ! Former revisions:
@@ -124 +124 @@
         ONLY: rad_net, rad_net_av, rad_lw_in, rad_lw_in_av, rad_lw_out,        &
               rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr, rad_lw_hr_av,          &
-              rad_lw_out_av, rad_sw_in, rad_sw_in_av, rad_sw_out,              &
-              rad_sw_out_av, rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr,         &
-              rad_sw_hr_av
+              rad_lw_out_av, rad_lw_out_change_0, rad_sw_in, rad_sw_in_av,     &
+              rad_sw_out, rad_sw_out_av, rad_sw_cs_hr, rad_sw_cs_hr_av,        &
+              rad_sw_hr, rad_sw_hr_av
 
     USE random_function_mod,                                                   &
@@ -348 +348 @@
 !--    First compare the version numbers
        READ ( 13 )  version_on_file
-       binary_version = '4.1'
+       binary_version = '4.2'
        IF ( TRIM( version_on_file ) /= TRIM( binary_version ) )  THEN
           WRITE( message_string, * ) 'version mismatch concerning data ',      &
@@ -877 +877 @@
                    rad_lw_out_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
                              tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
+
+                CASE ( 'rad_lw_out_change_0' )
+                   IF ( .NOT. ALLOCATED( rad_lw_out_change_0 ) )  THEN
+                      ALLOCATE( rad_lw_out_change_0(nysg:nyng,nxlg:nxrg) )
+                   ENDIF
+                   IF ( k == 1 )  READ ( 13 )  tmp_2d
+                   rad_lw_out_change_0(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)&
+                              = tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
 
                 CASE ( 'rad_lw_cs_hr' )