Changeset 3485 for palm/trunk/SOURCE/netcdf_data_input_mod.f90

Timestamp:

Nov 3, 2018 5:09:40 PM (6 years ago)

Author:

gronemeier

Message:

Add information about reference point to output files; corrected calculation of longitude values in case crs is not defined in input file

File:

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

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! - get central meridian from origin_lon if crs does not exist
+! - set default origin_lon to 0
+! 
+! 3483 2018-11-02 14:19:26Z raasch
 ! bugfix: misplaced directives for netCDF fixed
 ! 
@@ -362 +366 @@
        REAL(wp) ::  fill_w                  !< fill value for w
        REAL(wp) ::  latitude = 0.0_wp       !< latitude of the lower left corner
-       REAL(wp) ::  longitude = -3.0_wp     !< longitude of the lower left corner
+       REAL(wp) ::  longitude = 0.0_wp      !< longitude of the lower left corner
        REAL(wp) ::  origin_x = 500000.0_wp  !< UTM easting of the lower left corner
        REAL(wp) ::  origin_y = 0.0_wp       !< UTM northing of the lower left corner
@@ -958 +962 @@
                               coord_ref_sys%units,                             &
                               .FALSE., 'crs' )
+       ELSE
+!
+!--       Calculate central meridian from origin_lon
+          coord_ref_sys%longitude_of_central_meridian = &
+             CEILING( input_file_atts%origin_lon / 6.0_wp ) * 6.0_wp - 3.0_wp
        ENDIF
 !