Context Navigation

orchideeoasis.f90 @ 7746

Last change on this file since 7746 was 4646, checked in by josefine.ghattas, 7 years ago

Ticket #185 :
Created new module containing time variables and subroutines to calculate them. Variables in the time module are public and can be used elswhere in the model. They should not be modified outside the time module.

Note that now each time step have a time interval. xx_start or xx_end values for the date can be used. Previously the date values correspondend to the end of the interval. See description in module time for more information.

These changes do not make any difference in results for couled mode with LMDZ or dim2_driver offline driver.

These modifications also corrects the problem with do_slow related to orchideedriver, ticket #327

File size: 40.1 KB

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

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: branches/publications/ORCHIDEE_DFv1.0_site/src_oasisdriver/orchideeoasis.f90 @ 7746

Download in other formats: