Changeset 600 for palm/trunk/SOURCE/netcdf.f90

Timestamp:

Nov 24, 2010 4:10:51 PM (13 years ago)

Author:

raasch

Message:

New:
---

Changed:

---

Parameters moved from d3par to inipar: call_psolver_at_all_substeps,
cfl_factor, cycle_mg, mg_cycles,mg_switch_to_pe0_level, ngsrb, nsor,
omega_sor, prandtl_number, psolver, rayleigh_damping_factor,
rayleigh_damping_height, residual_limit (parin, read_var_list, write_var_list)

Due to this change, in routine skip_var_list (end of file read_var_list.f90),
variable ldum is replaced by cdum(LEN=1), because otherwise read errors (too
few data on file) would appear due to one of the additional parameters
(cycle_mg, which is a string of one single character) which are now stored
on the restart file.

Weblink to error message database changed to new trac server (message)

Errors:

---

Bugfix concerning check of cross-section levels on netcdf-files to be
extended (xz,yz) (netcdf)

Bugfix in opening of cross section netcdf-files (parallel opening with
netcdf4 only works for netcdf_data_format > 2) (check_open)

Default values of surface_scalarflux and surface_waterflux changed from 0.0
to 9999999.9. Giving the parameter the default values means, that the
respective surface fluxes are calculated by the MO-relations, so the old default
value did not allow to set the surface fluxes to zero explicitly.

File:

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

palm/trunk/SOURCE/netcdf.f90

@@ -7 +7 @@
 ! Current revisions:
 ! ------------------
-!
+! Bugfix concerning check of cross-section levels on netcdf-files to be
+! extended (xz,yz)
 !
 ! Former revisions:
@@ -2241 +2242 @@
 
 !
-!--       Write yv_xz data 
+!--       Write yv_xz data
           DO  i = 1, ns
              IF( section(i,2) == -1 )  THEN
@@ -2387 +2388 @@
           DO  i = 1, ns
              IF ( section(i,2) /= -1 )  THEN
-                IF ( ( section(i,2) * dy ) /= netcdf_data(i) )  THEN
+                IF ( ( ( section(i,2) + 0.5 ) * dy ) /= netcdf_data(i) )  THEN
                    message_string = 'NetCDF file for cross-sections ' //     &
                                TRIM( var ) // ' from previous run found,' // &
@@ -2971 +2972 @@
           DO  i = 1, ns
              IF ( section(i,3) /= -1 )  THEN
-                IF ( ( section(i,3) * dx ) /= netcdf_data(i) )  THEN
+                IF ( ( ( section(i,3) + 0.5 ) * dx ) /= netcdf_data(i) )  THEN
                    message_string = 'NetCDF file for cross-sections ' //    &
                               TRIM( var ) // ' from previous run found,' // &