[7541] | 1 | ! ================================================================================================================================= |
---|
| 2 | ! PROGRAM : orchideedriveroasis |
---|
| 3 | ! |
---|
| 4 | ! CONTACT : jan.polcher@lmd.jussieu.fr |
---|
| 5 | ! |
---|
| 6 | ! LICENCE : IPSL (2016) |
---|
| 7 | ! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC |
---|
| 8 | ! |
---|
| 9 | !>\BRIEF This is the main program for the driver which gets the forcing data from OASIS. It is used as an interface |
---|
| 10 | !! for WRF for instance. It allows to have another domain decomposition than the atmosphere. Particularly useful |
---|
| 11 | !! for coupling to atmospheric models not decomposed by the "apple method". |
---|
| 12 | !! |
---|
| 13 | !!\n DESCRIPTION: This program only organises the data and calls sechiba_main after having received the data from OASIS. |
---|
| 14 | !! The main work for the OASIS interface is done in orchoasis_tools.f90 module. |
---|
| 15 | !! Call the various modules to get the forcing data and provide it to SECHIBA. The only complexity |
---|
| 16 | !! is setting-up the domain decomposition and distributing the grid information. |
---|
| 17 | !! The code is parallel from tip to toe using the domain decomposition inherited from LMDZ. |
---|
| 18 | !! |
---|
| 19 | !! RECENT CHANGE(S): None |
---|
| 20 | !! |
---|
| 21 | !! REFERENCE(S) : |
---|
| 22 | !! |
---|
| 23 | !! SVN : |
---|
| 24 | !! $HeadURL: $ |
---|
| 25 | !! $Date: $ |
---|
| 26 | !! $Revision: $ |
---|
| 27 | !! \n |
---|
| 28 | !_ ================================================================================================================================ |
---|
| 29 | PROGRAM orchideeoasis |
---|
| 30 | !--------------------------------------------------------------------- |
---|
| 31 | !- |
---|
| 32 | !- |
---|
| 33 | !--------------------------------------------------------------------- |
---|
| 34 | USE defprec |
---|
| 35 | USE netcdf |
---|
| 36 | ! |
---|
| 37 | ! |
---|
| 38 | USE ioipsl_para |
---|
| 39 | USE mod_orchidee_para |
---|
| 40 | ! |
---|
| 41 | USE grid |
---|
| 42 | USE timer |
---|
| 43 | |
---|
| 44 | USE globgrd |
---|
| 45 | USE orchoasis_tools |
---|
| 46 | ! |
---|
| 47 | USE sechiba |
---|
| 48 | USE control |
---|
| 49 | USE ioipslctrl |
---|
| 50 | ! |
---|
| 51 | USE mod_oasis |
---|
| 52 | !- |
---|
| 53 | IMPLICIT NONE |
---|
| 54 | !- |
---|
| 55 | INCLUDE 'mpif.h' |
---|
| 56 | !- |
---|
| 57 | CHARACTER(LEN=80) :: gridfilename |
---|
| 58 | CHARACTER(LEN=8) :: model_guess |
---|
| 59 | INTEGER(i_std) :: iim_glo, jjm_glo, file_id |
---|
| 60 | !- |
---|
| 61 | INTEGER(i_std) :: nbseg |
---|
| 62 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:) :: lon_glo, lat_glo, area_glo |
---|
| 63 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:) :: mask_glo |
---|
| 64 | INTEGER(i_std), ALLOCATABLE, DIMENSION(:,:) :: maskinv_glo |
---|
| 65 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:,:,:) :: corners_glo |
---|
| 66 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:,:) :: corners_lon, corners_lat |
---|
| 67 | INTEGER(i_std) :: nbindex_g, kjpindex |
---|
| 68 | INTEGER(i_std), ALLOCATABLE, DIMENSION(:) :: kindex, kindex_g |
---|
| 69 | ! |
---|
| 70 | ! Variables for the global grid available on all procs and used |
---|
| 71 | ! to fill the ORCHIDEE variable on the root_proc |
---|
| 72 | ! |
---|
| 73 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:) :: lalo_glo |
---|
| 74 | REAL(r_std), ALLOCATABLE, DIMENSION(:) :: contfrac_glo |
---|
| 75 | CHARACTER(LEN=20) :: calendar |
---|
| 76 | !- |
---|
| 77 | !- Variables local to each processors. |
---|
| 78 | !- |
---|
| 79 | INTEGER(i_std) :: i, j, ik, nbdt, first_point |
---|
| 80 | INTEGER(i_std) :: itau, itau_offset, itau_sechiba |
---|
| 81 | REAL(r_std) :: date0, date0_shifted, dt, julian, julian0 |
---|
| 82 | INTEGER(i_std) :: rest_id, rest_id_stom |
---|
| 83 | INTEGER(i_std) :: hist_id, hist2_id, hist_id_stom, hist_id_stom_IPCC |
---|
| 84 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:) :: lalo_loc |
---|
| 85 | INTEGER(i_std) :: iim, jjm |
---|
| 86 | REAL(r_std), ALLOCATABLE, DIMENSION(:,:) :: lon, lat |
---|
| 87 | !- |
---|
| 88 | !- input fields |
---|
| 89 | !- |
---|
| 90 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: u !! Lowest level wind speed |
---|
| 91 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: v !! Lowest level wind speed |
---|
| 92 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: zlev_uv !! Height of first layer |
---|
| 93 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: zlev_tq !! Height of first layer |
---|
| 94 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: qair !! Lowest level specific humidity |
---|
| 95 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: precip_rain !! Rain precipitation |
---|
| 96 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: precip_snow !! Snow precipitation |
---|
| 97 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: lwdown !! Down-welling long-wave flux |
---|
| 98 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: swdown !! Downwelling surface short-wave flux |
---|
| 99 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: sinang !! cosine of solar zenith angle |
---|
| 100 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: temp_air !! Air temperature in Kelvin |
---|
| 101 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: epot_air !! Air potential energy |
---|
| 102 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: ccanopy !! CO2 concentration in the canopy |
---|
| 103 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: petAcoef !! Coeficients A from the PBL resolution |
---|
| 104 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: peqAcoef !! One for T and another for q |
---|
| 105 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: petBcoef !! Coeficients B from the PBL resolution |
---|
| 106 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: peqBcoef !! One for T and another for q |
---|
| 107 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: cdrag !! Cdrag |
---|
| 108 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: pb !! Lowest level pressure |
---|
| 109 | !- |
---|
| 110 | !- output fields |
---|
| 111 | !- |
---|
| 112 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: z0 !! Surface roughness |
---|
| 113 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: coastalflow !! Diffuse flow of water into the ocean (m^3/dt) |
---|
| 114 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: riverflow !! Largest rivers flowing into the ocean (m^3/dt) |
---|
| 115 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: tsol_rad !! Radiative surface temperature |
---|
| 116 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: vevapp !! Total of evaporation |
---|
| 117 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: temp_sol_new !! New soil temperature |
---|
| 118 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: qsurf !! Surface specific humidity |
---|
| 119 | REAL(r_std), ALLOCATABLE, DIMENSION (:,:) :: albedo !! Albedo |
---|
| 120 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: fluxsens !! Sensible chaleur flux |
---|
| 121 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: fluxlat !! Latent chaleur flux |
---|
| 122 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: emis !! Emissivity |
---|
| 123 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: netco2 !! netco2flux: Sum CO2 flux over PFTs |
---|
| 124 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: carblu !! fco2_lu: Land Cover Change CO2 flux |
---|
| 125 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: carbwh !! fco2_wh: Wood harvest CO2 flux |
---|
| 126 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: carbha !! fco2_ha: Crop harvest CO2 flux |
---|
| 127 | REAL(r_std), ALLOCATABLE, DIMENSION (:,:) :: veget !! Fraction of vegetation type (unitless, 0-1) |
---|
| 128 | REAL(r_std), ALLOCATABLE, DIMENSION (:,:) :: lai !! Leaf area index (m^2 m^{-2} |
---|
| 129 | REAL(r_std), ALLOCATABLE, DIMENSION (:,:) :: height !! Vegetation Height (m) |
---|
| 130 | !- |
---|
| 131 | !- Declarations for OASIS |
---|
| 132 | !- |
---|
| 133 | CHARACTER(LEN=6) :: comp_name = 'orchid' |
---|
| 134 | INTEGER(i_std) :: loc_size, glo_size ! number of pes used |
---|
| 135 | INTEGER(i_std) :: loc_rank, glo_rank ! rank of current pe |
---|
| 136 | INTEGER(i_std) :: localComm, LOCAL_OASIS_COMM ! local MPI communicator and Initialized |
---|
| 137 | INTEGER(i_std) :: comp_id ! component identification |
---|
| 138 | INTEGER(i_std) :: ierror, flag, oasis_info |
---|
| 139 | CHARACTER(LEN=3) :: chout |
---|
| 140 | CHARACTER(LEN=4) :: oasis_gridname |
---|
| 141 | CHARACTER(LEN=30) :: diagfilename |
---|
| 142 | !- |
---|
| 143 | !- |
---|
| 144 | !- |
---|
| 145 | REAL(r_std) :: atmco2 |
---|
| 146 | REAL(r_std), ALLOCATABLE, DIMENSION (:) :: u_tq, v_tq, swnet |
---|
| 147 | LOGICAL :: lrestart_read = .TRUE. !! Logical for _restart_ file to read |
---|
| 148 | LOGICAL :: lrestart_write = .FALSE. !! Logical for _restart_ file to write' |
---|
| 149 | ! |
---|
| 150 | ! Time variables |
---|
| 151 | ! |
---|
| 152 | LOGICAL, PARAMETER :: timemeasure=.TRUE. |
---|
| 153 | REAL(r_std) :: waitput_cputime=0.0, waitget_cputime=0.0, orchidee_cputime=0.0 |
---|
| 154 | REAL(r_std) :: waitput_walltime=0.0, waitget_walltime=0.0, orchidee_walltime=0.0 |
---|
| 155 | ! |
---|
| 156 | ! Print point |
---|
| 157 | ! |
---|
| 158 | !! REAL(r_std), DIMENSION(2) :: testpt=(/44.8,-25.3/) |
---|
| 159 | !! REAL(r_std), DIMENSION(2) :: testpt=(/44.8,-18.3/) |
---|
| 160 | !! REAL(r_std), DIMENSION(2) :: testpt=(/-60.25,-5.25/) |
---|
| 161 | !! REAL(r_std), DIMENSION(2) :: testpt=(/46.7,10.3/) |
---|
| 162 | !! REAL(r_std), DIMENSION(2) :: testpt=(/0.25,49.25/) |
---|
| 163 | ! Case when no ouput is desired. |
---|
| 164 | REAL(r_std), DIMENSION(2) :: testpt=(/9999.99,9999.99/) |
---|
| 165 | INTEGER(i_std) :: ktest |
---|
| 166 | ! |
---|
| 167 | !- |
---|
| 168 | !--------------------------------------------------------------------------------------- |
---|
| 169 | !- |
---|
| 170 | !- Define MPI communicator and set-up OASIS |
---|
| 171 | !- |
---|
| 172 | !--------------------------------------------------------------------------------------- |
---|
| 173 | !- |
---|
| 174 | ! |
---|
| 175 | CALL oasis_init_comp(comp_id, comp_name, ierror) |
---|
| 176 | CALL oasis_get_localcomm (LOCAL_OASIS_COMM, ierror) |
---|
| 177 | ! |
---|
| 178 | ! Set parallel processing in ORCHIDEE |
---|
| 179 | ! |
---|
| 180 | CALL Init_orchidee_para(LOCAL_OASIS_COMM) |
---|
| 181 | !==================================================================================== |
---|
| 182 | ! |
---|
| 183 | ! Start timer now that the paralelisation is started. |
---|
| 184 | ! |
---|
| 185 | IF ( timemeasure ) THEN |
---|
| 186 | CALL init_timer |
---|
| 187 | CALL start_timer(timer_global) |
---|
| 188 | CALL start_timer(timer_mpi) |
---|
| 189 | ENDIF |
---|
| 190 | !- |
---|
| 191 | ! |
---|
| 192 | !--------------------------------------------------------------------------------------- |
---|
| 193 | !- |
---|
| 194 | !- Print some diagnostics for this main of ORCHIDEE |
---|
| 195 | !- |
---|
| 196 | !--------------------------------------------------------------------------------------- |
---|
| 197 | !- |
---|
| 198 | CALL MPI_COMM_RANK(MPI_COMM_WORLD, glo_rank, ierror) |
---|
| 199 | CALL MPI_COMM_RANK(LOCAL_OASIS_COMM, loc_rank, ierror) |
---|
| 200 | CALL MPI_COMM_SIZE(MPI_COMM_WORLD, glo_size, ierror) |
---|
| 201 | CALL MPI_COMM_SIZE(LOCAL_OASIS_COMM, loc_size, ierror) |
---|
| 202 | !- |
---|
| 203 | ! |
---|
| 204 | WRITE (numout,*) '-----------------------------------------------------------' |
---|
| 205 | WRITE (numout,*) '-------------------- comp_name= ',comp_name,' ---------------' |
---|
| 206 | WRITE (numout,*) '-----------------------------------------------------------' |
---|
| 207 | WRITE (numout,*) TRIM(comp_name), ' Running with reals compiled as kind =',r_std |
---|
| 208 | WRITE (numout,*) 'I am component ', TRIM(comp_name), ' local rank :', loc_rank, & |
---|
| 209 | & " Global rank :", glo_rank |
---|
| 210 | WRITE (numout,*) 'CPUs we are using for the main :', loc_size, ' Global number :', glo_size |
---|
| 211 | WRITE (numout,*) 'Local and global MPI communicators :', LOCAL_OASIS_COMM, MPI_COMM_WORLD |
---|
| 212 | WRITE (numout,*) '----------------------------------------------------------' |
---|
| 213 | CALL flush(numout) |
---|
| 214 | ! |
---|
| 215 | !--------------------------------------------------------------------------------------- |
---|
| 216 | !- |
---|
| 217 | !- Start the getconf processes |
---|
| 218 | !- |
---|
| 219 | !--------------------------------------------------------------------------------------- |
---|
| 220 | !- |
---|
| 221 | !! CALL getin_name("run.def") |
---|
| 222 | !- |
---|
| 223 | !Config Key = GRID_FILE |
---|
| 224 | !Config Desc = Name of file containing the forcing data |
---|
| 225 | !Config If = [-] |
---|
| 226 | !Config Def = grid_file.nc |
---|
| 227 | !Config Help = This is the name of the file from which we will read |
---|
| 228 | !Config or write into it the description of the grid from |
---|
| 229 | !Config the forcing file. |
---|
| 230 | !Config compliant. |
---|
| 231 | !Config Units = [FILE] |
---|
| 232 | !- |
---|
| 233 | gridfilename='grid_file.nc' |
---|
| 234 | CALL getin_p('GRID_FILE', gridfilename) |
---|
| 235 | !- |
---|
| 236 | !- Get some basic variables from the run.def |
---|
| 237 | !- |
---|
| 238 | atmco2=350. |
---|
| 239 | CALL getin_p('ATM_CO2',atmco2) |
---|
| 240 | !--------------------------------------------------------------------------------------- |
---|
| 241 | !- |
---|
| 242 | !- Get the grid on all processors. |
---|
| 243 | !- |
---|
| 244 | !--------------------------------------------------------------------------------------- |
---|
| 245 | !- |
---|
| 246 | ! |
---|
| 247 | CALL globgrd_getdomsz(gridfilename, iim_glo, jjm_glo, nbindex_g, model_guess, file_id) |
---|
| 248 | nbseg = 4 |
---|
| 249 | ! |
---|
| 250 | !- |
---|
| 251 | !- Allocation of memory |
---|
| 252 | !- variables over the entire grid (thus in x,y) |
---|
| 253 | ALLOCATE(lon_glo(iim_glo, jjm_glo)) |
---|
| 254 | ALLOCATE(lat_glo(iim_glo, jjm_glo)) |
---|
| 255 | ALLOCATE(mask_glo(iim_glo, jjm_glo)) |
---|
| 256 | ALLOCATE(area_glo(iim_glo, jjm_glo)) |
---|
| 257 | ALLOCATE(corners_glo(iim_glo, jjm_glo, nbseg, 2)) |
---|
| 258 | ! |
---|
| 259 | ! Gathered variables |
---|
| 260 | ALLOCATE(kindex_g(nbindex_g)) |
---|
| 261 | ALLOCATE(contfrac_glo(nbindex_g)) |
---|
| 262 | !- |
---|
| 263 | !- |
---|
| 264 | !- |
---|
| 265 | CALL globgrd_getgrid(file_id, iim_glo, jjm_glo, nbindex_g, model_guess, & |
---|
| 266 | & lon_glo, lat_glo, mask_glo, area_glo, corners_glo,& |
---|
| 267 | & kindex_g, contfrac_glo, calendar) |
---|
| 268 | !- |
---|
| 269 | !- Set the calendar and get some information |
---|
| 270 | !- |
---|
| 271 | CALL ioconf_calendar(calendar) |
---|
| 272 | CALL ioget_calendar(one_year, one_day) |
---|
| 273 | !- |
---|
| 274 | !- get the time period for the run |
---|
| 275 | !- |
---|
| 276 | CALL orchoasis_time(date0, dt, nbdt) |
---|
| 277 | !- |
---|
| 278 | !- lalo needs to be created before going into the parallel region |
---|
| 279 | !- |
---|
| 280 | ALLOCATE(lalo_glo(nbindex_g,2)) |
---|
| 281 | DO ik=1,nbindex_g |
---|
| 282 | ! |
---|
| 283 | j = ((kindex_g(ik)-1)/iim_glo)+1 |
---|
| 284 | i = (kindex_g(ik)-(j-1)*iim_glo) |
---|
| 285 | ! |
---|
| 286 | IF ( i > iim_glo .OR. j > jjm_glo ) THEN |
---|
| 287 | WRITE(100+mpi_rank,*) "Error in the indexing (ik, kindex, i, j) : ", ik, kindex(ik), i, j |
---|
| 288 | STOP "ERROR in orchideeoasis" |
---|
| 289 | ENDIF |
---|
| 290 | ! |
---|
| 291 | lalo_glo(ik,1) = lat_glo(i,j) |
---|
| 292 | lalo_glo(ik,2) = lon_glo(i,j) |
---|
| 293 | ! |
---|
| 294 | ENDDO |
---|
| 295 | ! |
---|
| 296 | WRITE(100+mpi_rank,*) "Rank", mpi_rank, " Before parallel region All land points : ", nbindex_g |
---|
| 297 | WRITE(100+mpi_rank,*) "Rank", mpi_rank, " from ", iim_glo, " point in Lon. and ", jjm_glo, "in Lat." |
---|
| 298 | !- |
---|
| 299 | !--------------------------------------------------------------------------------------- |
---|
| 300 | !- |
---|
| 301 | !- Initialise the ORCHIDEE domain on the procs (longitude, latitude, indices) |
---|
| 302 | !- |
---|
| 303 | !--------------------------------------------------------------------------------------- |
---|
| 304 | !- |
---|
| 305 | !- init_data_para also transfers kindex_g to index_g (the variable used in ORCHIDEE) |
---|
| 306 | !- |
---|
| 307 | CALL grid_set_glo(iim_glo, jjm_glo, nbindex_g) |
---|
| 308 | CALL grid_allocate_glo(nbseg) |
---|
| 309 | ! Copy the list of indexes of land points into index_g used by ORCHIDEE and then broacast to all |
---|
| 310 | ! processors |
---|
| 311 | CALL bcast(nbindex_g) |
---|
| 312 | IF ( is_root_prc) index_g = kindex_g |
---|
| 313 | CALL bcast(index_g) |
---|
| 314 | ! |
---|
| 315 | WRITE(numout,*) "Rank", mpi_rank, "Into Init_orchidee_data_para_driver with ", nbindex_g,index_g(1) |
---|
| 316 | ! |
---|
| 317 | CALL Init_orchidee_data_para_driver(nbindex_g,index_g) |
---|
| 318 | CALL init_ioipsl_para |
---|
| 319 | ! |
---|
| 320 | WRITE(numout,*) "Rank", mpi_rank, "After init_data_para global size : ", nbp_glo, SIZE(index_g), iim_g, iim_glo, jjm_g, jjm_glo |
---|
| 321 | WRITE(numout,'("After init_data_para local : ij_nb, jj_nb",2I4)') iim_glo, jj_nb |
---|
| 322 | ! |
---|
| 323 | ! Allocate grid on the local processor |
---|
| 324 | ! |
---|
| 325 | IF ( model_guess == "regular") THEN |
---|
| 326 | CALL grid_init (nbp_loc, nbseg, regular_lonlat, "ForcingGrid") |
---|
| 327 | ELSE IF ( model_guess == "WRF") THEN |
---|
| 328 | CALL grid_init (nbp_loc, nbseg, regular_xy, "WRFGrid") |
---|
| 329 | ELSE |
---|
| 330 | CALL ipslerr(3, "orchideeoasis", "The grid found in the GRID_FILE is not supported by ORCHIDEE", "", "") |
---|
| 331 | ENDIF |
---|
| 332 | ! |
---|
| 333 | ! Transfer the global grid variables to the ORCHIDEE version on the root proc |
---|
| 334 | ! *_glo -> *_g |
---|
| 335 | ! Variables *_g were allocated with the CALL init_grid |
---|
| 336 | ! |
---|
| 337 | IF ( is_root_prc) THEN |
---|
| 338 | ! |
---|
| 339 | lalo_g(:,:) = lalo_glo(:,:) |
---|
| 340 | lon_g(:,:) = lon_glo(:,:) |
---|
| 341 | lat_g(:,:) = lat_glo(:,:) |
---|
| 342 | ! |
---|
| 343 | ENDIF |
---|
| 344 | !- |
---|
| 345 | ! |
---|
| 346 | ! Set the local dimensions of the fields |
---|
| 347 | ! |
---|
| 348 | iim = iim_glo |
---|
| 349 | jjm = jj_nb |
---|
| 350 | kjpindex = nbp_loc |
---|
| 351 | ! |
---|
| 352 | WRITE(numout,*) mpi_rank, "DIMENSIONS of grid on processor : iim, jjm, kjpindex = ", iim, jjm, kjpindex |
---|
| 353 | ! |
---|
| 354 | ! Allocate the local arrays we need : |
---|
| 355 | ! |
---|
| 356 | ALLOCATE(lon(iim,jjm), lat(iim,jjm)) |
---|
| 357 | ALLOCATE(kindex(kjpindex)) |
---|
| 358 | ! |
---|
| 359 | lon=lon_glo(:,jj_para_begin(mpi_rank):jj_para_end(mpi_rank)) |
---|
| 360 | lat=lat_glo(:,jj_para_begin(mpi_rank):jj_para_end(mpi_rank)) |
---|
| 361 | ! |
---|
| 362 | ! |
---|
| 363 | ! Redistribute the indeces on all procs (apple distribution of land points) |
---|
| 364 | ! |
---|
| 365 | CALL bcast(lon_g) |
---|
| 366 | CALL bcast(lat_g) |
---|
| 367 | CALL scatter(index_g, kindex) |
---|
| 368 | ! |
---|
| 369 | ! |
---|
| 370 | ! Apply the offset needed so that kindex refers to the index of the land point |
---|
| 371 | ! on the current region, i.e. the local lon lat domain. |
---|
| 372 | ! |
---|
| 373 | kindex(1:kjpindex)=kindex(1:kjpindex)-(jj_begin-1)*iim_glo |
---|
| 374 | ! |
---|
| 375 | ! This routine transforms the global grid into a series of polygons for all land |
---|
| 376 | ! points identified by index_g. |
---|
| 377 | ! |
---|
| 378 | CALL grid_stuff(nbindex_g, iim_g, jjm_g, lon_g, lat_g, index_g, contfrac_glo) |
---|
| 379 | ! |
---|
| 380 | ! Distribute the global lalo to the local processor level lalo |
---|
| 381 | ! |
---|
| 382 | ALLOCATE(lalo_loc(kjpindex,2)) |
---|
| 383 | CALL scatter(lalo_glo, lalo_loc) |
---|
| 384 | lalo(:,:) = lalo_loc(:,:) |
---|
| 385 | ! |
---|
| 386 | !- |
---|
| 387 | !--------------------------------------------------------------------------------------- |
---|
| 388 | !- |
---|
| 389 | !- Allocate the space for the per processor coupling variables |
---|
| 390 | !- |
---|
| 391 | !--------------------------------------------------------------------------------------- |
---|
| 392 | !- |
---|
| 393 | ALLOCATE(zlev_tq(kjpindex), zlev_uv(kjpindex)) |
---|
| 394 | ALLOCATE(u(kjpindex), v(kjpindex), pb(kjpindex)) |
---|
| 395 | ALLOCATE(temp_air(kjpindex)) |
---|
| 396 | ALLOCATE(qair(kjpindex)) |
---|
| 397 | ALLOCATE(precip_rain(kjpindex), precip_snow(kjpindex)) |
---|
| 398 | ALLOCATE(swdown(kjpindex), lwdown(kjpindex), sinang(kjpindex)) |
---|
| 399 | ! |
---|
| 400 | ALLOCATE(epot_air(kjpindex), ccanopy(kjpindex), cdrag(kjpindex), swnet(kjpindex)) |
---|
| 401 | ! |
---|
| 402 | ALLOCATE(petAcoef(kjpindex), peqAcoef(kjpindex), petBcoef(kjpindex), peqBcoef(kjpindex)) |
---|
| 403 | ALLOCATE(u_tq(kjpindex), v_tq(kjpindex)) |
---|
| 404 | ! |
---|
| 405 | ALLOCATE(vevapp(kjpindex), fluxsens(kjpindex), fluxlat(kjpindex), coastalflow(kjpindex), riverflow(kjpindex)) |
---|
| 406 | ALLOCATE(netco2(kjpindex), carblu(kjpindex)) |
---|
| 407 | ALLOCATE(carbwh(kjpindex), carbha(kjpindex)) |
---|
| 408 | ALLOCATE(tsol_rad(kjpindex), temp_sol_new(kjpindex), qsurf(kjpindex)) |
---|
| 409 | ALLOCATE(albedo(kjpindex,2), emis(kjpindex), z0(kjpindex)) |
---|
| 410 | ALLOCATE(veget(kjpindex,nvm)) |
---|
| 411 | ALLOCATE(lai(kjpindex,nvm)) |
---|
| 412 | ALLOCATE(height(kjpindex,nvm)) |
---|
| 413 | ! |
---|
| 414 | WRITE(numout,*) "Rank", mpi_rank, "Domain size per proc !:", iim_glo, jjm_glo, kjpindex, SIZE(kindex) |
---|
| 415 | WRITE(numout,*) "Rank", mpi_rank, "In parallel region land index starts at : ", kindex(1) |
---|
| 416 | ! |
---|
| 417 | CALL orchoasis_time(date0, dt, nbdt) |
---|
| 418 | ! |
---|
| 419 | CALL control_initialize |
---|
| 420 | dt_sechiba = dt |
---|
| 421 | ! |
---|
| 422 | itau = 0 |
---|
| 423 | ! |
---|
| 424 | CALL ioipslctrl_restini(itau, date0, dt, rest_id, rest_id_stom, itau_offset, date0_shifted) |
---|
| 425 | ! |
---|
| 426 | ! Get the date of the first time step |
---|
| 427 | ! |
---|
| 428 | julian = date0 + 0.5*(dt/one_day) |
---|
| 429 | CALL ju2ymds (julian, year, month, day, sec) |
---|
| 430 | ! |
---|
| 431 | in_julian = itau2date(itau, date0, dt) |
---|
| 432 | CALL ymds2ju (year,1,1,zero, julian0) |
---|
| 433 | julian_diff = in_julian-julian0 |
---|
| 434 | |
---|
| 435 | CALL xios_orchidee_init( MPI_COMM_ORCH, & |
---|
| 436 | date0, year, month, day, julian_diff, & |
---|
| 437 | lon, lat ) |
---|
| 438 | ! |
---|
| 439 | itau_sechiba = itau + itau_offset |
---|
| 440 | ! |
---|
| 441 | WRITE(numout,*) "orchoasis_history :", iim, jjm, SHAPE(lon) |
---|
| 442 | !- Initialize IOIPSL sechiba output files |
---|
| 443 | CALL ioipslctrl_history(iim, jjm, lon, lat, kindex, kjpindex, itau_sechiba, & |
---|
| 444 | date0, dt, hist_id, hist2_id, hist_id_stom, hist_id_stom_IPCC) |
---|
| 445 | WRITE(numout,*) "HISTORY : Define for ", itau, date0, dt |
---|
| 446 | ! |
---|
| 447 | !- |
---|
| 448 | !- |
---|
| 449 | !--------------------------------------------------------------------------------------- |
---|
| 450 | !- |
---|
| 451 | !- Send the grid to OASIS ... only on the root processor |
---|
| 452 | !- |
---|
| 453 | !--------------------------------------------------------------------------------------- |
---|
| 454 | !- |
---|
| 455 | IF ( is_root_prc ) THEN |
---|
| 456 | ! |
---|
| 457 | ! TOCOMPLETE - Put here OASIS grid, corner, areas and mask writing calls ! |
---|
| 458 | ! |
---|
| 459 | oasis_gridname = model_guess(1:4) |
---|
| 460 | CALL oasis_start_grids_writing(flag) |
---|
| 461 | ! |
---|
| 462 | CALL oasis_write_grid(oasis_gridname, iim_glo, jjm_glo, lon_glo, lat_glo) |
---|
| 463 | ! |
---|
| 464 | ALLOCATE(corners_lon(iim_glo, jjm_glo, nbseg), corners_lat(iim_glo, jjm_glo, nbseg)) |
---|
| 465 | corners_lon(:,:,:) = corners_glo(:,:,:,1) |
---|
| 466 | corners_lat(:,:,:) = corners_glo(:,:,:,2) |
---|
| 467 | CALL oasis_write_corner(oasis_gridname, iim_glo, jjm_glo, nbseg, corners_lon, corners_lat) |
---|
| 468 | ! |
---|
| 469 | ALLOCATE(maskinv_glo(iim_glo, jjm_glo)) |
---|
| 470 | DO i=1,iim_glo |
---|
| 471 | DO j=1,jjm_glo |
---|
| 472 | IF (mask_glo(i,j) == 0) THEN |
---|
| 473 | maskinv_glo(i,j) = 1 |
---|
| 474 | ELSE |
---|
| 475 | maskinv_glo(i,j) = 0 |
---|
| 476 | ENDIF |
---|
| 477 | ENDDO |
---|
| 478 | ENDDO |
---|
| 479 | CALL oasis_write_mask(oasis_gridname, iim_glo, jjm_glo, maskinv_glo) |
---|
| 480 | ! |
---|
| 481 | CALL oasis_write_area(oasis_gridname, iim_glo, jjm_glo, area_glo) |
---|
| 482 | ! |
---|
| 483 | CALL oasis_terminate_grids_writing() |
---|
| 484 | ! |
---|
| 485 | IF (model_guess == "WRF") THEN |
---|
| 486 | oasis_gridname = "twrf" |
---|
| 487 | ! |
---|
| 488 | CALL oasis_start_grids_writing(flag) |
---|
| 489 | ! |
---|
| 490 | CALL oasis_write_grid(oasis_gridname, iim_glo, jjm_glo, lon_glo, lat_glo) |
---|
| 491 | ! |
---|
| 492 | CALL oasis_write_corner(oasis_gridname, iim_glo, jjm_glo, nbseg, corners_lon, corners_lat) |
---|
| 493 | ! |
---|
| 494 | CALL oasis_write_mask(oasis_gridname, iim_glo, jjm_glo, maskinv_glo) |
---|
| 495 | ! |
---|
| 496 | CALL oasis_write_area(oasis_gridname, iim_glo, jjm_glo, area_glo) |
---|
| 497 | ! |
---|
| 498 | CALL oasis_terminate_grids_writing() |
---|
| 499 | ENDIF |
---|
| 500 | DEALLOCATE(corners_lon, corners_lat) |
---|
| 501 | ! |
---|
| 502 | ENDIF |
---|
| 503 | ! |
---|
| 504 | call flush(numout) |
---|
| 505 | !- |
---|
| 506 | !--------------------------------------------------------------------------------------- |
---|
| 507 | !- |
---|
| 508 | !- Declare the variables to be exchanged through OASIS |
---|
| 509 | !- |
---|
| 510 | !--------------------------------------------------------------------------------------- |
---|
| 511 | !- |
---|
| 512 | CALL orchoasis_defvar(mpi_rank, kjpindex) |
---|
| 513 | ! |
---|
| 514 | !- |
---|
| 515 | !--------------------------------------------------------------------------------------- |
---|
| 516 | !- |
---|
| 517 | !- Get a first set of forcing data |
---|
| 518 | !- |
---|
| 519 | !--------------------------------------------------------------------------------------- |
---|
| 520 | !- |
---|
| 521 | itau = 0 |
---|
| 522 | julian = date0 + (itau+0.5)*(dt/one_day) |
---|
| 523 | CALL ju2ymds (julian, year, month, day, sec) |
---|
| 524 | !- |
---|
| 525 | !--------------------------------------------------------------------------------------- |
---|
| 526 | !- |
---|
| 527 | !- Going into the time loop |
---|
| 528 | !- |
---|
| 529 | !--------------------------------------------------------------------------------------- |
---|
| 530 | !- |
---|
| 531 | ! |
---|
| 532 | DO itau=1,nbdt |
---|
| 533 | ! |
---|
| 534 | julian = date0 + (itau-0.5)*(dt/one_day) |
---|
| 535 | ! Needed for STOMATE but should be taken from the arguments of SECHIBA |
---|
| 536 | in_julian = itau2date(itau, date0, dt) |
---|
| 537 | ! |
---|
| 538 | CALL ju2ymds (julian, year, month, day, sec) |
---|
| 539 | CALL ymds2ju (year,1,1,zero, julian0) |
---|
| 540 | julian_diff = in_julian-julian0 |
---|
| 541 | ! |
---|
| 542 | ! |
---|
| 543 | IF ( itau == nbdt ) lrestart_write = .TRUE. |
---|
| 544 | ! |
---|
| 545 | !--------------------------------------------------------------------------------------- |
---|
| 546 | !- |
---|
| 547 | !- Get the variables from OASIS |
---|
| 548 | !- |
---|
| 549 | !--------------------------------------------------------------------------------------- |
---|
| 550 | ! |
---|
| 551 | ! OASIS get call are blocking so they need to be done once all is finsihed |
---|
| 552 | ! |
---|
| 553 | CALL orchoasis_getvar(itau-1, dt, kjpindex, kindex - (ii_begin - 1), & |
---|
| 554 | & zlev_tq, zlev_uv, temp_air, qair, & |
---|
| 555 | & precip_rain, precip_snow, swnet, lwdown, sinang, u, v, pb, cdrag) |
---|
| 556 | ! |
---|
| 557 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, temp_air, "RECEIVED Air temperature") |
---|
| 558 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, qair, "RECEIVED Air humidity") |
---|
| 559 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, precip_rain*one_day, "RECEIVED Rainfall") |
---|
| 560 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, precip_snow*one_day, "RECEIVED Snowfall") |
---|
| 561 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, swnet, "RECEIVED net solar") |
---|
| 562 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, lwdown, "RECEIVED lwdown") |
---|
| 563 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, u, "RECEIVED East-ward wind") |
---|
| 564 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, v, "RECEIVED North-ward wind") |
---|
| 565 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, pb, "RECEIVED surface pressure") |
---|
| 566 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, zlev_uv, "RECEIVED UV height") |
---|
| 567 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, zlev_tq, "RECEIVED TQ height") |
---|
| 568 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, sinang, "RECEIVED sinang") |
---|
| 569 | ! |
---|
| 570 | ! |
---|
| 571 | ! Adaptation of the forcing data to SECHIBA's needs |
---|
| 572 | ! |
---|
| 573 | ! Contrary to what the documentation says, ORCHIDEE expects surface pressure in hPa. |
---|
| 574 | pb(:) = pb(:)/100. |
---|
| 575 | epot_air(:) = cp_air*temp_air(:)+cte_grav*zlev_tq(:) |
---|
| 576 | ccanopy(:) = atmco2 |
---|
| 577 | ! |
---|
| 578 | petBcoef(:) = epot_air(:) |
---|
| 579 | peqBcoef(:) = qair(:) |
---|
| 580 | petAcoef(:) = zero |
---|
| 581 | peqAcoef(:) = zero |
---|
| 582 | ! |
---|
| 583 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, z0, & |
---|
| 584 | & "Z0 before compute", ktest) |
---|
| 585 | ! |
---|
| 586 | ! Interpolate the wind (which is at hight zlev_uv) to the same height |
---|
| 587 | ! as the temperature and humidity (at zlev_tq). |
---|
| 588 | ! |
---|
| 589 | u_tq(:) = u(:)*LOG(zlev_tq(:)/z0(:))/LOG(zlev_uv(:)/z0(:)) |
---|
| 590 | v_tq(:) = v(:)*LOG(zlev_tq(:)/z0(:))/LOG(zlev_uv(:)/z0(:)) |
---|
| 591 | ! |
---|
| 592 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, u_tq, "USED East-ward wind") |
---|
| 593 | ! |
---|
| 594 | IF ( itau .NE. 1 ) THEN |
---|
| 595 | IF ( timemeasure ) THEN |
---|
| 596 | waitget_cputime = waitget_cputime + Get_cpu_Time(timer_global) |
---|
| 597 | waitget_walltime = waitget_walltime + Get_real_Time(timer_global) |
---|
| 598 | CALL stop_timer(timer_global) |
---|
| 599 | CALL start_timer(timer_global) |
---|
| 600 | ENDIF |
---|
| 601 | ENDIF |
---|
| 602 | ! |
---|
| 603 | !--------------------------------------------------------------------------------------- |
---|
| 604 | !- |
---|
| 605 | !- IF first time step : Call to SECHIBA_initialize to set-up ORCHIDEE before doing an actual call |
---|
| 606 | !- which will provide the first fluxes. |
---|
| 607 | !- |
---|
| 608 | !--------------------------------------------------------------------------------------- |
---|
| 609 | ! |
---|
| 610 | itau_sechiba = itau+itau_offset |
---|
| 611 | ! |
---|
| 612 | ! Update the calendar in xios by sending the new time step |
---|
| 613 | CALL xios_orchidee_update_calendar(itau_sechiba) |
---|
| 614 | ! |
---|
| 615 | IF ( itau == 1 ) THEN |
---|
| 616 | ! |
---|
| 617 | IF ( timemeasure ) THEN |
---|
| 618 | WRITE(numout,*) '------> CPU Time for start-up of main : ',Get_cpu_Time(timer_global) |
---|
| 619 | WRITE(numout,*) '------> Real Time for start-up of main : ',Get_real_Time(timer_global) |
---|
| 620 | CALL stop_timer(timer_global) |
---|
| 621 | CALL start_timer(timer_global) |
---|
| 622 | ENDIF |
---|
| 623 | ! |
---|
| 624 | CALL sechiba_initialize( & |
---|
| 625 | itau_sechiba, iim*jjm, kjpindex, kindex, & |
---|
| 626 | lalo_loc, contfrac, neighbours, resolution, zlev_tq, & |
---|
| 627 | u_tq, v_tq, qair, temp_air, & |
---|
| 628 | petAcoef, peqAcoef, petBcoef, peqBcoef, & |
---|
| 629 | precip_rain, precip_snow, lwdown, swnet, swdown, & |
---|
| 630 | pb, rest_id, hist_id, hist2_id, & |
---|
| 631 | rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & |
---|
| 632 | coastalflow, riverflow, tsol_rad, vevapp, qsurf, & |
---|
| 633 | z0, albedo, fluxsens, fluxlat, emis, & |
---|
| 634 | temp_sol_new, cdrag) |
---|
| 635 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, temp_sol_new, "Init temp_sol_new") |
---|
| 636 | ! |
---|
| 637 | ! Use the obtained albedo to diagnose the Downward Solar |
---|
| 638 | ! |
---|
| 639 | swdown(:) = swnet(:)/(1.-(albedo(:,1)+albedo(:,2))/2.) |
---|
| 640 | ! |
---|
| 641 | lrestart_read = .FALSE. |
---|
| 642 | ! |
---|
| 643 | CALL histwrite_p(hist_id, 'LandPoints', itau+1, (/ REAL(kindex) /), kjpindex, kindex) |
---|
| 644 | CALL histwrite_p(hist_id, 'Areas', itau+1, area, kjpindex, kindex) |
---|
| 645 | CALL histwrite_p(hist_id, 'Contfrac', itau+1, contfrac, kjpindex, kindex) |
---|
| 646 | ! |
---|
| 647 | IF ( timemeasure ) THEN |
---|
| 648 | WRITE(numout,*) '------> CPU Time for set-up of ORCHIDEE : ',Get_cpu_Time(timer_global) |
---|
| 649 | WRITE(numout,*) '------> Real Time for set-up of ORCHIDEE : ',Get_real_Time(timer_global) |
---|
| 650 | CALL stop_timer(timer_global) |
---|
| 651 | CALL start_timer(timer_global) |
---|
| 652 | ENDIF |
---|
| 653 | ! |
---|
| 654 | ENDIF |
---|
| 655 | ! |
---|
| 656 | !--------------------------------------------------------------------------------------- |
---|
| 657 | !- |
---|
| 658 | !- Call to SECHIBA |
---|
| 659 | !- |
---|
| 660 | !--------------------------------------------------------------------------------------- |
---|
| 661 | ! |
---|
| 662 | CALL sechiba_main (itau_sechiba, iim*jjm, kjpindex, kindex, & |
---|
| 663 | & lrestart_read, lrestart_write, & |
---|
| 664 | & lalo_loc, contfrac, neighbours, resolution, & |
---|
| 665 | ! First level conditions |
---|
| 666 | & zlev_tq, u_tq, v_tq, qair, temp_air, epot_air, ccanopy, & |
---|
| 667 | ! Variables for the implicit coupling |
---|
| 668 | & cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, & |
---|
| 669 | ! Rain, snow, radiation and surface pressure |
---|
| 670 | & precip_rain ,precip_snow, lwdown, swnet, swdown, sinang, pb, & |
---|
| 671 | ! Output : Fluxes |
---|
| 672 | & vevapp, fluxsens, fluxlat, coastalflow, riverflow, & |
---|
| 673 | ! CO2 fluxes |
---|
| 674 | & netco2, carblu, carbwh, carbha, & |
---|
| 675 | ! Surface temperatures and surface properties |
---|
| 676 | & tsol_rad, temp_sol_new, qsurf, albedo, emis, z0, & |
---|
| 677 | ! Vegeation, lai and vegetation height |
---|
| 678 | & veget, lai, height, & |
---|
| 679 | ! File ids |
---|
| 680 | & rest_id, hist_id, hist2_id, rest_id_stom, hist_id_stom, hist_id_stom_IPCC) |
---|
| 681 | ! |
---|
| 682 | ! Use the obtained albedo to diagnose the Downward Solar |
---|
| 683 | ! |
---|
| 684 | swdown(:) = swnet(:)/(1.-(albedo(:,1)+albedo(:,2))/2.) |
---|
| 685 | ! |
---|
| 686 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, temp_sol_new, "Produced temp_sol_new") |
---|
| 687 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, fluxsens, "Produced fluxsens") |
---|
| 688 | CALL orchoasis_printpoint(julian, testpt(1), testpt(2), kjpindex, lalo_loc, fluxlat, "Produced fluxlat") |
---|
| 689 | ! |
---|
| 690 | IF ( timemeasure ) THEN |
---|
| 691 | orchidee_cputime = orchidee_cputime + Get_cpu_Time(timer_global) |
---|
| 692 | orchidee_walltime = orchidee_walltime + Get_real_Time(timer_global) |
---|
| 693 | CALL stop_timer(timer_global) |
---|
| 694 | CALL start_timer(timer_global) |
---|
| 695 | ENDIF |
---|
| 696 | ! |
---|
| 697 | !--------------------------------------------------------------------------------------- |
---|
| 698 | ! Send the ORCHIDEE output back to the driver |
---|
| 699 | !--------------------------------------------------------------------------------------- |
---|
| 700 | ! |
---|
| 701 | CALL orchoasis_putvar(itau, dt, kjpindex, kindex - (ii_begin - 1),& |
---|
| 702 | & vevapp, fluxsens, fluxlat, coastalflow, riverflow, & |
---|
| 703 | & netco2, carblu, tsol_rad, temp_sol_new, qsurf, albedo, emis, z0) |
---|
| 704 | ! |
---|
| 705 | IF ( timemeasure ) THEN |
---|
| 706 | waitput_cputime = waitput_cputime + Get_cpu_Time(timer_global) |
---|
| 707 | waitput_walltime = waitput_walltime + Get_real_Time(timer_global) |
---|
| 708 | CALL stop_timer(timer_global) |
---|
| 709 | CALL start_timer(timer_global) |
---|
| 710 | ENDIF |
---|
| 711 | ! |
---|
| 712 | !--------------------------------------------------------------------------------------- |
---|
| 713 | !- |
---|
| 714 | !- Write diagnostics |
---|
| 715 | !- |
---|
| 716 | !--------------------------------------------------------------------------------------- |
---|
| 717 | ! |
---|
| 718 | !--------------------------------------------------------------------------------------- |
---|
| 719 | !- |
---|
| 720 | !- Write diagnostics |
---|
| 721 | !- |
---|
| 722 | !--------------------------------------------------------------------------------------- |
---|
| 723 | ! |
---|
| 724 | CALL xios_orchidee_send_field("LandPoints" ,(/ ( REAL(ik), ik=1,kjpindex ) /)) |
---|
| 725 | CALL xios_orchidee_send_field("Areas", area) |
---|
| 726 | CALL xios_orchidee_send_field("Contfrac",contfrac) |
---|
| 727 | CALL xios_orchidee_send_field("evap",vevapp*one_day/dt_sechiba) |
---|
| 728 | CALL xios_orchidee_send_field("evap_alma",vevapp/dt_sechiba) |
---|
| 729 | CALL xios_orchidee_send_field("coastalflow",coastalflow/dt_sechiba) |
---|
| 730 | CALL xios_orchidee_send_field("riverflow",riverflow/dt_sechiba) |
---|
| 731 | CALL xios_orchidee_send_field("temp_sol_C",temp_sol_new-ZeroCelsius) |
---|
| 732 | CALL xios_orchidee_send_field("temp_sol_K",temp_sol_new) |
---|
| 733 | CALL xios_orchidee_send_field("fluxsens",fluxsens) |
---|
| 734 | CALL xios_orchidee_send_field("fluxlat",fluxlat) |
---|
| 735 | CALL xios_orchidee_send_field("tair",temp_air) |
---|
| 736 | CALL xios_orchidee_send_field("qair",qair) |
---|
| 737 | CALL xios_orchidee_send_field("swnet",swnet) |
---|
| 738 | CALL xios_orchidee_send_field("swdown",swdown) |
---|
| 739 | ! zpb in hPa, output in Pa |
---|
| 740 | CALL xios_orchidee_send_field("Psurf",pb*100.) |
---|
| 741 | ! |
---|
| 742 | IF ( .NOT. almaoutput ) THEN |
---|
| 743 | ! |
---|
| 744 | ! ORCHIDEE INPUT variables |
---|
| 745 | ! |
---|
| 746 | CALL histwrite_p (hist_id, 'swdown', itau_sechiba, swdown, kjpindex, kindex) |
---|
| 747 | CALL histwrite_p (hist_id, 'tair', itau_sechiba, temp_air, kjpindex, kindex) |
---|
| 748 | CALL histwrite_p (hist_id, 'qair', itau_sechiba, qair, kjpindex, kindex) |
---|
| 749 | CALL histwrite_p (hist_id, 'evap', itau_sechiba, vevapp, kjpindex, kindex) |
---|
| 750 | CALL histwrite_p (hist_id, 'coastalflow',itau_sechiba, coastalflow, kjpindex, kindex) |
---|
| 751 | CALL histwrite_p (hist_id, 'riverflow',itau_sechiba, riverflow, kjpindex, kindex) |
---|
| 752 | ! |
---|
| 753 | CALL histwrite_p (hist_id, 'temp_sol', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 754 | CALL histwrite_p (hist_id, 'tsol_max', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 755 | CALL histwrite_p (hist_id, 'tsol_min', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 756 | CALL histwrite_p (hist_id, 'fluxsens', itau_sechiba, fluxsens, kjpindex, kindex) |
---|
| 757 | CALL histwrite_p (hist_id, 'fluxlat', itau_sechiba, fluxlat, kjpindex, kindex) |
---|
| 758 | CALL histwrite_p (hist_id, 'swnet', itau_sechiba, swnet, kjpindex, kindex) |
---|
| 759 | CALL histwrite_p (hist_id, 'alb_vis', itau_sechiba, albedo(:,1), kjpindex, kindex) |
---|
| 760 | CALL histwrite_p (hist_id, 'alb_nir', itau_sechiba, albedo(:,2), kjpindex, kindex) |
---|
| 761 | ! |
---|
| 762 | IF ( hist2_id > 0 ) THEN |
---|
| 763 | CALL histwrite_p (hist2_id, 'swdown', itau_sechiba, swdown, kjpindex, kindex) |
---|
| 764 | CALL histwrite_p (hist2_id, 'tair', itau_sechiba, temp_air, kjpindex, kindex) |
---|
| 765 | CALL histwrite_p (hist2_id, 'qair', itau_sechiba, qair, kjpindex, kindex) |
---|
| 766 | ! |
---|
| 767 | CALL histwrite_p (hist2_id, 'evap', itau_sechiba, vevapp, kjpindex, kindex) |
---|
| 768 | CALL histwrite_p (hist2_id, 'coastalflow',itau_sechiba, coastalflow, kjpindex, kindex) |
---|
| 769 | CALL histwrite_p (hist2_id, 'riverflow',itau_sechiba, riverflow, kjpindex, kindex) |
---|
| 770 | ! |
---|
| 771 | CALL histwrite_p (hist2_id, 'temp_sol', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 772 | CALL histwrite_p (hist2_id, 'tsol_max', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 773 | CALL histwrite_p (hist2_id, 'tsol_min', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 774 | CALL histwrite_p (hist2_id, 'fluxsens', itau_sechiba, fluxsens, kjpindex, kindex) |
---|
| 775 | CALL histwrite_p (hist2_id, 'fluxlat', itau_sechiba, fluxlat, kjpindex, kindex) |
---|
| 776 | CALL histwrite_p (hist2_id, 'swnet', itau_sechiba, swnet, kjpindex, kindex) |
---|
| 777 | ! |
---|
| 778 | CALL histwrite_p (hist2_id, 'alb_vis', itau_sechiba, albedo(:,1), kjpindex, kindex) |
---|
| 779 | CALL histwrite_p (hist2_id, 'alb_nir', itau_sechiba, albedo(:,2), kjpindex, kindex) |
---|
| 780 | ENDIF |
---|
| 781 | ELSE |
---|
| 782 | ! |
---|
| 783 | ! Input variables |
---|
| 784 | ! |
---|
| 785 | CALL histwrite_p (hist_id, 'SinAng', itau_sechiba, sinang, kjpindex, kindex) |
---|
| 786 | CALL histwrite_p (hist_id, 'LWdown', itau_sechiba, lwdown, kjpindex, kindex) |
---|
| 787 | CALL histwrite_p (hist_id, 'SWdown', itau_sechiba, swdown, kjpindex, kindex) |
---|
| 788 | CALL histwrite_p (hist_id, 'Tair', itau_sechiba, temp_air, kjpindex, kindex) |
---|
| 789 | CALL histwrite_p (hist_id, 'Qair', itau_sechiba, qair, kjpindex, kindex) |
---|
| 790 | CALL histwrite_p (hist_id, 'SurfP', itau_sechiba, pb, kjpindex, kindex) |
---|
| 791 | CALL histwrite_p (hist_id, 'Windu', itau_sechiba, u_tq, kjpindex, kindex) |
---|
| 792 | CALL histwrite_p (hist_id, 'Windv', itau_sechiba, v_tq, kjpindex, kindex) |
---|
| 793 | ! |
---|
| 794 | CALL histwrite_p (hist_id, 'Evap', itau_sechiba, vevapp, kjpindex, kindex) |
---|
| 795 | CALL histwrite_p (hist_id, 'SWnet', itau_sechiba, swnet, kjpindex, kindex) |
---|
| 796 | CALL histwrite_p (hist_id, 'Qh', itau_sechiba, fluxsens, kjpindex, kindex) |
---|
| 797 | CALL histwrite_p (hist_id, 'Qle', itau_sechiba, fluxlat, kjpindex, kindex) |
---|
| 798 | CALL histwrite_p (hist_id, 'AvgSurfT', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 799 | CALL histwrite_p (hist_id, 'RadT', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 800 | ! |
---|
| 801 | ! There is a mess with the units passed to the coupler. To be checked with Marc |
---|
| 802 | ! |
---|
| 803 | IF ( river_routing ) THEN |
---|
| 804 | CALL histwrite_p (hist_id, 'CoastalFlow',itau_sechiba, coastalflow, kjpindex, kindex) |
---|
| 805 | CALL histwrite_p (hist_id, 'RiverFlow',itau_sechiba, riverflow, kjpindex, kindex) |
---|
| 806 | ENDIF |
---|
| 807 | ! |
---|
| 808 | IF ( hist2_id > 0 ) THEN |
---|
| 809 | CALL histwrite_p (hist2_id, 'Evap', itau_sechiba, vevapp, kjpindex, kindex) |
---|
| 810 | CALL histwrite_p (hist2_id, 'SWnet', itau_sechiba, swnet, kjpindex, kindex) |
---|
| 811 | CALL histwrite_p (hist2_id, 'Qh', itau_sechiba, fluxsens, kjpindex, kindex) |
---|
| 812 | CALL histwrite_p (hist2_id, 'Qle', itau_sechiba, fluxlat, kjpindex, kindex) |
---|
| 813 | CALL histwrite_p (hist2_id, 'AvgSurfT', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 814 | CALL histwrite_p (hist2_id, 'RadT', itau_sechiba, temp_sol_new, kjpindex, kindex) |
---|
| 815 | ENDIF |
---|
| 816 | ENDIF |
---|
| 817 | ! |
---|
| 818 | ! |
---|
| 819 | ! |
---|
| 820 | ENDDO |
---|
| 821 | ! |
---|
| 822 | !- |
---|
| 823 | !--------------------------------------------------------------------------------------- |
---|
| 824 | !- |
---|
| 825 | !- Get time and close IOIPSL, OASIS and MPI |
---|
| 826 | !- |
---|
| 827 | !--------------------------------------------------------------------------------------- |
---|
| 828 | !- |
---|
| 829 | ! |
---|
| 830 | CALL histclo |
---|
| 831 | IF(is_root_prc) THEN |
---|
| 832 | CALL restclo |
---|
| 833 | CALL getin_dump |
---|
| 834 | !- |
---|
| 835 | DEALLOCATE(maskinv_glo) |
---|
| 836 | ! |
---|
| 837 | ENDIF |
---|
| 838 | !- |
---|
| 839 | !- Deallocate all variables and reset init flags |
---|
| 840 | !- |
---|
| 841 | CALL orchideeoasis_clear() |
---|
| 842 | ! |
---|
| 843 | WRITE(numout,*) "MPI finalized" |
---|
| 844 | !- |
---|
| 845 | IF ( timemeasure ) THEN |
---|
| 846 | WRITE(numout,*) '------> Total CPU Time waiting to get forcing : ',waitget_cputime |
---|
| 847 | WRITE(numout,*) '------> Total Real Time waiting to get forcing : ',waitget_walltime |
---|
| 848 | WRITE(numout,*) '------> Total CPU Time for ORCHIDEE : ', orchidee_cputime |
---|
| 849 | WRITE(numout,*) '------> Total Real Time for ORCHIDEE : ', orchidee_walltime |
---|
| 850 | WRITE(numout,*) '------> Total CPU Time waiting to put fluxes : ',waitput_cputime |
---|
| 851 | WRITE(numout,*) '------> Total Real Time waiting to put fluxes : ',waitput_walltime |
---|
| 852 | WRITE(numout,*) '------> Total without MPI : CPU Time : ', Get_cpu_Time(timer_mpi) |
---|
| 853 | WRITE(numout,*) '------> Total without MPI : Real Time : ', Get_real_Time(timer_mpi) |
---|
| 854 | CALL stop_timer(timer_global) |
---|
| 855 | CALL stop_timer(timer_mpi) |
---|
| 856 | ENDIF |
---|
| 857 | !- |
---|
| 858 | CALL oasis_terminate(ierror) |
---|
| 859 | ! |
---|
| 860 | WRITE(numout,*) "OASIS terminated" |
---|
| 861 | !- |
---|
| 862 | !- |
---|
| 863 | CONTAINS |
---|
| 864 | |
---|
| 865 | !! ================================================================================================================================ |
---|
| 866 | !! SUBROUTINE : orchideeoasis_clear |
---|
| 867 | !! |
---|
| 868 | !>\BRIEF Clear orchideeoasis |
---|
| 869 | !! |
---|
| 870 | !! DESCRIPTION : Deallocate memory and reset initialization variables to there original values |
---|
| 871 | !! |
---|
| 872 | !_ ================================================================================================================================ |
---|
| 873 | |
---|
| 874 | SUBROUTINE orchideeoasis_clear |
---|
| 875 | |
---|
| 876 | !- Deallocate all variables existing on all procs (list still incomplete) |
---|
| 877 | |
---|
| 878 | IF ( ALLOCATED(lon_glo) ) DEALLOCATE(lon_glo) |
---|
| 879 | IF ( ALLOCATED(lat_glo) ) DEALLOCATE(lat_glo) |
---|
| 880 | IF ( ALLOCATED(mask_glo) ) DEALLOCATE(mask_glo) |
---|
| 881 | IF ( ALLOCATED(area_glo) ) DEALLOCATE(area_glo) |
---|
| 882 | IF ( ALLOCATED(corners_glo) ) DEALLOCATE(corners_glo) |
---|
| 883 | IF ( ALLOCATED(kindex_g) ) DEALLOCATE(kindex_g) |
---|
| 884 | IF ( ALLOCATED(contfrac_glo) ) DEALLOCATE(contfrac_glo) |
---|
| 885 | IF ( ALLOCATED(lalo_glo) ) DEALLOCATE(lalo_glo) |
---|
| 886 | IF ( ALLOCATED(lon) ) DEALLOCATE(lon) |
---|
| 887 | IF ( ALLOCATED(lat) ) DEALLOCATE(lat) |
---|
| 888 | IF ( ALLOCATED(kindex) ) DEALLOCATE(kindex) |
---|
| 889 | IF ( ALLOCATED(lalo_loc) ) DEALLOCATE(lalo_loc) |
---|
| 890 | IF ( ALLOCATED(zlev_tq) ) DEALLOCATE(zlev_tq) |
---|
| 891 | IF ( ALLOCATED(zlev_uv) ) DEALLOCATE(zlev_uv) |
---|
| 892 | IF ( ALLOCATED(u) ) DEALLOCATE(u) |
---|
| 893 | IF ( ALLOCATED(v) ) DEALLOCATE(v) |
---|
| 894 | IF ( ALLOCATED(pb) ) DEALLOCATE(pb) |
---|
| 895 | IF ( ALLOCATED(temp_air) ) DEALLOCATE(temp_air) |
---|
| 896 | IF ( ALLOCATED(qair) ) DEALLOCATE(qair) |
---|
| 897 | IF ( ALLOCATED(precip_rain) ) DEALLOCATE(precip_rain) |
---|
| 898 | IF ( ALLOCATED(precip_snow) ) DEALLOCATE(precip_snow) |
---|
| 899 | IF ( ALLOCATED(swdown) ) DEALLOCATE(swdown) |
---|
| 900 | IF ( ALLOCATED(swnet) ) DEALLOCATE(swnet) |
---|
| 901 | IF ( ALLOCATED(lwdown) ) DEALLOCATE(lwdown) |
---|
| 902 | IF ( ALLOCATED(sinang) ) DEALLOCATE(sinang) |
---|
| 903 | IF ( ALLOCATED(epot_air) ) DEALLOCATE(epot_air) |
---|
| 904 | IF ( ALLOCATED(ccanopy) ) DEALLOCATE(ccanopy) |
---|
| 905 | IF ( ALLOCATED(cdrag) ) DEALLOCATE(cdrag) |
---|
| 906 | IF ( ALLOCATED(swnet) ) DEALLOCATE(swnet) |
---|
| 907 | IF ( ALLOCATED(petAcoef) ) DEALLOCATE(petAcoef) |
---|
| 908 | IF ( ALLOCATED(peqAcoef) ) DEALLOCATE(peqAcoef) |
---|
| 909 | IF ( ALLOCATED(petBcoef) ) DEALLOCATE(petBcoef) |
---|
| 910 | IF ( ALLOCATED(peqBcoef) ) DEALLOCATE(peqBcoef) |
---|
| 911 | IF ( ALLOCATED(u_tq) ) DEALLOCATE(u_tq) |
---|
| 912 | IF ( ALLOCATED(v_tq) ) DEALLOCATE(v_tq) |
---|
| 913 | IF ( ALLOCATED(vevapp) ) DEALLOCATE(vevapp) |
---|
| 914 | IF ( ALLOCATED(fluxsens) ) DEALLOCATE(fluxsens) |
---|
| 915 | IF ( ALLOCATED(fluxlat) ) DEALLOCATE(fluxlat) |
---|
| 916 | IF ( ALLOCATED(coastalflow) ) DEALLOCATE(coastalflow) |
---|
| 917 | IF ( ALLOCATED(riverflow) ) DEALLOCATE(riverflow) |
---|
| 918 | IF ( ALLOCATED(netco2) ) DEALLOCATE(netco2) |
---|
| 919 | IF ( ALLOCATED(carblu) ) DEALLOCATE(carblu) |
---|
| 920 | IF ( ALLOCATED(carbwh) ) DEALLOCATE(carbwh) |
---|
| 921 | IF ( ALLOCATED(carbha) ) DEALLOCATE(carbha) |
---|
| 922 | IF ( ALLOCATED(tsol_rad) ) DEALLOCATE(tsol_rad) |
---|
| 923 | IF ( ALLOCATED(temp_sol_new) ) DEALLOCATE(temp_sol_new) |
---|
| 924 | IF ( ALLOCATED(qsurf) ) DEALLOCATE(qsurf) |
---|
| 925 | IF ( ALLOCATED(albedo) ) DEALLOCATE(albedo) |
---|
| 926 | IF ( ALLOCATED(emis) ) DEALLOCATE(emis) |
---|
| 927 | IF ( ALLOCATED(z0) ) DEALLOCATE(z0) |
---|
| 928 | |
---|
| 929 | CALL sechiba_clear() |
---|
| 930 | WRITE(numout,*) "Memory deallocated" |
---|
| 931 | |
---|
| 932 | END SUBROUTINE orchideeoasis_clear |
---|
| 933 | |
---|
| 934 | END PROGRAM orchideeoasis |
---|
| 935 | |
---|