sfc_climo_gen 1.14.0
Loading...
Searching...
No Matches
driver.F90
Go to the documentation of this file.
1
5
24 program driver
25
26 use model_grid
27 use source_grid
28 use program_setup
29 use utils
30 use esmf
31
32 implicit none
33
34 integer :: rc
35 integer :: localpet
36 integer :: npets
37
38 type(esmf_vm) :: vm
39 type(esmf_regridmethod_flag) :: method
40
41!-------------------------------------------------------------------------
42! Initialize MPI and ESMF environments.
43!-------------------------------------------------------------------------
44
45 call mpi_init(rc)
46
47 print*,"- INITIALIZE ESMF"
48 call esmf_initialize(rc=rc)
49 if(esmf_logfounderror(rctocheck=rc,msg=esmf_logerr_passthru,line=__line__,file=__file__)) &
50 call error_handler("INITIALIZING ESMF", rc)
51
52 print*,"- CALL VMGetGlobal"
53 call esmf_vmgetglobal(vm, rc=rc)
54 if(esmf_logfounderror(rctocheck=rc,msg=esmf_logerr_passthru,line=__line__,file=__file__)) &
55 call error_handler("IN VMGetGlobal.", rc)
56
57 print*,"- CALL VMGet"
58 call esmf_vmget(vm, localpet=localpet, petcount=npets, rc=rc)
59 if(esmf_logfounderror(rctocheck=rc,msg=esmf_logerr_passthru,line=__line__,file=__file__)) &
60 call error_handler("IN VMGet.", rc)
61
62 print*,'- NPETS IS ',npets
63 print*,'- LOCAL PET ',localpet
64
65!-------------------------------------------------------------------------
66! Program set up step.
67!-------------------------------------------------------------------------
68
69 call read_setup_namelist(localpet)
70
71!-------------------------------------------------------------------------
72! Define fv3 model grid.
73!-------------------------------------------------------------------------
74
75 call define_model_grid(localpet, npets)
76
77!-------------------------------------------------------------------------
78! Interpolate fields. Vegetation type must be done first as
79! it defines which points are permanent land ice.
80!-------------------------------------------------------------------------
81
82 call define_source_grid(localpet, npets, input_vegetation_type_file)
83 if (fract_vegsoil_type) then
84 print*,'- WILL OUTPUT VEGETATION TYPE FRACTION AND DOMINANT CATEGORY.'
85 call interp_frac_cats(localpet, input_vegetation_type_file)
86 else
87 print*,'- WILL OUTPUT DOMINANT VEGETATION TYPE.'
88 method=esmf_regridmethod_nearest_stod
89 call interp(localpet, method, input_vegetation_type_file)
90 endif
91 call source_grid_cleanup
92
93! Snow free albedo
94
95 if (trim(input_snowfree_albedo_file) /= "NULL") then
96 call define_source_grid(localpet, npets, input_snowfree_albedo_file)
97 method=esmf_regridmethod_bilinear
98 if (trim(snowfree_albedo_method)=="conserve") method=esmf_regridmethod_conserve
99 call interp(localpet, method, input_snowfree_albedo_file)
100 call source_grid_cleanup
101 endif
102
103! Maximum snow albedo
104
105 if (trim(input_maximum_snow_albedo_file) /= "NULL") then
106 call define_source_grid(localpet, npets, input_maximum_snow_albedo_file)
107 method=esmf_regridmethod_bilinear
108 if (trim(maximum_snow_albedo_method)=="conserve") method=esmf_regridmethod_conserve
109 call interp(localpet, method, input_maximum_snow_albedo_file)
110 call source_grid_cleanup
111 endif
112
113! FACSF - fractional coverage for strong zenith angle
114! dependant albedo.
115
116 if (trim(input_facsf_file) /= "NULL") then
117 call define_source_grid(localpet, npets, input_facsf_file)
118 method=esmf_regridmethod_bilinear
119 call interp(localpet, method, input_facsf_file)
120 call source_grid_cleanup
121 endif
122
123! Soil substrate temperature
124
125 if (trim(input_substrate_temperature_file) /= "NULL") then
126 call define_source_grid(localpet, npets, input_substrate_temperature_file)
127 method=esmf_regridmethod_bilinear
128 call interp(localpet, method, input_substrate_temperature_file)
129 call source_grid_cleanup
130 endif
131
132! Slope type
133
134 if (trim(input_slope_type_file) /= "NULL") then
135 call define_source_grid(localpet, npets, input_slope_type_file)
136 method=esmf_regridmethod_nearest_stod
137 call interp(localpet, method, input_slope_type_file)
138 call source_grid_cleanup
139 endif
140
141! Soil type
142
143 if (trim(input_soil_type_file) /= "NULL") then
144 call define_source_grid(localpet, npets, input_soil_type_file)
145 if (fract_vegsoil_type) then
146 print*,'- WILL OUTPUT SOIL TYPE FRACTION AND DOMINANT CATEGORY.'
147 call interp_frac_cats(localpet, input_soil_type_file)
148 else
149 print*,'- WILL OUTPUT DOMINANT SOIL TYPE.'
150 method=esmf_regridmethod_nearest_stod
151 call interp(localpet, method, input_soil_type_file)
152 endif
153 call source_grid_cleanup
154 endif
155
156! Soil color
157
158 if (trim(input_soil_color_file) /= "NULL") then
159 call define_source_grid(localpet, npets, input_soil_color_file)
160 method=esmf_regridmethod_nearest_stod
161 call interp(localpet, method, input_soil_color_file)
162 call source_grid_cleanup
163 endif
164
165
166! Vegetation greenness
167
168 if (trim(input_vegetation_greenness_file) /= "NULL") then
169 call define_source_grid(localpet, npets, input_vegetation_greenness_file)
170 method=esmf_regridmethod_bilinear
171 if (trim(vegetation_greenness_method)=="conserve") method=esmf_regridmethod_conserve
172 call interp(localpet, method, input_vegetation_greenness_file)
173 call source_grid_cleanup
174 endif
175
176! Leaf Area Index
177
178 if (trim(input_leaf_area_index_file) /= "NULL") then
179 call define_source_grid(localpet, npets, input_leaf_area_index_file)
180 method=esmf_regridmethod_bilinear
181 if (trim(leaf_area_index_method)=="conserve") method=esmf_regridmethod_conserve
182 call interp(localpet, method, input_leaf_area_index_file)
183 call source_grid_cleanup
184 endif
185
186 call model_grid_cleanup
187
188 print*,"- CALL ESMF_finalize"
189 call esmf_finalize(endflag=esmf_end_keepmpi, rc=rc)
190
191 call mpi_finalize(rc)
192
193 print*,'- DONE.'
194 stop
195
196 end program driver
driver
program driver
Reads static surface data on a global lat/lon grid, interpolates the data to the fv3 model grid,...
Definition driver.F90:24
interp
subroutine interp(localpet, method, input_file)
Read the input source data and interpolate it to the model grid.
Definition interp.F90:14
interp_frac_cats
subroutine interp_frac_cats(localpet, input_file)
Read the input source data and interpolate it to the model grid.
Definition interp_frac_cats.F90:15
model_grid
This module defines the model grid.
Definition model_grid.F90:10
model_grid::define_model_grid
subroutine, public define_model_grid(localpet, npets)
Define model grid.
Definition model_grid.F90:62
model_grid::model_grid_cleanup
subroutine, public model_grid_cleanup
Model grid cleanup.
Definition model_grid.F90:424
program_setup
Set up program execution.
Definition program_setup.f90:17
program_setup::input_substrate_temperature_file
character(len=500), public input_substrate_temperature_file
File containing input soil substrate temperature data.
Definition program_setup.f90:27
program_setup::input_leaf_area_index_file
character(len=500), public input_leaf_area_index_file
File containing input leaf area index data.
Definition program_setup.f90:23
program_setup::fract_vegsoil_type
logical, public fract_vegsoil_type
When true, output the percentage of each soil and vegetation type category, and the dominant category...
Definition program_setup.f90:51
program_setup::input_maximum_snow_albedo_file
character(len=500), public input_maximum_snow_albedo_file
File containing input maximum snow albedo data.
Definition program_setup.f90:28
program_setup::maximum_snow_albedo_method
character(len=50), public maximum_snow_albedo_method
Interpolation method for max snow albedo.
Definition program_setup.f90:41
program_setup::input_vegetation_greenness_file
character(len=500), public input_vegetation_greenness_file
File containing input vegetation greenness data.
Definition program_setup.f90:34
program_setup::input_snowfree_albedo_file
character(len=500), public input_snowfree_albedo_file
File containing input snow-free albedo data.
Definition program_setup.f90:29
program_setup::leaf_area_index_method
character(len=50), public leaf_area_index_method
Interpolation method for leaf area index.
Definition program_setup.f90:39
program_setup::input_slope_type_file
character(len=500), public input_slope_type_file
File containing input slope type data.
Definition program_setup.f90:30
program_setup::input_facsf_file
character(len=500), public input_facsf_file
File containing input fractional coverage data for strong zenith angle dependent albedo.
Definition program_setup.f90:24
program_setup::input_soil_type_file
character(len=500), public input_soil_type_file
File containing input soil type data.
Definition program_setup.f90:31
program_setup::read_setup_namelist
subroutine, public read_setup_namelist(localpet)
Read program setup namelist.
Definition program_setup.f90:64
program_setup::input_vegetation_type_file
character(len=500), public input_vegetation_type_file
File containing input vegetation type data.
Definition program_setup.f90:33
program_setup::snowfree_albedo_method
character(len=50), public snowfree_albedo_method
Interpolation method for snowfree albedo.
Definition program_setup.f90:43
program_setup::input_soil_color_file
character(len=500), public input_soil_color_file
File containing input soil color data.
Definition program_setup.f90:32
program_setup::vegetation_greenness_method
character(len=50), public vegetation_greenness_method
Interpolation method for vegetation greenness.
Definition program_setup.f90:45
source_grid
Read grid specs, date information and land/sea mask for the source data that will be interpolated to ...
Definition source_grid.F90:12
source_grid::source_grid_cleanup
subroutine, public source_grid_cleanup
Free up memory associated with this module.
Definition source_grid.F90:404
source_grid::define_source_grid
subroutine, public define_source_grid(localpet, npets, input_file)
Defines esmf grid object for source grid.
Definition source_grid.F90:52
utils
Utilities.
Definition utils.f90:8
utils::error_handler
subroutine, public error_handler(string, rc)
Handle errors.
Definition utils.f90:49