Context Navigation

cpl_mod.F90 @ 4104

Last change on this file since 4104 was 1274, checked in by Laurent Fairhead, 15 years ago

possibilite de passer le module du stress du vent (activation ou pas via la namcouple Oasis)
petite correction OpenMP

Property svn:eol-style set to native
Property svn:executable set to *
Property svn:keywords set to Author Date Id Revision

File size: 54.8 KB

Line
1	!
2	MODULE cpl_mod
3	!
4	! This module excahanges and transforms all fields that should be recieved or sent to
5	! coupler. The transformation of the fields are done from the grid 1D-array in phylmd
6	! to the regular 2D grid accepted by the coupler. Cumulation of the fields for each
7	! timestep is done in here.
8	!
9	! Each type of surface that recevie fields from the coupler have a subroutine named
10	! cpl_receive_XXX_fields and each surface that have fields to be sent to the coupler
11	! have a subroutine named cpl_send_XXX_fields.
12	!
13	!*************************************************************************************
14
15	! Use statements
16	!*************************************************************************************
17	USE dimphy, ONLY : klon
18	USE mod_phys_lmdz_para
19	USE ioipsl
20	USE iophy
21
22	! The module oasis is always used. Without the cpp key CPP_COUPLE only the parameters
23	! in the module are compiled and not the subroutines.
24	USE oasis
25	USE write_field_phy
26
27	! Global attributes
28	!*************************************************************************************
29	IMPLICIT NONE
30	PRIVATE
31
32	! All subroutine are public except cpl_send_all
33	PUBLIC :: cpl_init, cpl_receive_frac, cpl_receive_ocean_fields, cpl_receive_seaice_fields, &
34	cpl_send_ocean_fields, cpl_send_seaice_fields, cpl_send_land_fields, &
35	cpl_send_landice_fields, gath2cpl
36
37
38	! Declaration of module variables
39	!*************************************************************************************
40	! variable for coupling period
41	INTEGER, SAVE :: nexca
42	!$OMP THREADPRIVATE(nexca)
43
44	! variables for cumulating fields during a coupling periode :
45	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_sols, cpl_nsol, cpl_rain
46	!$OMP THREADPRIVATE(cpl_sols,cpl_nsol,cpl_rain)
47	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_snow, cpl_evap, cpl_tsol
48	!$OMP THREADPRIVATE(cpl_snow,cpl_evap,cpl_tsol)
49	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_fder, cpl_albe, cpl_taux, cpl_tauy
50	!$OMP THREADPRIVATE(cpl_fder,cpl_albe,cpl_taux,cpl_tauy)
51	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_windsp
52	!$OMP THREADPRIVATE(cpl_windsp)
53	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_taumod
54	!$OMP THREADPRIVATE(cpl_taumod)
55	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_atm_co2
56	!$OMP THREADPRIVATE(cpl_atm_co2)
57	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_rriv2D, cpl_rcoa2D, cpl_rlic2D
58	!$OMP THREADPRIVATE(cpl_rriv2D,cpl_rcoa2D,cpl_rlic2D)
59
60	! variables read from coupler :
61	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: read_sst ! sea surface temperature
62	!$OMP THREADPRIVATE(read_sst)
63	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: read_sit ! sea ice temperature
64	!$OMP THREADPRIVATE(read_sit)
65	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: read_sic ! sea ice fraction
66	!$OMP THREADPRIVATE(read_sic)
67	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: read_alb_sic ! albedo at sea ice
68	!$OMP THREADPRIVATE(read_alb_sic)
69	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: read_u0, read_v0 ! ocean surface current
70	!$OMP THREADPRIVATE(read_u0,read_v0)
71	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: read_co2 ! ocean co2 flux
72	!$OMP THREADPRIVATE(read_co2)
73	INTEGER, ALLOCATABLE, DIMENSION(:), SAVE :: unity
74	!$OMP THREADPRIVATE(unity)
75	INTEGER, SAVE :: nidct, nidcs
76	!$OMP THREADPRIVATE(nidct,nidcs)
77
78	! variables to be sent to the coupler
79	REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: cpl_sols2D, cpl_nsol2D, cpl_rain2D
80	!$OMP THREADPRIVATE(cpl_sols2D, cpl_nsol2D, cpl_rain2D)
81	REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: cpl_snow2D, cpl_evap2D, cpl_tsol2D
82	!$OMP THREADPRIVATE(cpl_snow2D, cpl_evap2D, cpl_tsol2D)
83	REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: cpl_fder2D, cpl_albe2D
84	!$OMP THREADPRIVATE(cpl_fder2D, cpl_albe2D)
85	REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: cpl_taux2D, cpl_tauy2D
86	!$OMP THREADPRIVATE(cpl_taux2D, cpl_tauy2D)
87	REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: cpl_taumod2D
88	!$OMP THREADPRIVATE(cpl_taumod2D)
89	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_windsp2D
90	!$OMP THREADPRIVATE(cpl_windsp2D)
91	REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: cpl_atm_co22D
92	!$OMP THREADPRIVATE(cpl_atm_co22D)
93
94	CONTAINS
95	!
96	!************************************************************************************
97	!
98	SUBROUTINE cpl_init(dtime, rlon, rlat)
99	USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, fco2_ocn_day
100
101	INCLUDE "dimensions.h"
102	INCLUDE "indicesol.h"
103	INCLUDE "control.h"
104	INCLUDE "temps.h"
105	INCLUDE "iniprint.h"
106
107	! Input arguments
108	!*************************************************************************************
109	REAL, INTENT(IN) :: dtime
110	REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat
111
112	! Local variables
113	!*************************************************************************************
114	INTEGER :: error, sum_error, ig, i
115	INTEGER :: jf, nhoridct
116	INTEGER :: nhoridcs
117	INTEGER :: idtime
118	INTEGER :: idayref
119	INTEGER :: npas ! only for OASIS2
120	REAL :: zjulian
121	REAL, DIMENSION(iim,jjm+1) :: zx_lon, zx_lat
122	CHARACTER(len = 20) :: modname = 'cpl_init'
123	CHARACTER(len = 80) :: abort_message
124	CHARACTER(len=80) :: clintocplnam, clfromcplnam
125
126	!*************************************************************************************
127	! Calculate coupling period
128	!
129	!*************************************************************************************
130
131	npas = itaufin/ iphysiq
132	nexca = 86400 / dtime
133	WRITE(lunout,*)' ##### Ocean couple #####'
134	WRITE(lunout,*)' Valeurs des pas de temps'
135	WRITE(lunout,*)' npas = ', npas
136	WRITE(lunout,*)' nexca = ', nexca
137
138	!*************************************************************************************
139	! Allocate variables
140	!
141	!*************************************************************************************
142	error = 0
143	sum_error = 0
144
145	ALLOCATE(unity(klon), stat = error)
146	sum_error = sum_error + error
147	ALLOCATE(cpl_sols(klon,2), stat = error)
148	sum_error = sum_error + error
149	ALLOCATE(cpl_nsol(klon,2), stat = error)
150	sum_error = sum_error + error
151	ALLOCATE(cpl_rain(klon,2), stat = error)
152	sum_error = sum_error + error
153	ALLOCATE(cpl_snow(klon,2), stat = error)
154	sum_error = sum_error + error
155	ALLOCATE(cpl_evap(klon,2), stat = error)
156	sum_error = sum_error + error
157	ALLOCATE(cpl_tsol(klon,2), stat = error)
158	sum_error = sum_error + error
159	ALLOCATE(cpl_fder(klon,2), stat = error)
160	sum_error = sum_error + error
161	ALLOCATE(cpl_albe(klon,2), stat = error)
162	sum_error = sum_error + error
163	ALLOCATE(cpl_taux(klon,2), stat = error)
164	sum_error = sum_error + error
165	ALLOCATE(cpl_tauy(klon,2), stat = error)
166	sum_error = sum_error + error
167	ALLOCATE(cpl_windsp(klon,2), stat = error)
168	sum_error = sum_error + error
169	ALLOCATE(cpl_taumod(klon,2), stat = error)
170	sum_error = sum_error + error
171	ALLOCATE(cpl_rriv2D(iim,jj_nb), stat=error)
172	sum_error = sum_error + error
173	ALLOCATE(cpl_rcoa2D(iim,jj_nb), stat=error)
174	sum_error = sum_error + error
175	ALLOCATE(cpl_rlic2D(iim,jj_nb), stat=error)
176	sum_error = sum_error + error
177	ALLOCATE(read_sst(iim, jj_nb), stat = error)
178	sum_error = sum_error + error
179	ALLOCATE(read_sic(iim, jj_nb), stat = error)
180	sum_error = sum_error + error
181	ALLOCATE(read_sit(iim, jj_nb), stat = error)
182	sum_error = sum_error + error
183	ALLOCATE(read_alb_sic(iim, jj_nb), stat = error)
184	sum_error = sum_error + error
185	ALLOCATE(read_u0(iim, jj_nb), stat = error)
186	sum_error = sum_error + error
187	ALLOCATE(read_v0(iim, jj_nb), stat = error)
188	sum_error = sum_error + error
189
190	IF (carbon_cycle_cpl) THEN
191	ALLOCATE(read_co2(iim, jj_nb), stat = error)
192	sum_error = sum_error + error
193	ALLOCATE(cpl_atm_co2(klon,2), stat = error)
194	sum_error = sum_error + error
195
196	! Allocate variable in carbon_cycle_mod
197	ALLOCATE(fco2_ocn_day(klon), stat = error)
198	sum_error = sum_error + error
199	END IF
200
201	IF (sum_error /= 0) THEN
202	abort_message='Pb allocation variables couplees'
203	CALL abort_gcm(modname,abort_message,1)
204	ENDIF
205	!*************************************************************************************
206	! Initialize the allocated varaibles
207	!
208	!*************************************************************************************
209	DO ig = 1, klon
210	unity(ig) = ig
211	ENDDO
212
213	!*************************************************************************************
214	! Initialize coupling
215	!
216	!*************************************************************************************
217	idtime = INT(dtime)
218	#ifdef CPP_COUPLE
219	CALL inicma
220	#endif
221
222	!*************************************************************************************
223	! initialize NetCDF output
224	!
225	!*************************************************************************************
226	IF (is_sequential) THEN
227	idayref = day_ini
228	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
229	CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
230	DO i = 1, iim
231	zx_lon(i,1) = rlon(i+1)
232	zx_lon(i,jjm+1) = rlon(i+1)
233	ENDDO
234	CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
235	clintocplnam="cpl_atm_tauflx"
236	CALL histbeg(clintocplnam, iim,zx_lon(:,1),jjm+1,zx_lat(1,:),&
237	1,iim,1,jjm+1, itau_phy,zjulian,dtime,nhoridct,nidct)
238	! no vertical axis
239	CALL histdef(nidct, 'tauxe','tauxe', &
240	"-",iim, jjm+1, nhoridct, 1, 1, 1, -99, 32, "inst", dtime,dtime)
241	CALL histdef(nidct, 'tauyn','tauyn', &
242	"-",iim, jjm+1, nhoridct, 1, 1, 1, -99, 32, "inst", dtime,dtime)
243	CALL histdef(nidct, 'tmp_lon','tmp_lon', &
244	"-",iim, jjm+1, nhoridct, 1, 1, 1, -99, 32, "inst", dtime,dtime)
245	CALL histdef(nidct, 'tmp_lat','tmp_lat', &
246	"-",iim, jjm+1, nhoridct, 1, 1, 1, -99, 32, "inst", dtime,dtime)
247	DO jf=1,maxsend
248	IF (infosend(i)%action) THEN
249	CALL histdef(nidct, infosend(i)%name ,infosend(i)%name , &
250	"-",iim, jjm+1, nhoridct, 1, 1, 1, -99, 32, "inst", dtime,dtime)
251	ENDIF
252	END DO
253	CALL histend(nidct)
254	CALL histsync(nidct)
255
256	clfromcplnam="cpl_atm_sst"
257	CALL histbeg(clfromcplnam, iim,zx_lon(:,1),jjm+1,zx_lat(1,:),1,iim,1,jjm+1, &
258	0,zjulian,dtime,nhoridcs,nidcs)
259	! no vertical axis
260	DO jf=1,maxrecv
261	IF (inforecv(i)%action) THEN
262	CALL histdef(nidcs,inforecv(i)%name ,inforecv(i)%name , &
263	"-",iim, jjm+1, nhoridcs, 1, 1, 1, -99, 32, "inst", dtime,dtime)
264	ENDIF
265	END DO
266	CALL histend(nidcs)
267	CALL histsync(nidcs)
268
269	ENDIF ! is_sequential
270
271	END SUBROUTINE cpl_init
272
273	!
274	!*************************************************************************************
275	!
276
277	SUBROUTINE cpl_receive_frac(itime, dtime, pctsrf, is_modified)
278	! This subroutine receives from coupler for both ocean and seaice
279	! 4 fields : read_sst, read_sic, read_sit and read_alb_sic.
280	! The new sea-ice-land-landice fraction is returned. The others fields
281	! are stored in this module.
282	USE surface_data
283	USE phys_state_var_mod, ONLY : rlon, rlat
284	USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
285
286	INCLUDE "indicesol.h"
287	INCLUDE "temps.h"
288	INCLUDE "iniprint.h"
289	INCLUDE "YOMCST.h"
290	INCLUDE "dimensions.h"
291
292	! Arguments
293	!************************************************************************************
294	INTEGER, INTENT(IN) :: itime
295	REAL, INTENT(IN) :: dtime
296	REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: pctsrf
297	LOGICAL, INTENT(OUT) :: is_modified
298
299	! Local variables
300	!************************************************************************************
301	INTEGER :: j, i, time_sec
302	INTEGER :: itau_w
303	INTEGER, DIMENSION(iim*(jjm+1)) :: ndexcs
304	CHARACTER(len = 20) :: modname = 'cpl_receive_frac'
305	CHARACTER(len = 80) :: abort_message
306	REAL, DIMENSION(klon) :: read_sic1D
307	REAL, DIMENSION(iim,jj_nb,maxrecv) :: tab_read_flds
308	REAL, DIMENSION(klon,nbsrf) :: pctsrf_old
309	REAL, DIMENSION(klon_mpi) :: rlon_mpi, rlat_mpi
310	REAL, DIMENSION(iim, jj_nb) :: tmp_lon, tmp_lat
311	REAL, DIMENSION(iim, jj_nb) :: tmp_r0
312
313	!*************************************************************************************
314	! Start calculation
315	! Get fields from coupler
316	!
317	!*************************************************************************************
318
319	is_modified=.FALSE.
320
321	! Check if right moment to receive from coupler
322	IF (MOD(itime, nexca) == 1) THEN
323	is_modified=.TRUE.
324
325	time_sec=(itime-1)*dtime
326	#ifdef CPP_COUPLE
327	!$OMP MASTER
328	CALL fromcpl(time_sec, tab_read_flds)
329	!$OMP END MASTER
330	#endif
331
332	! NetCDF output of received fields
333	IF (is_sequential) THEN
334	ndexcs(:) = 0
335	itau_w = itau_phy + itime
336	DO i = 1, maxrecv
337	IF (inforecv(i)%action) THEN
338	CALL histwrite(nidcs,inforecv(i)%name,itau_w,tab_read_flds(:,:,i),iim*(jjm+1),ndexcs)
339	ENDIF
340	END DO
341	ENDIF
342
343
344	! Save each field in a 2D array.
345	!$OMP MASTER
346	read_sst(:,:) = tab_read_flds(:,:,idr_sisutw) ! Sea surface temperature
347	read_sic(:,:) = tab_read_flds(:,:,idr_icecov) ! Sea ice concentration
348	read_alb_sic(:,:) = tab_read_flds(:,:,idr_icealw) ! Albedo at sea ice
349	read_sit(:,:) = tab_read_flds(:,:,idr_icetem) ! Sea ice temperature
350	!$OMP END MASTER
351
352	IF (cpl_current) THEN
353
354	! Transform the longitudes and latitudes on 2D arrays
355	CALL gather_omp(rlon,rlon_mpi)
356	CALL gather_omp(rlat,rlat_mpi)
357	!$OMP MASTER
358	CALL Grid1DTo2D_mpi(rlon_mpi,tmp_lon)
359	CALL Grid1DTo2D_mpi(rlat_mpi,tmp_lat)
360
361	! Transform the currents from cartesian to spheric coordinates
362	! tmp_r0 should be zero
363	CALL geo2atm(iim, jj_nb, tab_read_flds(:,:,idr_curenx), &
364	tab_read_flds(:,:,idr_cureny), tab_read_flds(:,:,idr_curenz), &
365	tmp_lon, tmp_lat, &
366	read_u0(:,:), read_v0(:,:), tmp_r0(:,:))
367	!$OMP END MASTER
368
369	ELSE
370	read_u0(:,:) = 0.
371	read_v0(:,:) = 0.
372	ENDIF
373
374	IF (carbon_cycle_cpl) THEN
375	!$OMP MASTER
376	read_co2(:,:) = tab_read_flds(:,:,idr_oceco2) ! CO2 flux
377	!$OMP END MASTER
378	ENDIF
379
380	!*************************************************************************************
381	! Transform seaice fraction (read_sic : ocean-seaice mask) into global
382	! fraction (pctsrf : ocean-seaice-land-landice mask)
383	!
384	!*************************************************************************************
385	CALL cpl2gath(read_sic, read_sic1D, klon, unity)
386
387	pctsrf_old(:,:) = pctsrf(:,:)
388	DO i = 1, klon
389	! treatment only of points with ocean and/or seaice
390	! old land-ocean mask can not be changed
391	IF (pctsrf_old(i,is_oce) + pctsrf_old(i,is_sic) > 0.) THEN
392	pctsrf(i,is_sic) = (pctsrf_old(i,is_oce) + pctsrf_old(i,is_sic)) &
393	* read_sic1D(i)
394	pctsrf(i,is_oce) = (pctsrf_old(i,is_oce) + pctsrf_old(i,is_sic)) &
395	- pctsrf(i,is_sic)
396	ENDIF
397	ENDDO
398
399	END IF ! if time to receive
400
401	END SUBROUTINE cpl_receive_frac
402
403	!
404	!*************************************************************************************
405	!
406
407	SUBROUTINE cpl_receive_ocean_fields(knon, knindex, tsurf_new, u0_new, v0_new)
408	!
409	! This routine returns the field for the ocean that has been read from the coupler
410	! (done earlier with cpl_receive_frac). The field is the temperature.
411	! The temperature is transformed into 1D array with valid points from index 1 to knon.
412	!
413	USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, fco2_ocn_day
414	INCLUDE "indicesol.h"
415
416	! Input arguments
417	!*************************************************************************************
418	INTEGER, INTENT(IN) :: knon
419	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
420
421	! Output arguments
422	!*************************************************************************************
423	REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new
424	REAL, DIMENSION(klon), INTENT(OUT) :: u0_new
425	REAL, DIMENSION(klon), INTENT(OUT) :: v0_new
426
427	! Local variables
428	!*************************************************************************************
429	INTEGER :: i
430	INTEGER, DIMENSION(klon) :: index
431	REAL, DIMENSION(klon) :: sic_new
432
433	!*************************************************************************************
434	! Transform read_sst into compressed 1D variable tsurf_new
435	!
436	!*************************************************************************************
437	CALL cpl2gath(read_sst, tsurf_new, knon, knindex)
438	CALL cpl2gath(read_sic, sic_new, knon, knindex)
439	CALL cpl2gath(read_u0, u0_new, knon, knindex)
440	CALL cpl2gath(read_v0, v0_new, knon, knindex)
441
442	!*************************************************************************************
443	! Transform read_co2 into uncompressed 1D variable fco2_ocn_day added directly in
444	! the module carbon_cycle_mod
445	!
446	!*************************************************************************************
447	IF (carbon_cycle_cpl) THEN
448	DO i=1,klon
449	index(i)=i
450	END DO
451	CALL cpl2gath(read_co2, fco2_ocn_day, klon, index)
452	END IF
453
454	!*************************************************************************************
455	! The fields received from the coupler have to be weighted with the fraction of ocean
456	! in relation to the total sea-ice+ocean
457	!
458	!*************************************************************************************
459	DO i=1, knon
460	tsurf_new(i) = tsurf_new(i)/(1. - sic_new(i))
461	END DO
462
463	END SUBROUTINE cpl_receive_ocean_fields
464
465	!
466	!*************************************************************************************
467	!
468
469	SUBROUTINE cpl_receive_seaice_fields(knon, knindex, &
470	tsurf_new, alb_new, u0_new, v0_new)
471	!
472	! This routine returns the fields for the seaice that have been read from the coupler
473	! (done earlier with cpl_receive_frac). These fields are the temperature and
474	! albedo at sea ice surface and fraction of sea ice.
475	! The fields are transformed into 1D arrays with valid points from index 1 to knon.
476	!
477
478	! Input arguments
479	!*************************************************************************************
480	INTEGER, INTENT(IN) :: knon
481	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
482
483	! Output arguments
484	!*************************************************************************************
485	REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new
486	REAL, DIMENSION(klon), INTENT(OUT) :: alb_new
487	REAL, DIMENSION(klon), INTENT(OUT) :: u0_new
488	REAL, DIMENSION(klon), INTENT(OUT) :: v0_new
489
490	! Local variables
491	!*************************************************************************************
492	INTEGER :: i
493	REAL, DIMENSION(klon) :: sic_new
494
495	!*************************************************************************************
496	! Transform fields read from coupler from 2D into compressed 1D variables
497	!
498	!*************************************************************************************
499	CALL cpl2gath(read_sit, tsurf_new, knon, knindex)
500	CALL cpl2gath(read_alb_sic, alb_new, knon, knindex)
501	CALL cpl2gath(read_sic, sic_new, knon, knindex)
502	CALL cpl2gath(read_u0, u0_new, knon, knindex)
503	CALL cpl2gath(read_v0, v0_new, knon, knindex)
504
505	!*************************************************************************************
506	! The fields received from the coupler have to be weighted with the sea-ice
507	! concentration (in relation to the total sea-ice + ocean).
508	!
509	!*************************************************************************************
510	DO i= 1, knon
511	tsurf_new(i) = tsurf_new(i) / sic_new(i)
512	alb_new(i) = alb_new(i) / sic_new(i)
513	END DO
514
515	END SUBROUTINE cpl_receive_seaice_fields
516
517	!
518	!*************************************************************************************
519	!
520
521	SUBROUTINE cpl_send_ocean_fields(itime, knon, knindex, &
522	swdown, lwdown, fluxlat, fluxsens, &
523	precip_rain, precip_snow, evap, tsurf, fder, albsol, taux, tauy, windsp)
524	!
525	! This subroutine cumulates some fields for each time-step during a coupling
526	! period. At last time-step in a coupling period the fields are transformed to the
527	! grid accepted by the coupler. No sending to the coupler will be done from here
528	! (it is done in cpl_send_seaice_fields).
529	!
530	USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, co2_send
531	INCLUDE "indicesol.h"
532	INCLUDE "dimensions.h"
533
534	! Input arguments
535	!*************************************************************************************
536	INTEGER, INTENT(IN) :: itime
537	INTEGER, INTENT(IN) :: knon
538	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
539	REAL, DIMENSION(klon), INTENT(IN) :: swdown, lwdown
540	REAL, DIMENSION(klon), INTENT(IN) :: fluxlat, fluxsens
541	REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow
542	REAL, DIMENSION(klon), INTENT(IN) :: evap, tsurf, fder, albsol
543	REAL, DIMENSION(klon), INTENT(IN) :: taux, tauy, windsp
544
545	! Local variables
546	!*************************************************************************************
547	INTEGER :: cpl_index, ig
548	INTEGER :: error, sum_error
549	CHARACTER(len = 25) :: modname = 'cpl_send_ocean_fields'
550	CHARACTER(len = 80) :: abort_message
551
552	!*************************************************************************************
553	! Start calculation
554	! The ocean points are saved with second array index=1
555	!
556	!*************************************************************************************
557	cpl_index = 1
558
559	!*************************************************************************************
560	! Reset fields to zero in the beginning of a new coupling period
561	!
562	!*************************************************************************************
563	IF (MOD(itime, nexca) == 1) THEN
564	cpl_sols(1:knon,cpl_index) = 0.0
565	cpl_nsol(1:knon,cpl_index) = 0.0
566	cpl_rain(1:knon,cpl_index) = 0.0
567	cpl_snow(1:knon,cpl_index) = 0.0
568	cpl_evap(1:knon,cpl_index) = 0.0
569	cpl_tsol(1:knon,cpl_index) = 0.0
570	cpl_fder(1:knon,cpl_index) = 0.0
571	cpl_albe(1:knon,cpl_index) = 0.0
572	cpl_taux(1:knon,cpl_index) = 0.0
573	cpl_tauy(1:knon,cpl_index) = 0.0
574	cpl_windsp(1:knon,cpl_index) = 0.0
575	cpl_taumod(1:knon,cpl_index) = 0.0
576	IF (carbon_cycle_cpl) cpl_atm_co2(1:knon,cpl_index) = 0.0
577	ENDIF
578
579	!*************************************************************************************
580	! Cumulate at each time-step
581	!
582	!*************************************************************************************
583	DO ig = 1, knon
584	cpl_sols(ig,cpl_index) = cpl_sols(ig,cpl_index) + &
585	swdown(ig) / FLOAT(nexca)
586	cpl_nsol(ig,cpl_index) = cpl_nsol(ig,cpl_index) + &
587	(lwdown(ig) + fluxlat(ig) +fluxsens(ig)) / FLOAT(nexca)
588	cpl_rain(ig,cpl_index) = cpl_rain(ig,cpl_index) + &
589	precip_rain(ig) / FLOAT(nexca)
590	cpl_snow(ig,cpl_index) = cpl_snow(ig,cpl_index) + &
591	precip_snow(ig) / FLOAT(nexca)
592	cpl_evap(ig,cpl_index) = cpl_evap(ig,cpl_index) + &
593	evap(ig) / FLOAT(nexca)
594	cpl_tsol(ig,cpl_index) = cpl_tsol(ig,cpl_index) + &
595	tsurf(ig) / FLOAT(nexca)
596	cpl_fder(ig,cpl_index) = cpl_fder(ig,cpl_index) + &
597	fder(ig) / FLOAT(nexca)
598	cpl_albe(ig,cpl_index) = cpl_albe(ig,cpl_index) + &
599	albsol(ig) / FLOAT(nexca)
600	cpl_taux(ig,cpl_index) = cpl_taux(ig,cpl_index) + &
601	taux(ig) / FLOAT(nexca)
602	cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) + &
603	tauy(ig) / FLOAT(nexca)
604	cpl_windsp(ig,cpl_index) = cpl_windsp(ig,cpl_index) + &
605	windsp(ig) / FLOAT(nexca)
606	cpl_taumod(ig,cpl_index) = cpl_taumod(ig,cpl_index) + &
607	SQRT ( taux(ig)taux(ig)+tauy(ig)tauy(ig) ) / FLOAT (nexca)
608
609	IF (carbon_cycle_cpl) THEN
610	cpl_atm_co2(ig,cpl_index) = cpl_atm_co2(ig,cpl_index) + &
611	co2_send(knindex(ig))/ FLOAT(nexca)
612	END IF
613	ENDDO
614
615	!*************************************************************************************
616	! If the time-step corresponds to the end of coupling period the
617	! fields are transformed to the 2D grid.
618	! No sending to the coupler (it is done from cpl_send_seaice_fields).
619	!
620	!*************************************************************************************
621	IF (MOD(itime, nexca) == 0) THEN
622
623	IF (.NOT. ALLOCATED(cpl_sols2D)) THEN
624	sum_error = 0
625	ALLOCATE(cpl_sols2D(iim,jj_nb,2), stat=error)
626	sum_error = sum_error + error
627	ALLOCATE(cpl_nsol2D(iim,jj_nb,2), stat=error)
628	sum_error = sum_error + error
629	ALLOCATE(cpl_rain2D(iim,jj_nb,2), stat=error)
630	sum_error = sum_error + error
631	ALLOCATE(cpl_snow2D(iim,jj_nb,2), stat=error)
632	sum_error = sum_error + error
633	ALLOCATE(cpl_evap2D(iim,jj_nb,2), stat=error)
634	sum_error = sum_error + error
635	ALLOCATE(cpl_tsol2D(iim,jj_nb,2), stat=error)
636	sum_error = sum_error + error
637	ALLOCATE(cpl_fder2D(iim,jj_nb,2), stat=error)
638	sum_error = sum_error + error
639	ALLOCATE(cpl_albe2D(iim,jj_nb,2), stat=error)
640	sum_error = sum_error + error
641	ALLOCATE(cpl_taux2D(iim,jj_nb,2), stat=error)
642	sum_error = sum_error + error
643	ALLOCATE(cpl_tauy2D(iim,jj_nb,2), stat=error)
644	sum_error = sum_error + error
645	ALLOCATE(cpl_windsp2D(iim,jj_nb), stat=error)
646	sum_error = sum_error + error
647	ALLOCATE(cpl_taumod2D(iim,jj_nb,2), stat=error)
648	sum_error = sum_error + error
649
650	IF (carbon_cycle_cpl) THEN
651	ALLOCATE(cpl_atm_co22D(iim,jj_nb), stat=error)
652	sum_error = sum_error + error
653	END IF
654
655	IF (sum_error /= 0) THEN
656	abort_message='Pb allocation variables couplees pour l''ecriture'
657	CALL abort_gcm(modname,abort_message,1)
658	ENDIF
659	ENDIF
660
661
662	CALL gath2cpl(cpl_sols(:,cpl_index), cpl_sols2D(:,:,cpl_index), &
663	knon, knindex)
664
665	CALL gath2cpl(cpl_nsol(:,cpl_index), cpl_nsol2D(:,:,cpl_index), &
666	knon, knindex)
667
668	CALL gath2cpl(cpl_rain(:,cpl_index), cpl_rain2D(:,:,cpl_index), &
669	knon, knindex)
670
671	CALL gath2cpl(cpl_snow(:,cpl_index), cpl_snow2D(:,:,cpl_index), &
672	knon, knindex)
673
674	CALL gath2cpl(cpl_evap(:,cpl_index), cpl_evap2D(:,:,cpl_index), &
675	knon, knindex)
676
677	! cpl_tsol2D(:,:,:) not used!
678	CALL gath2cpl(cpl_tsol(:,cpl_index), cpl_tsol2D(:,:, cpl_index), &
679	knon, knindex)
680
681	! cpl_fder2D(:,:,1) not used, only cpl_fder(:,:,2)!
682	CALL gath2cpl(cpl_fder(:,cpl_index), cpl_fder2D(:,:,cpl_index), &
683	knon, knindex)
684
685	! cpl_albe2D(:,:,:) not used!
686	CALL gath2cpl(cpl_albe(:,cpl_index), cpl_albe2D(:,:,cpl_index), &
687	knon, knindex)
688
689	CALL gath2cpl(cpl_taux(:,cpl_index), cpl_taux2D(:,:,cpl_index), &
690	knon, knindex)
691
692	CALL gath2cpl(cpl_tauy(:,cpl_index), cpl_tauy2D(:,:,cpl_index), &
693	knon, knindex)
694
695	CALL gath2cpl(cpl_windsp(:,cpl_index), cpl_windsp2D(:,:), &
696	knon, knindex)
697
698	CALL gath2cpl(cpl_taumod(:,cpl_index), cpl_taumod2D(:,:,cpl_index), &
699	knon, knindex)
700
701	IF (carbon_cycle_cpl) &
702	CALL gath2cpl(cpl_atm_co2(:,cpl_index), cpl_atm_co22D(:,:), knon, knindex)
703	ENDIF
704
705	END SUBROUTINE cpl_send_ocean_fields
706
707	!
708	!*************************************************************************************
709	!
710
711	SUBROUTINE cpl_send_seaice_fields(itime, dtime, knon, knindex, &
712	pctsrf, lafin, rlon, rlat, &
713	swdown, lwdown, fluxlat, fluxsens, &
714	precip_rain, precip_snow, evap, tsurf, fder, albsol, taux, tauy)
715	!
716	! This subroutine cumulates some fields for each time-step during a coupling
717	! period. At last time-step in a coupling period the fields are transformed to the
718	! grid accepted by the coupler. All fields for all types of surfaces are sent to
719	! the coupler.
720	!
721	USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
722	INCLUDE "indicesol.h"
723	INCLUDE "dimensions.h"
724
725	! Input arguments
726	!*************************************************************************************
727	INTEGER, INTENT(IN) :: itime
728	INTEGER, INTENT(IN) :: knon
729	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
730	REAL, INTENT(IN) :: dtime
731	REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat
732	REAL, DIMENSION(klon), INTENT(IN) :: swdown, lwdown
733	REAL, DIMENSION(klon), INTENT(IN) :: fluxlat, fluxsens
734	REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow
735	REAL, DIMENSION(klon), INTENT(IN) :: evap, tsurf, fder
736	REAL, DIMENSION(klon), INTENT(IN) :: albsol, taux, tauy
737	REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf
738	LOGICAL, INTENT(IN) :: lafin
739
740	! Local variables
741	!*************************************************************************************
742	INTEGER :: cpl_index, ig
743	INTEGER :: error, sum_error
744	CHARACTER(len = 25) :: modname = 'cpl_send_seaice_fields'
745	CHARACTER(len = 80) :: abort_message
746	REAL, DIMENSION(klon) :: cpl_fder_tmp
747
748	!*************************************************************************************
749	! Start calulation
750	! The sea-ice points are saved with second array index=2
751	!
752	!*************************************************************************************
753	cpl_index = 2
754
755	!*************************************************************************************
756	! Reset fields to zero in the beginning of a new coupling period
757	!
758	!*************************************************************************************
759	IF (MOD(itime, nexca) == 1) THEN
760	cpl_sols(1:knon,cpl_index) = 0.0
761	cpl_nsol(1:knon,cpl_index) = 0.0
762	cpl_rain(1:knon,cpl_index) = 0.0
763	cpl_snow(1:knon,cpl_index) = 0.0
764	cpl_evap(1:knon,cpl_index) = 0.0
765	cpl_tsol(1:knon,cpl_index) = 0.0
766	cpl_fder(1:knon,cpl_index) = 0.0
767	cpl_albe(1:knon,cpl_index) = 0.0
768	cpl_taux(1:knon,cpl_index) = 0.0
769	cpl_tauy(1:knon,cpl_index) = 0.0
770	cpl_taumod(1:knon,cpl_index) = 0.0
771	ENDIF
772
773	!*************************************************************************************
774	! Cumulate at each time-step
775	!
776	!*************************************************************************************
777	DO ig = 1, knon
778	cpl_sols(ig,cpl_index) = cpl_sols(ig,cpl_index) + &
779	swdown(ig) / FLOAT(nexca)
780	cpl_nsol(ig,cpl_index) = cpl_nsol(ig,cpl_index) + &
781	(lwdown(ig) + fluxlat(ig) +fluxsens(ig)) / FLOAT(nexca)
782	cpl_rain(ig,cpl_index) = cpl_rain(ig,cpl_index) + &
783	precip_rain(ig) / FLOAT(nexca)
784	cpl_snow(ig,cpl_index) = cpl_snow(ig,cpl_index) + &
785	precip_snow(ig) / FLOAT(nexca)
786	cpl_evap(ig,cpl_index) = cpl_evap(ig,cpl_index) + &
787	evap(ig) / FLOAT(nexca)
788	cpl_tsol(ig,cpl_index) = cpl_tsol(ig,cpl_index) + &
789	tsurf(ig) / FLOAT(nexca)
790	cpl_fder(ig,cpl_index) = cpl_fder(ig,cpl_index) + &
791	fder(ig) / FLOAT(nexca)
792	cpl_albe(ig,cpl_index) = cpl_albe(ig,cpl_index) + &
793	albsol(ig) / FLOAT(nexca)
794	cpl_taux(ig,cpl_index) = cpl_taux(ig,cpl_index) + &
795	taux(ig) / FLOAT(nexca)
796	cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) + &
797	tauy(ig) / FLOAT(nexca)
798	cpl_taumod(ig,cpl_index) = cpl_taumod(ig,cpl_index) + &
799	SQRT ( taux(ig)taux(ig)+tauy(ig)tauy(ig) ) / FLOAT(nexca)
800	ENDDO
801
802	!*************************************************************************************
803	! If the time-step corresponds to the end of coupling period the
804	! fields are transformed to the 2D grid and all fields are sent to coupler.
805	!
806	!*************************************************************************************
807	IF (MOD(itime, nexca) == 0) THEN
808	IF (.NOT. ALLOCATED(cpl_sols2D)) THEN
809	sum_error = 0
810	ALLOCATE(cpl_sols2D(iim,jj_nb,2), stat=error)
811	sum_error = sum_error + error
812	ALLOCATE(cpl_nsol2D(iim,jj_nb,2), stat=error)
813	sum_error = sum_error + error
814	ALLOCATE(cpl_rain2D(iim,jj_nb,2), stat=error)
815	sum_error = sum_error + error
816	ALLOCATE(cpl_snow2D(iim,jj_nb,2), stat=error)
817	sum_error = sum_error + error
818	ALLOCATE(cpl_evap2D(iim,jj_nb,2), stat=error)
819	sum_error = sum_error + error
820	ALLOCATE(cpl_tsol2D(iim,jj_nb,2), stat=error)
821	sum_error = sum_error + error
822	ALLOCATE(cpl_fder2D(iim,jj_nb,2), stat=error)
823	sum_error = sum_error + error
824	ALLOCATE(cpl_albe2D(iim,jj_nb,2), stat=error)
825	sum_error = sum_error + error
826	ALLOCATE(cpl_taux2D(iim,jj_nb,2), stat=error)
827	sum_error = sum_error + error
828	ALLOCATE(cpl_tauy2D(iim,jj_nb,2), stat=error)
829	sum_error = sum_error + error
830	ALLOCATE(cpl_windsp2D(iim,jj_nb), stat=error)
831	sum_error = sum_error + error
832	ALLOCATE(cpl_taumod2D(iim,jj_nb,2), stat=error)
833	sum_error = sum_error + error
834
835	IF (carbon_cycle_cpl) THEN
836	ALLOCATE(cpl_atm_co22D(iim,jj_nb), stat=error)
837	sum_error = sum_error + error
838	END IF
839
840	IF (sum_error /= 0) THEN
841	abort_message='Pb allocation variables couplees pour l''ecriture'
842	CALL abort_gcm(modname,abort_message,1)
843	ENDIF
844	ENDIF
845
846	CALL gath2cpl(cpl_sols(:,cpl_index), cpl_sols2D(:,:,cpl_index), &
847	knon, knindex)
848
849	CALL gath2cpl(cpl_nsol(:,cpl_index), cpl_nsol2D(:,:,cpl_index), &
850	knon, knindex)
851
852	CALL gath2cpl(cpl_rain(:,cpl_index), cpl_rain2D(:,:,cpl_index), &
853	knon, knindex)
854
855	CALL gath2cpl(cpl_snow(:,cpl_index), cpl_snow2D(:,:,cpl_index), &
856	knon, knindex)
857
858	CALL gath2cpl(cpl_evap(:,cpl_index), cpl_evap2D(:,:,cpl_index), &
859	knon, knindex)
860
861	! cpl_tsol2D(:,:,:) not used!
862	CALL gath2cpl(cpl_tsol(:,cpl_index), cpl_tsol2D(:,:, cpl_index), &
863	knon, knindex)
864
865	! Set default value and decompress before gath2cpl
866	cpl_fder_tmp(:) = -20.
867	DO ig = 1, knon
868	cpl_fder_tmp(knindex(ig))=cpl_fder(ig,cpl_index)
869	END DO
870	CALL gath2cpl(cpl_fder_tmp(:), cpl_fder2D(:,:,cpl_index), &
871	klon, unity)
872
873	! cpl_albe2D(:,:,:) not used!
874	CALL gath2cpl(cpl_albe(:,cpl_index), cpl_albe2D(:,:,cpl_index), &
875	knon, knindex)
876
877	CALL gath2cpl(cpl_taux(:,cpl_index), cpl_taux2D(:,:,cpl_index), &
878	knon, knindex)
879
880	CALL gath2cpl(cpl_tauy(:,cpl_index), cpl_tauy2D(:,:,cpl_index), &
881	knon, knindex)
882
883	CALL gath2cpl(cpl_taumod(:,cpl_index), cpl_taumod2D(:,:,cpl_index), &
884	knon, knindex)
885
886	! Send all fields
887	CALL cpl_send_all(itime, dtime, pctsrf, lafin, rlon, rlat)
888	ENDIF
889
890	END SUBROUTINE cpl_send_seaice_fields
891
892	!
893	!*************************************************************************************
894	!
895
896	SUBROUTINE cpl_send_land_fields(itime, knon, knindex, rriv_in, rcoa_in)
897	!
898	! This subroutine cumulates some fields for each time-step during a coupling
899	! period. At last time-step in a coupling period the fields are transformed to the
900	! grid accepted by the coupler. No sending to the coupler will be done from here
901	! (it is done in cpl_send_seaice_fields).
902	!
903	INCLUDE "dimensions.h"
904
905	! Input arguments
906	!*************************************************************************************
907	INTEGER, INTENT(IN) :: itime
908	INTEGER, INTENT(IN) :: knon
909	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
910	REAL, DIMENSION(klon), INTENT(IN) :: rriv_in
911	REAL, DIMENSION(klon), INTENT(IN) :: rcoa_in
912
913	! Local variables
914	!*************************************************************************************
915	REAL, DIMENSION(iim,jj_nb) :: rriv2D
916	REAL, DIMENSION(iim,jj_nb) :: rcoa2D
917
918	!*************************************************************************************
919	! Rearrange fields in 2D variables
920	! First initialize to zero to avoid unvalid points causing problems
921	!
922	!*************************************************************************************
923	!$OMP MASTER
924	rriv2D(:,:) = 0.0
925	rcoa2D(:,:) = 0.0
926	!$OMP END MASTER
927	CALL gath2cpl(rriv_in, rriv2D, knon, knindex)
928	CALL gath2cpl(rcoa_in, rcoa2D, knon, knindex)
929
930	!*************************************************************************************
931	! Reset cumulated fields to zero in the beginning of a new coupling period
932	!
933	!*************************************************************************************
934	IF (MOD(itime, nexca) == 1) THEN
935	!$OMP MASTER
936	cpl_rriv2D(:,:) = 0.0
937	cpl_rcoa2D(:,:) = 0.0
938	!$OMP END MASTER
939	ENDIF
940
941	!*************************************************************************************
942	! Cumulate : Following fields should be cumulated at each time-step
943	!
944	!*************************************************************************************
945	!$OMP MASTER
946	cpl_rriv2D(:,:) = cpl_rriv2D(:,:) + rriv2D(:,:) / FLOAT(nexca)
947	cpl_rcoa2D(:,:) = cpl_rcoa2D(:,:) + rcoa2D(:,:) / FLOAT(nexca)
948	!$OMP END MASTER
949
950	END SUBROUTINE cpl_send_land_fields
951
952	!
953	!*************************************************************************************
954	!
955
956	SUBROUTINE cpl_send_landice_fields(itime, knon, knindex, rlic_in)
957	! This subroutine cumulates the field for melting ice for each time-step
958	! during a coupling period. This routine will not send to coupler. Sending
959	! will be done in cpl_send_seaice_fields.
960	!
961
962	INCLUDE "dimensions.h"
963
964	! Input varibales
965	!*************************************************************************************
966	INTEGER, INTENT(IN) :: itime
967	INTEGER, INTENT(IN) :: knon
968	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
969	REAL, DIMENSION(klon), INTENT(IN) :: rlic_in
970
971	! Local varibales
972	!*************************************************************************************
973	REAL, DIMENSION(iim,jj_nb) :: rlic2D
974
975	!*************************************************************************************
976	! Rearrange field in a 2D variable
977	! First initialize to zero to avoid unvalid points causing problems
978	!
979	!*************************************************************************************
980	!$OMP MASTER
981	rlic2D(:,:) = 0.0
982	!$OMP END MASTER
983	CALL gath2cpl(rlic_in, rlic2D, knon, knindex)
984
985	!*************************************************************************************
986	! Reset field to zero in the beginning of a new coupling period
987	!
988	!*************************************************************************************
989	IF (MOD(itime, nexca) == 1) THEN
990	!$OMP MASTER
991	cpl_rlic2D(:,:) = 0.0
992	!$OMP END MASTER
993	ENDIF
994
995	!*************************************************************************************
996	! Cumulate : Melting ice should be cumulated at each time-step
997	!
998	!*************************************************************************************
999	!$OMP MASTER
1000	cpl_rlic2D(:,:) = cpl_rlic2D(:,:) + rlic2D(:,:) / FLOAT(nexca)
1001	!$OMP END MASTER
1002
1003	END SUBROUTINE cpl_send_landice_fields
1004
1005	!
1006	!*************************************************************************************
1007	!
1008
1009	SUBROUTINE cpl_send_all(itime, dtime, pctsrf, lafin, rlon, rlat)
1010	! This routine will send fields for all different surfaces to the coupler.
1011	! This subroutine should be executed after calculations by the last surface(sea-ice),
1012	! all calculations at the different surfaces have to be done before.
1013	!
1014	USE surface_data
1015	USE carbon_cycle_mod, ONLY : carbon_cycle_cpl
1016	! Some includes
1017	!*************************************************************************************
1018	INCLUDE "indicesol.h"
1019	INCLUDE "temps.h"
1020	INCLUDE "dimensions.h"
1021
1022	! Input arguments
1023	!*************************************************************************************
1024	INTEGER, INTENT(IN) :: itime
1025	REAL, INTENT(IN) :: dtime
1026	REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat
1027	REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf
1028	LOGICAL, INTENT(IN) :: lafin
1029
1030	! Local variables
1031	!*************************************************************************************
1032	INTEGER :: error, sum_error, j
1033	INTEGER :: itau_w
1034	INTEGER :: time_sec
1035	INTEGER, DIMENSION(iim*(jjm+1)) :: ndexct
1036	REAL :: Up, Down
1037	REAL, DIMENSION(iim, jj_nb) :: tmp_lon, tmp_lat
1038	REAL, DIMENSION(iim, jj_nb, 4) :: pctsrf2D
1039	REAL, DIMENSION(iim, jj_nb) :: deno
1040	CHARACTER(len = 20) :: modname = 'cpl_send_all'
1041	CHARACTER(len = 80) :: abort_message
1042
1043	! Variables with fields to coupler
1044	REAL, DIMENSION(iim, jj_nb) :: tmp_taux
1045	REAL, DIMENSION(iim, jj_nb) :: tmp_tauy
1046	REAL, DIMENSION(iim, jj_nb) :: tmp_calv
1047	! Table with all fields to send to coupler
1048	REAL, DIMENSION(iim, jj_nb, maxsend) :: tab_flds
1049	REAL, DIMENSION(klon_mpi) :: rlon_mpi, rlat_mpi
1050
1051	#ifdef CPP_MPI
1052	INCLUDE 'mpif.h'
1053	INTEGER, DIMENSION(MPI_STATUS_SIZE) :: status
1054	#endif
1055
1056	! End definitions
1057	!*************************************************************************************
1058
1059
1060
1061	!*************************************************************************************
1062	! All fields are stored in a table tab_flds(:,:,:)
1063	! First store the fields which are already on the right format
1064	!
1065	!*************************************************************************************
1066	!$OMP MASTER
1067	tab_flds(:,:,ids_windsp) = cpl_windsp2D(:,:)
1068	tab_flds(:,:,ids_shfice) = cpl_sols2D(:,:,2)
1069	tab_flds(:,:,ids_nsfice) = cpl_nsol2D(:,:,2)
1070	tab_flds(:,:,ids_dflxdt) = cpl_fder2D(:,:,2)
1071
1072	IF (version_ocean=='nemo') THEN
1073	tab_flds(:,:,ids_liqrun) = cpl_rriv2D(:,:) + cpl_rcoa2D(:,:)
1074	IF (carbon_cycle_cpl) tab_flds(:,:,ids_atmco2)=cpl_atm_co22D(:,:)
1075	ELSE IF (version_ocean=='opa8') THEN
1076	tab_flds(:,:,ids_shfoce) = cpl_sols2D(:,:,1)
1077	tab_flds(:,:,ids_nsfoce) = cpl_nsol2D(:,:,1)
1078	tab_flds(:,:,ids_icevap) = cpl_evap2D(:,:,2)
1079	tab_flds(:,:,ids_ocevap) = cpl_evap2D(:,:,1)
1080	tab_flds(:,:,ids_runcoa) = cpl_rcoa2D(:,:)
1081	tab_flds(:,:,ids_rivflu) = cpl_rriv2D(:,:)
1082	END IF
1083
1084	!*************************************************************************************
1085	! Transform the fraction of sub-surfaces from 1D to 2D array
1086	!
1087	!*************************************************************************************
1088	pctsrf2D(:,:,:) = 0.
1089	!$OMP END MASTER
1090	CALL gath2cpl(pctsrf(:,is_oce), pctsrf2D(:,:,is_oce), klon, unity)
1091	CALL gath2cpl(pctsrf(:,is_sic), pctsrf2D(:,:,is_sic), klon, unity)
1092	CALL gath2cpl(pctsrf(:,is_lic), pctsrf2D(:,:,is_lic), klon, unity)
1093
1094	!*************************************************************************************
1095	! Calculate the average calving per latitude
1096	! Store calving in tab_flds(:,:,19)
1097	!
1098	!*************************************************************************************
1099	IF (is_omp_root) THEN
1100
1101	DO j = 1, jj_nb
1102	tmp_calv(:,j) = DOT_PRODUCT (cpl_rlic2D(1:iim,j), &
1103	pctsrf2D(1:iim,j,is_lic)) / REAL(iim)
1104	ENDDO
1105
1106
1107	IF (is_parallel) THEN
1108	IF (.NOT. is_north_pole) THEN
1109	#ifdef CPP_MPI
1110	CALL MPI_RECV(Up,1,MPI_REAL_LMDZ,mpi_rank-1,1234,COMM_LMDZ_PHY,status,error)
1111	CALL MPI_SEND(tmp_calv(1,1),1,MPI_REAL_LMDZ,mpi_rank-1,1234,COMM_LMDZ_PHY,error)
1112	#endif
1113	ENDIF
1114
1115	IF (.NOT. is_south_pole) THEN
1116	#ifdef CPP_MPI
1117	CALL MPI_SEND(tmp_calv(1,jj_nb),1,MPI_REAL_LMDZ,mpi_rank+1,1234,COMM_LMDZ_PHY,error)
1118	CALL MPI_RECV(down,1,MPI_REAL_LMDZ,mpi_rank+1,1234,COMM_LMDZ_PHY,status,error)
1119	#endif
1120	ENDIF
1121
1122	IF (.NOT. is_north_pole .AND. ii_begin /=1) THEN
1123	Up=Up+tmp_calv(iim,1)
1124	tmp_calv(:,1)=Up
1125	ENDIF
1126
1127	IF (.NOT. is_south_pole .AND. ii_end /= iim) THEN
1128	Down=Down+tmp_calv(1,jj_nb)
1129	tmp_calv(:,jj_nb)=Down
1130	ENDIF
1131	ENDIF
1132
1133	tab_flds(:,:,ids_calvin) = tmp_calv(:,:)
1134
1135	!*************************************************************************************
1136	! Calculate total flux for snow, rain and wind with weighted addition using the
1137	! fractions of ocean and seaice.
1138	!
1139	!*************************************************************************************
1140	! fraction oce+seaice
1141	deno = pctsrf2D(:,:,is_oce) + pctsrf2D(:,:,is_sic)
1142
1143	IF (version_ocean=='nemo') THEN
1144	tab_flds(:,:,ids_shftot) = 0.0
1145	tab_flds(:,:,ids_nsftot) = 0.0
1146	tab_flds(:,:,ids_totrai) = 0.0
1147	tab_flds(:,:,ids_totsno) = 0.0
1148	tab_flds(:,:,ids_toteva) = 0.0
1149	tab_flds(:,:,ids_taumod) = 0.0
1150
1151	tmp_taux(:,:) = 0.0
1152	tmp_tauy(:,:) = 0.0
1153	! For all valid grid cells containing some fraction of ocean or sea-ice
1154	WHERE ( deno(:,:) /= 0 )
1155	tmp_taux = cpl_taux2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1156	cpl_taux2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1157	tmp_tauy = cpl_tauy2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1158	cpl_tauy2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1159
1160	tab_flds(:,:,ids_shftot) = cpl_sols2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1161	cpl_sols2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1162	tab_flds(:,:,ids_nsftot) = cpl_nsol2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1163	cpl_nsol2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1164	tab_flds(:,:,ids_totrai) = cpl_rain2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1165	cpl_rain2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1166	tab_flds(:,:,ids_totsno) = cpl_snow2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1167	cpl_snow2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1168	tab_flds(:,:,ids_toteva) = cpl_evap2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1169	cpl_evap2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1170	tab_flds(:,:,ids_taumod) = cpl_taumod2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1171	cpl_taumod2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1172
1173	ENDWHERE
1174
1175	tab_flds(:,:,ids_icevap) = cpl_evap2D(:,:,2)
1176
1177	ELSE IF (version_ocean=='opa8') THEN
1178	! Store fields for rain and snow in tab_flds(:,:,15) and tab_flds(:,:,16)
1179	tab_flds(:,:,ids_totrai) = 0.0
1180	tab_flds(:,:,ids_totsno) = 0.0
1181	tmp_taux(:,:) = 0.0
1182	tmp_tauy(:,:) = 0.0
1183	! For all valid grid cells containing some fraction of ocean or sea-ice
1184	WHERE ( deno(:,:) /= 0 )
1185	tab_flds(:,:,ids_totrai) = cpl_rain2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1186	cpl_rain2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1187	tab_flds(:,:,ids_totsno) = cpl_snow2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1188	cpl_snow2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1189
1190	tmp_taux = cpl_taux2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1191	cpl_taux2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1192	tmp_tauy = cpl_tauy2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) + &
1193	cpl_tauy2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
1194	ENDWHERE
1195	END IF
1196
1197	ENDIF ! is_omp_root
1198
1199	!*************************************************************************************
1200	! Transform the wind components from local atmospheric 2D coordinates to geocentric
1201	! 3D coordinates.
1202	! Store the resulting wind components in tab_flds(:,:,1:6)
1203	!*************************************************************************************
1204
1205	! Transform the longitudes and latitudes on 2D arrays
1206
1207	CALL gather_omp(rlon,rlon_mpi)
1208	CALL gather_omp(rlat,rlat_mpi)
1209	!$OMP MASTER
1210	CALL Grid1DTo2D_mpi(rlon_mpi,tmp_lon)
1211	CALL Grid1DTo2D_mpi(rlat_mpi,tmp_lat)
1212	!$OMP END MASTER
1213
1214	IF (is_sequential) THEN
1215	IF (is_north_pole) tmp_lon(:,1) = tmp_lon(:,2)
1216	IF (is_south_pole) tmp_lon(:,jjm+1) = tmp_lon(:,jjm)
1217	ENDIF
1218
1219	! NetCDF output of the wind before transformation of coordinate system
1220	IF (is_sequential) THEN
1221	ndexct(:) = 0
1222	itau_w = itau_phy + itime
1223	CALL histwrite(nidct,'tauxe',itau_w,tmp_taux,iim*(jjm+1),ndexct)
1224	CALL histwrite(nidct,'tauyn',itau_w,tmp_tauy,iim*(jjm+1),ndexct)
1225	CALL histwrite(nidct,'tmp_lon',itau_w,tmp_lon,iim*(jjm+1),ndexct)
1226	CALL histwrite(nidct,'tmp_lat',itau_w,tmp_lat,iim*(jjm+1),ndexct)
1227	ENDIF
1228
1229	! Transform the wind from spherical atmospheric 2D coordinates to geocentric
1230	! cartesian 3D coordinates
1231	!$OMP MASTER
1232	CALL atm2geo (iim, jj_nb, tmp_taux, tmp_tauy, tmp_lon, tmp_lat, &
1233	tab_flds(:,:,ids_tauxxu), tab_flds(:,:,ids_tauyyu), tab_flds(:,:,ids_tauzzu) )
1234
1235	tab_flds(:,:,ids_tauxxv) = tab_flds(:,:,ids_tauxxu)
1236	tab_flds(:,:,ids_tauyyv) = tab_flds(:,:,ids_tauyyu)
1237	tab_flds(:,:,ids_tauzzv) = tab_flds(:,:,ids_tauzzu)
1238	!$OMP END MASTER
1239
1240	!*************************************************************************************
1241	! NetCDF output of all fields just before sending to coupler.
1242	!
1243	!*************************************************************************************
1244	IF (is_sequential) THEN
1245	DO j=1,maxsend
1246	IF (infosend(j)%action) CALL histwrite(nidct,infosend(j)%name, itau_w, &
1247	tab_flds(:,:,j),iim*(jjm+1),ndexct)
1248	ENDDO
1249	ENDIF
1250	!*************************************************************************************
1251	! Send the table of all fields
1252	!
1253	!*************************************************************************************
1254	time_sec=(itime-1)*dtime
1255	#ifdef CPP_COUPLE
1256	!$OMP MASTER
1257	CALL intocpl(time_sec, lafin, tab_flds(:,:,:))
1258	!$OMP END MASTER
1259	#endif
1260
1261	!*************************************************************************************
1262	! Finish with some dellocate
1263	!
1264	!*************************************************************************************
1265	sum_error=0
1266	DEALLOCATE(cpl_sols2D, cpl_nsol2D, cpl_rain2D, cpl_snow2D, stat=error )
1267	sum_error = sum_error + error
1268	DEALLOCATE(cpl_evap2D, cpl_tsol2D, cpl_fder2D, cpl_albe2D, stat=error )
1269	sum_error = sum_error + error
1270	DEALLOCATE(cpl_taux2D, cpl_tauy2D, cpl_windsp2D, cpl_taumod2D, stat=error )
1271	sum_error = sum_error + error
1272
1273	IF (carbon_cycle_cpl) THEN
1274	DEALLOCATE(cpl_atm_co22D, stat=error )
1275	sum_error = sum_error + error
1276	END IF
1277
1278	IF (sum_error /= 0) THEN
1279	abort_message='Pb in deallocation of cpl_xxxx2D coupling variables'
1280	CALL abort_gcm(modname,abort_message,1)
1281	ENDIF
1282
1283	END SUBROUTINE cpl_send_all
1284	!
1285	!*************************************************************************************
1286	!
1287	SUBROUTINE cpl2gath(champ_in, champ_out, knon, knindex)
1288	USE mod_phys_lmdz_para
1289	! Cette routine transforme un champs de la grille 2D recu du coupleur sur la grille
1290	! 'gathered' (la grille physiq comprime).
1291	!
1292	!
1293	! input:
1294	! champ_in champ sur la grille 2D
1295	! knon nombre de points dans le domaine a traiter
1296	! knindex index des points de la surface a traiter
1297	!
1298	! output:
1299	! champ_out champ sur la grille 'gatherd'
1300	!
1301	INCLUDE "dimensions.h"
1302
1303	! Input
1304	INTEGER, INTENT(IN) :: knon
1305	REAL, DIMENSION(iim,jj_nb), INTENT(IN) :: champ_in
1306	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
1307
1308	! Output
1309	REAL, DIMENSION(klon_mpi), INTENT(OUT) :: champ_out
1310
1311	! Local
1312	INTEGER :: i, ig
1313	REAL, DIMENSION(klon_mpi) :: temp_mpi
1314	REAL, DIMENSION(klon) :: temp_omp
1315
1316	!*************************************************************************************
1317	!
1318
1319
1320	! Transform from 2 dimensions (iim,jj_nb) to 1 dimension (klon)
1321	!$OMP MASTER
1322	CALL Grid2Dto1D_mpi(champ_in,temp_mpi)
1323	!$OMP END MASTER
1324
1325	CALL scatter_omp(temp_mpi,temp_omp)
1326
1327	! Compress from klon to knon
1328	DO i = 1, knon
1329	ig = knindex(i)
1330	champ_out(i) = temp_omp(ig)
1331	ENDDO
1332
1333	END SUBROUTINE cpl2gath
1334	!
1335	!*************************************************************************************
1336	!
1337	SUBROUTINE gath2cpl(champ_in, champ_out, knon, knindex)
1338	USE mod_phys_lmdz_para
1339	! Cette routine ecrit un champ 'gathered' sur la grille 2D pour le passer
1340	! au coupleur.
1341	!
1342	! input:
1343	! champ_in champ sur la grille gathere
1344	! knon nombre de points dans le domaine a traiter
1345	! knindex index des points de la surface a traiter
1346	!
1347	! output:
1348	! champ_out champ sur la grille 2D
1349	!
1350	INCLUDE "dimensions.h"
1351
1352	! Input arguments
1353	!*************************************************************************************
1354	INTEGER, INTENT(IN) :: knon
1355	REAL, DIMENSION(klon), INTENT(IN) :: champ_in
1356	INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
1357
1358	! Output arguments
1359	!*************************************************************************************
1360	REAL, DIMENSION(iim,jj_nb), INTENT(OUT) :: champ_out
1361
1362	! Local variables
1363	!*************************************************************************************
1364	INTEGER :: i, ig
1365	REAL, DIMENSION(klon) :: temp_omp
1366	REAL, DIMENSION(klon_mpi) :: temp_mpi
1367	!*************************************************************************************
1368
1369	! Decompress from knon to klon
1370	temp_omp = 0.
1371	DO i = 1, knon
1372	ig = knindex(i)
1373	temp_omp(ig) = champ_in(i)
1374	ENDDO
1375
1376	! Transform from 1 dimension (klon) to 2 dimensions (iim,jj_nb)
1377	CALL gather_omp(temp_omp,temp_mpi)
1378
1379	!$OMP MASTER
1380	CALL Grid1Dto2D_mpi(temp_mpi,champ_out)
1381
1382	IF (is_north_pole) champ_out(:,1)=temp_mpi(1)
1383	IF (is_south_pole) champ_out(:,jj_nb)=temp_mpi(klon)
1384	!$OMP END MASTER
1385
1386	END SUBROUTINE gath2cpl
1387	!
1388	!*************************************************************************************
1389	!
1390	END MODULE cpl_mod
1391

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

Download in other formats:

Original Format