AirTemperatureDailyMeanInit Subroutine
public subroutine AirTemperatureDailyMeanInit(ini, mask, dtMeteo, tstart, dem_loaded)
Initialize air temperature
Arguments
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(IniList), | intent(in) | :: | ini | |||
| type(grid_integer), | intent(in) | :: | mask |
defines interpolation extent |
||
| integer(kind=short), | intent(in) | :: | dtMeteo |
deltat of meteo data reading |
||
| type(DateTime), | intent(in) | :: | tstart |
initial time |
||
| logical, | intent(in) | :: | dem_loaded |
true if dem has been loaded |
Variables
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| character(len=1000), | public | :: | filename | ||||
| type(DateTime), | public | :: | lapseTime | ||||
| type(grid_real), | public | :: | meteoTemp |
Source Code
SUBROUTINE AirTemperatureDailyMeanInit & ! ( ini, mask, dtMeteo, tstart, dem_loaded ) IMPLICIT NONE TYPE (IniList), INTENT(IN) :: ini TYPE (grid_integer), INTENT(IN) :: mask !!defines interpolation extent INTEGER (KIND = short), INTENT(IN) :: dtMeteo !! deltat of meteo data reading TYPE (DateTime), INTENT(IN) :: tstart !!initial time LOGICAL , INTENT (IN) :: dem_loaded !! true if dem has been loaded !local declarations CHARACTER (LEN = 1000) :: filename TYPE (DateTime) :: lapseTime TYPE (grid_real) :: meteoTemp !-------------------------end of declarations---------------------------------- !set initial time timeNew = tstart !initialize grid CALL NewGrid (temperatureAirDailyMean, mask, 0.) CALL NewGrid (grid_devst, mask, 0.) !set deltat dtTemperatureDailyMean = day !check dt is multiple of dtMeteo IF (.NOT.(MOD(dtTemperatureDailyMean,dtMeteo) == 0)) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'dt is not multiple of dtMeteo') END IF !set valid threshold IF (KeyIsPresent ('valid-threshold', ini, & section = 'temperature-daily-mean') ) THEN valid_prcn = IniReadReal ('valid-threshold', ini, & section = 'temperature-daily-mean') ELSE ! set to default = 1.0 CALL Catch ('info', 'AirTemperatureDailyMean', & 'valid-threshold not defined, set to 1') valid_prcn = 1. END IF !set cell-size cellsizeInterpolation = mask % cellsize !interpolation-assignment method IF (KeyIsPresent('interpolation-assignment', ini, & section = 'temperature-daily-mean') ) THEN interpolationMethod_assignment = IniReadInt ('interpolation-assignment', ini, & section = 'temperature-daily-mean') ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'interpolation-assignment missing in meteo configuration file') END IF !set interpolation method IF (interpolationMethod_assignment == 1) THEN interpolationMethod = IniReadInt ('interpolation', ini, & section = 'temperature-daily-mean') CALL SetSpecificProperties (interpolationMethod, ini, mask) ELSE !read map interpolationMethod = -1 CALL GridByIni (ini, interpolationMethod_map, 'temperature-daily-mean', & "interpolation") !scan for interpolation methods included interpolationMethod_vector = 0 DO i = 1, interpolationMethod_map % idim DO j = 1, interpolationMethod_map % jdim IF (interpolationMethod_map % mat(i,j) /= & interpolationMethod_map % nodata) THEN interpolationMethod_vector & (interpolationMethod_map % mat (i,j)) = & interpolationMethod_map % mat (i,j) END IF END DO END DO DO i = 1, 3 CALL SetSpecificProperties (interpolationMethod_vector (i), ini, mask) END DO END IF !scale factor and offset IF (KeyIsPresent('offset', ini, section = 'temperature-daily-mean')) THEN offset_value = IniReadReal ('offset', ini, section = 'temperature-daily-mean') ELSE offset_value = MISSING_DEF_REAL END IF IF (KeyIsPresent('scale-factor', ini, section = 'temperature-daily-mean')) THEN scale_factor = IniReadReal ('scale-factor', ini, section = 'temperature-daily-mean') ELSE scale_factor = MISSING_DEF_REAL END IF !set power_idw IF (KeyIsPresent('idw-power', ini, section = 'temperature-daily-mean')) THEN idw_power = IniReadReal ('idw-power', ini, section = 'temperature-daily-mean') ELSE !set default value idw_power = 2. END IF !file filename = IniReadString ('file', ini, section = 'temperature-daily-mean') IF (interpolationMethod_assignment == 1 .AND. & interpolationMethod == 0) THEN !data are stored in net-cdf file !store net-cdf filename fileGrid = filename IF ( KeyIsPresent ('variable', ini, section = 'temperature-daily-mean') ) THEN temperatureAirDailyMean % var_name = IniReadString ('variable', & ini, section = 'temperature-daily-mean') !read grid and store as temporary file CALL NewGrid (layer = meteoTemp, fileName = TRIM(fileGrid), & fileFormat = NET_CDF, & variable = TRIM(temperatureAirDailyMean % var_name) ) ELSE IF (KeyIsPresent ('standard_name', ini, & section = 'temperature-daily-mean') ) THEN temperatureAirDailyMean % standard_name = IniReadString ('standard_name', & ini, section = 'temperature-daily-mean') ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'standard_name or variable missing in section temeprature' ) END IF !set cellsize to zero. Optimal cellsize is automatically computed cellsizeInterpolation = 0. ELSE !open file containing local measurements fileunit = GetUnit () OPEN(fileunit,file = filename(1:LEN_TRIM(filename)),status='old') END IF !use elevation to drift interpolation IF ( KeyIsPresent ('elevation-drift', ini, section = 'temperature-daily-mean') ) THEN elevationDrift = IniReadInt ('elevation-drift', ini, section = 'temperature-daily-mean') ELSE elevationDrift = 0 !default, suppress drift END IF IF (elevationDrift == 1) THEN IF (interpolationMethod == 0 ) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'elevation drift cannot be applied when interpolation = 0 ') END IF !check if have been loaded by domain properties IF ( .NOT. dem_loaded) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'dem for elevation drift was not loaded ') END IF !need lapse-rate computation? IF (KeyIsPresent ('lapse-rate-computation', ini, & section = 'temperature-daily-mean') ) THEN lapse_rate_computation = IniReadInt ('lapse-rate-computation', ini, & section = 'temperature-daily-mean') ELSE !lapse-rate-computation is missing CALL Catch ('error', 'AirTemperatureDailyMean', & 'lapse-rate-computation missing in meteo configuration file') END IF IF ( lapse_rate_computation == 1 .AND. interpolationMethod == 0 ) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'lapse_rate_computation option requires interpolation different from 0') END IF IF (KeyIsPresent ('r2min', ini, & section = 'temperature-daily-mean') ) THEN r2min = IniReadInt ('r2min', ini, & section = 'temperature-daily-mean') ELSE !r2min is missing IF (lapse_rate_computation == 1) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'r2min missing in meteo configuration file') END IF END IF !set lapse rate assignment IF (KeyIsPresent ('lapse-rate-assignment', ini, & section = 'temperature-daily-mean') ) THEN lapse_rate_assignment = IniReadInt ('lapse-rate-assignment', ini, & section = 'temperature-daily-mean') ELSE !lapse-rate-assignment missing IF (lapse_rate_computation == 0) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'lapse-rate-assignment missing in meteo configuration file') ELSE lapse_rate_assignment = 1 END IF END IF !lapse rate map IF (lapse_rate_assignment == 2) THEN IF ( SubSectionIsPresent (subsection = 'lapse-rate-map', & section = 'temperature-daily-mean', iniDb = ini) ) THEN !check all keywords are defined IF (IniReadString ('format', ini, 'temperature-daily-mean',& 'lapse-rate-map') == 'net-cdf') THEN IF (KeyIsPresent ('time',ini, 'temperature-daily-mean', & 'lapse-rate-map')) THEN timeString = IniReadString ('time', ini, & 'temperature-daily-mean', 'lapse-rate-map') lapseTime = timeString IF (lapseTime > tstart) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'lapse-rate-map time greater than starting time' ) END IF ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'time keyword is missing in lapse-rate-map subsection' ) END IF END IF CALL GridByIni (ini, lapse_map, 'temperature-daily-mean', "lapse-rate-map") !check coordinate reference system and spatial resolution IF (.NOT. CRSisEqual (mask = mask, grid = lapse_map, & checkCells= .TRUE.)) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'wrong lapse-rate-map spatial reference system ') END IF ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'lapse-rate-map missing in meteo configuration file') END IF ELSE !build lapse map from scalar value IF ( KeyIsPresent ('lapse-rate-scalar', ini, & section = 'temperature-daily-mean') ) THEN lapse = IniReadReal ('lapse-rate-scalar', ini, & section = 'temperature-daily-mean') CALL NewGrid (lapse_map, mask) lapse_map = lapse ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'lapse-rate-scalar missing in meteo configuration file') END IF END IF END IF !grid exporting settings IF (KeyIsPresent ('export', ini, section = 'temperature-daily-mean') ) THEN export = IniReadInt ('export', ini, section = 'temperature-daily-mean') ELSE export = 0 END IF IF (export == 1) THEN !set export path name IF (KeyIsPresent ('export-path', ini, section = 'temperature-daily-mean') ) THEN export_path = IniReadString ('export-path', ini, section = 'temperature-daily-mean') !check if path ends with / or \ IF ( export_path (LEN_TRIM (export_path):LEN_TRIM (export_path)) /= '\' .AND. & export_path (LEN_TRIM (export_path):LEN_TRIM (export_path)) /= '/' ) THEN export_path (LEN_TRIM (export_path)+1:LEN_TRIM (export_path)+1) = '\' END IF IF (INDEX (export_path,'.') == 1) THEN !detected relative path, remove '.' export_path = export_path (2:LEN_TRIM(export_path)) !build absolute path export_path = TRIM(CurrentDir() ) // TRIM(export_path) END IF !check OS convention IF (GetOS () == WIN32) THEN export_path = ReplaceChar (export_path,'/','\') ELSE export_path = ReplaceChar (export_path,'\','/') END IF !check folder exists IF ( .NOT. DirExists (TRIM (export_path) ) ) THEN CALL Catch ('info', 'AirTemperatureDailyMean', & 'creating directory: ', argument = TRIM(export_path)) CALL DirNew (export_path) END IF ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'missing export-path') END IF !starting time IF (KeyIsPresent ('export-start', ini, section = 'temperature-daily-mean') ) THEN timeString = IniReadString ('export-start', ini, 'temperature-daily-mean') export_start = timeString ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'missing export-start') END IF !initialize timeNewExport timeNewExport = export_start !ending time IF (KeyIsPresent ('export-stop', ini, section = 'temperature-daily-mean') ) THEN timeString = IniReadString ('export-stop', ini, 'temperature-daily-mean') export_stop = timeString ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'missing export-start') END IF !export dt IF (KeyIsPresent ('export-dt', ini, section = 'temperature-daily-mean') ) THEN export_dt = IniReadInt ('export-dt', ini, section = 'temperature-daily-mean') ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'missing export-dt') END IF !coordinate reference system IF (KeyIsPresent ('export-epsg', ini, section = 'temperature-daily-mean') ) THEN export_epsg = IniReadInt ('export-epsg', ini, section = 'temperature-daily-mean') exportGridMapping = DecodeEPSG (export_epsg) IF (exportGridMapping == temperatureAirDailyMean % grid_mapping) THEN needConversion = .FALSE. !initialize grid for exporting with CRS as meteo variable CALL NewGrid (exportedGrid, temperatureAirDailyMean) CALL NewGrid (exportedGridVar, temperatureAirDailyMean) ELSE needConversion = .TRUE. !initialize grid for converting with required CRS exportedGrid % grid_mapping = DecodeEPSG (export_epsg) exportedGridVar % grid_mapping = DecodeEPSG (export_epsg) END IF ELSE CALL Catch ('info', 'AirTemperatureDailyMean', & 'export-epsg not defined, use default') needConversion = .FALSE. !initialize grid for exporting with CRS CALL NewGrid (exportedGrid, temperatureAirDailyMean) CALL NewGrid (exportedGridVar, temperatureAirDailyMean) END IF !export file format IF (KeyIsPresent ('export-format', ini, section = 'temperature-daily-mean') ) THEN export_format = IniReadInt ('export-format', ini, section = 'temperature-daily-mean') ELSE CALL Catch ('error', 'AirTemperatureDailyMean', & 'missing export-format') END IF IF (export_format == 3) THEN !grid map CALL SetLongName ( 'air_temperature', exportedGrid) CALL SetStandardName ( 'air_temperature', exportedGrid) CALL SetUnits ('degree_Celsius', exportedGrid) !this unit is for exporting grid, it is converted internally to m/s !if file exists, remove it export_file = TRIM(export_path) // 'temperature_daily_mean.nc' IF ( FileExists (export_file) ) THEN CALL FileDelete (export_file) END IF !variance of kriging CALL SetLongName ( 'air_temperature_variance', exportedGridVar) CALL SetStandardName ( 'air_temperature_variance', exportedGridVar) CALL SetUnits ('degree_Celsius', exportedGridVar) !this unit is for exporting grid !if file exists, remove it export_file_var = TRIM(export_path) // 'temperature_daily_mean_variance.nc' IF ( FileExists (export_file_var) ) THEN CALL FileDelete (export_file_var) END IF END IF END IF !complete initialization IF (interpolationMethod == 0) THEN !Get the dt of imported field. Assume dt is regular dtGrid = GetDtGrid (filename = fileGrid, checkRegular = .TRUE.) !check dt is multiple of dtGrid IF (.NOT.(MOD(dtTemperatureDailyMean,dtGrid) == 0)) THEN CALL Catch ('error', 'AirTemperatureDailyMean', & 'dt temperature daily mean is not multiple of dt of input grid') END IF ELSE !populate metadata CALL ReadMetadata (fileunit, thermometersDailyMean) !check spatial reference system IF ( .NOT. thermometersDailyMean % mapping == mask % grid_mapping) THEN CALL Catch ('info', 'AirTemperatureDailyMean', & 'converting coordinate of stations') !convert stations' coordinate point1 % system = thermometersDailyMean % mapping point2 % system = mask % grid_mapping thermometersDailyMean % mapping = mask % grid_mapping DO i = 1, thermometersDailyMean % countObs point1 % easting = thermometersDailyMean % obs (i) % xyz % easting point1 % northing = thermometersDailyMean % obs (i) % xyz % northing point1 % elevation = thermometersDailyMean % obs (i) % z CALL Convert (point1, point2) thermometersDailyMean % obs (i) % xyz % easting = point2 % easting thermometersDailyMean % obs (i) % xyz % northing = point2 % northing END DO END IF !set supplementary stations network for elevation drift IF (elevationDrift == 1) THEN stationsRefElev = thermometersDailyMean DO i = 1, thermometersDailyMean % countObs stationsRefElev % obs (i) % xyz % elevation = refElevation END DO END IF !allocate vectors if lapse rate has to be computed from data IF ( lapse_rate_computation == 1) THEN ALLOCATE (vectorT (thermometersDailyMean % countObs) ) ALLOCATE (vectorZ (thermometersDailyMean % countObs)) END IF END IF RETURN END SUBROUTINE AirTemperatureDailymeanInit