Changeset 600 for palm/trunk/SOURCE/check_open.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/check_open.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_open.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix in opening of cross section netcdf-files (parallel opening with
+! netcdf4 only works for netcdf_data_format > 2)
 !
 ! Former revisions:
@@ -723 +724 @@
 !
 !--          Open an existing NetCDF file for output
+             IF ( netcdf_data_format < 3 )  THEN
+                nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xy(av) )
+             ELSE
 #if defined( __netcdf4 )
-             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xy(av), &
-                                  COMM = comm2d, INFO = MPI_INFO_NULL )
-#else
-             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xy(av) )
-#endif
+                nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xy(av), &
+                                     COMM = comm2d, INFO = MPI_INFO_NULL )
+#endif
+             ENDIF
 
              CALL handle_netcdf_error( 'check_open', 20 )
@@ -817 +820 @@
           IF ( netcdf_extend )  THEN
 !
-!--          Open an existing NetCDF file for output 
+!--          Open an existing NetCDF file for output
+             IF ( netcdf_data_format < 3 )  THEN
+                nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xz(av) )
+             ELSE
 #if defined( __netcdf4 )
-             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xz(av), &
-                                  COMM = comm2d, INFO = MPI_INFO_NULL )
-#else
-             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xz(av) )
-#endif
+                nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xz(av), &
+                                     COMM = comm2d, INFO = MPI_INFO_NULL )
+#endif
+             ENDIF
+
              CALL handle_netcdf_error( 'check_open', 23 )
 !
@@ -911 +917 @@
           IF ( netcdf_extend )  THEN
 !
-!--          Open an existing NetCDF file for output 
+!--          Open an existing NetCDF file for output
+             IF ( netcdf_data_format < 3 )  THEN 
+                nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_yz(av) )
+             ELSE
 #if defined( __netcdf4 )
-             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_yz(av), &
-                                  COMM = comm2d, INFO = MPI_INFO_NULL )
-#else
-             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_yz(av) )
-#endif
+                nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_yz(av), &
+                                     COMM = comm2d, INFO = MPI_INFO_NULL )
+#endif
+             ENDIF
+
              CALL handle_netcdf_error( 'check_open', 26 )
 !