SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_types_glt.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 !GLT_LIC The GELATO model is a seaice model used in stand-alone or embedded mode.
6 !GLT_LIC It has been developed by Meteo-France. The holder of GELATO is Meteo-France.
7 !GLT_LIC
8 !GLT_LIC This software is governed by the CeCILL-C license under French law and biding
9 !GLT_LIC by the rules of distribution of free software. See the CeCILL-C_V1-en.txt
10 !GLT_LIC (English) and CeCILL-C_V1-fr.txt (French) for details. The CeCILL is a free
11 !GLT_LIC software license, explicitly compatible with the GNU GPL
12 !GLT_LIC (see http://www.gnu.org/licenses/license-list.en.html#CeCILL)
13 !GLT_LIC
14 !GLT_LIC The CeCILL-C licence agreement grants users the right to modify and re-use the
15 !GLT_LIC software governed by this free software license. The exercising of this right
16 !GLT_LIC is conditional upon the obligation to make available to the community the
17 !GLT_LIC modifications made to the source code of the software so as to contribute to
18 !GLT_LIC its evolution.
19 !GLT_LIC
20 !GLT_LIC In consideration of access to the source code and the rights to copy, modify
21 !GLT_LIC and redistribute granted by the license, users are provided only with a limited
22 !GLT_LIC warranty and the software's author, the holder of the economic rights, and the
23 !GLT_LIC successive licensors only have limited liability. In this respect, the risks
24 !GLT_LIC associated with loading, using, modifying and/or developing or reproducing the
25 !GLT_LIC software by the user are brought to the user's attention, given its Free
26 !GLT_LIC Software status, which may make it complicated to use, with the result that its
27 !GLT_LIC use is reserved for developers and experienced professionals having in-depth
28 !GLT_LIC computer knowledge. Users are therefore encouraged to load and test the
29 !GLT_LIC suitability of the software as regards their requirements in conditions enabling
30 !GLT_LIC the security of their systems and/or data to be ensured and, more generally, to
31 !GLT_LIC use and operate it in the same conditions of security.
32 !GLT_LIC
33 !GLT_LIC The GELATO sofware is cureently distibuted with the SURFEX software, available at
34 !GLT_LIC http://www.cnrm.meteo.fr/surfex. The fact that you download the software deemed that
35 !GLT_LIC you had knowledge of the CeCILL-C license and that you accept its terms.
36 !GLT_LIC Attempts to use this software in a way not complying with CeCILL-C license
37 !GLT_LIC may lead to prosecution.
38 !GLT_LIC
39 ! =======================================================================
40 ! ========================= MODULE modd_types_glt ===========================
41 ! =======================================================================
42 !
43 ! Goal:
44 ! -----
45 ! This module contains the definitions of all the structured
46 ! variables used in Gelato.
47 !
48 ! Created : 1996 (D. Salas y Melia)
49 ! Modified: 2009/08 (D. Salas y Melia)
50 ! Remove obsolete fields, add new types_glt used by Gelato's new
51 ! interface (allowing double or single physics), re-classify
52 ! Modified: 2010/06 (D. Salas y Melia)
53 ! New CMIP5 diagnostics
54 ! Modified: 2011/12 (A. Voldoire)
55 ! New diagnostics were added to compute water balance
56 ! Modified: 2012/01 (M. Chevallier)
57 ! Explicit melt pond parameterisation
58 ! Modified: 2012/11 (D. Salas y Melia)
59 ! Introduce a 'super-structure'. Add ice-ocean stress
60 ! components to t_tfl
61 ! Modified: 2013/07 (S. Senesi) : #ifdef in_surfex for dynamic arrays
62 ! in t_glt
63 ! -----------------------------------------------------------------------
64 !
66 !
68 !
69 !
70 ! 1. Structure defining the domain of study
71 ! ==========================================
72 !
73 TYPE t_dom
74 ! If tmk = 1, the grid T-point is ocean, else it is land.
75  INTEGER :: &
76  tmk
77 ! If umk = 1, the grid U-point is ocean, else it is land.
78  INTEGER :: &
79  umk
80 ! If vmk = 1, the grid V-point is ocean, else it is land.
81  INTEGER :: &
82  vmk
83 ! If imk = 1, the grid T-point is ocean, else it is land. Ghost points masked.
84  INTEGER :: &
85  imk
86 ! Values of i and j index in the global grid
87  INTEGER :: &
88  indi,indj
89 ! Grid point longitude (radians)
90  REAL :: &
91  lon
92 ! Grid point latitude (radians)
93  REAL :: &
94  lat
95 ! Grid cell length (m) in the X-direction
96  REAL :: &
97  dxc
98 ! Grid cell length (m) in the Y-direction
99  REAL :: &
100  dyc
101 ! Grid cell surface (=dxc*dyc), ghost points masked
102  REAL :: &
103  srf
104 END TYPE t_dom
105 !
106 !
107 !
108 ! 2. Structures defining the forcing data
109 ! ========================================
110 !
111 ! 2.1. Ocean mixed layer forcing
112 ! -------------------------------
113 !
114 ! .. This type contains all variables describing the state of the
115 ! ocean mixed layer, that influences sea ice and leads from the bottom.
116 !
117 TYPE t_mxl
118 !
119 ! Frazil flux from ocean, converted in equivalent W.m-2
120  REAL :: &
121  qml
122 ! Ocean heat flux, per m2 of marine surface (W.m-2)
123  REAL :: &
124  qoc
125 ! Temperature (K)
126  REAL :: &
127  tml
128 ! Freezing temperature (K)
129  REAL :: &
130  mlf
131 ! Salinity (psu)
132  REAL :: &
133  sml
134 ! SSH (m)
135  REAL :: &
136  ssh
137 ! Velocity (m.s-1)
138  REAL :: &
139  uml,vml
140 ! Mixed layer heat content wrt ice melting T / kg, i.e. Cp x (Tml-Tice_m)
141  REAL &
142  hco
143 END TYPE t_mxl
144 !
145 !
146 ! 2.2. Atmospheric forcing, independent of the nature of the surface
147 ! -------------------------------------------------------------------
148 !
149 ! .. Liquid and solid precipitation does not depend on the surface on
150 ! which they fall. We suppose here that the stress is not different on
151 ! a snow, water or bare ice surface, which is not totally the case in
152 ! reality.
153 !
154 TYPE t_atm
155 !
156 ! Liquid precipitations (m/time step)
157  REAL :: &
158  lip
159 ! Solid precipitations (m/time step)
160  REAL :: &
161  sop
162 ! Zonal wind stress (N.m-2)
163  REAL :: &
164  ztx
165 ! Meridional wind stress (N.m-2)
166  REAL :: &
167  mty
168 END TYPE t_atm
169 !
170 !
171 ! 2.3. Atmospheric forcing, dependent on the nature of the surface
172 ! -----------------------------------------------------------------
173 !
174 ! .. This structure may concern sea ice or open water. It contains:
175 ! - short wave absorbed flux
176 ! - non solar heat flux (i.e the sum of long wave, latent and
177 ! sensible heat fluxes)
178 ! - the derivative of this flux by surface temperature
179 ! - the evaporation rate
180 !
181 TYPE t_blk
182 ! Short wave flux absorbed by the surface (W.m-2)
183  REAL :: &
184  swa
185 ! Non solar flux affecting the surface (W.m-2)
186  REAL :: &
187  nsf
188 ! Derivative of non solar flux by the surface temperature (W.m-2.K-1)
189  REAL :: &
190  dfl
191 ! Evaporation rate (m.s-1)
192  REAL :: &
193  eva
194 END TYPE t_blk
195 !
196 !
197 !
198 ! 3. Structures describing sea ice dynamic and thermodynamic states
199 ! ==================================================================
200 !
201 ! 3.1. Sea ice dynamics (EVP, JFNK)
202 ! ----------------------------------
203 !
204 TYPE t_evp
205 ! (u,v) velocity component
206  REAL :: &
207  uvl,vvl
208 ! Stress tensor components
209  REAL :: &
210  oone,oose,oosw,oonw
211  REAL :: &
212  otne,otse,otsw,otnw
213  REAL :: &
214  ttne,ttse,ttsw,ttnw
215 END TYPE t_evp
216 !
217 TYPE t_jfn
218 ! (u,v) velocity components
219  REAL :: &
220  uvl,vvl
221 END TYPE t_jfn
222 !
223 !
224 ! 3.2. Sea ice thermodynamics - 0D variables
225 ! -------------------------------------------
226 !
227 TYPE t_sit
228 !
229 ! Boolean that is TRUE if ice class is empty
230  LOGICAL :: &
231  esi
232 ! Ice age
233  REAL :: &
234  age
235 ! Snow layer albedo
236  REAL :: &
237  asn
238 ! Areal fraction for each ice type
239  REAL :: &
240  fsi
241 ! Ice thickness for each ice type
242  REAL :: &
243  hsi
244 ! Ice salinity for each ice type
245  REAL :: &
246  ssi
247 ! Snow thickness over each ice type
248  REAL :: &
249  hsn
250 ! Density of snow over each ice type
251  REAL :: &
252  rsn
253 ! Ice or snow surface temperature
254  REAL :: &
255  tsf
256 ! Melt pond volume
257  REAL :: &
258  vmp
259 END TYPE t_sit
260 !
261 !
262 ! 3.3. Sea ice thermodynamics - 1D variables
263 ! -------------------------------------------
264 !
265 TYPE t_vtp
266 ! Vertical gltools_enthalpy profile in the slab (J.kg-1)
267  REAL :: &
268  ent
269 ! Vertical salinity profile (g.kg-1)
270  REAL :: &
271  vsp
272 END TYPE t_vtp
273 !
274 !
275 ! 3.4. Leads - 0D variables
276 ! --------------------------
277 !
278 ! .. May be removed in subsequent versions, unless more physics takes
279 ! places in leads
280 !
281 TYPE t_opw
282 ! Areal fraction covered by the lead
283  REAL :: &
284  fld
285 ! Lead temperature (K)
286  REAL :: &
287  tld
288 !
289 END TYPE t_opw
290 !
291 !
292 !
293 ! 4. Output structures
294 ! =====================
295 !
296 ! 4.1. Fluxes transmitted to the mixed layer
297 ! -------------------------------------------
298 !
299 ! NOTE: This is internal to Gelato (not involved in Gelato's interface
300 ! with the ocean model)
301 !
302 ! .. This structure contains all heat, non solar and solar short wave heat
303 ! fluxes that cross leads or sea ice to reach the mixed layer. The
304 ! interest of these fluxes lies mainly in the framework of ocean/sea ice
305 ! coupling, and the classifications of variables is the following
306 ! inside the type, if type_var%xyz
307 ! x = (l,t,s,c,w) are respectively short wave, non solar,
308 ! concentration/dilution equivalent water flux (also called virtual
309 ! water flux) and fresh water mass fluxes. Heat fluxes are in W.m-2 and
310 ! water virtual fluxes, water fluxes and salt fluxes are in kg.m-2.s-1.
311 ! y or z = (i,l,o) are respectively sea ice, lead and mixed layer.
312 ! y is the medium the flux originates from, and z is the destination
313 ! medium.
314 !
315 ! * Example : flx%wio is the fresh water flux from sea ice to the mixed
316 ! layer.
317 !
318 ! .. Note that this structure also contains the ice-ocean stress
319 ! components along x and y (respectively xio and yio)
320 !
321 TYPE t_tfl
322  REAL :: &
323  lio,llo
324  REAL :: &
325  tio,tlo
326  REAL :: &
327  sio ! Note that slo is always zero
328  REAL :: &
329  cio ! Note that clo would be identical to wlo
330  REAL :: &
331  wio,wlo
332  REAL :: &
333  xio,yio
334 END TYPE t_tfl
335 !
336 !
337 ! 4.2. Fluxes transmitted to the ocean model
338 ! -------------------------------------------
339 !
340 ! NOTE: In contrast with t_tfl, this type is transmitted to the ocean
341 ! model (not strictly internal to Gelato)
342 !
343 TYPE t_2oc
344 !
345 ! Net non solar heat flux (W.m-2)
346  REAL :: &
347  nsf
348 ! Net short wave heat flux (W.m-2)
349  REAL :: &
350  swa
351 ! Concentration/dilution flux (kg.m-2.s-1)
352  REAL :: &
353  cdf
354 ! Salt flux (kg.m-2.s-1)
355  REAL :: &
356  saf
357 ! Water flux (kg.m-2.s-1)
358  REAL :: &
359  wfl
360 ! Stress X-component (N.m-2)
361  REAL :: &
362  ztx
363 ! Stress Y-component (N.m-2)
364  REAL :: &
365  mty
366 ! Ice-Ocean friction velocity U* (m.s-2)
367  REAL :: &
368  ust
369 END TYPE t_2oc
370 !
371 !
372 ! 4.3. Surface description, sent to the atmosphere (through the ocean model ?)
373 ! -----------------------------------------------------------------------------
374 !
375 ! .. This structure includes a fraction, albedo and temperature. Two cases
376 ! arise:
377 ! - if single physics is used (only one input flux to be shared between
378 ! the different kinds of surface): the fraction should be the total
379 ! fraction of sea ice. The temperature and albedo should be weighted over
380 ! water and sea ice.
381 ! - if double physics is used, the fraction is stille the total fraction
382 ! of sea ice, but the temperature and albedo are now weighted over
383 ! sea ice only.
384 !
385 TYPE t_2at
386 !
387 ! Sea ice fraction (0-1)
388  REAL :: &
389  fsi
390 ! Albedo (0-1)
391  REAL :: &
392  alb
393 ! Temperature (K)
394  REAL :: &
395  tsf
396 END TYPE t_2at
397 !
398 !
399 ! 4.4. Energy and water budget
400 ! -----------------------------
401 !
402 TYPE t_bud
403 !
404 ! Sea ice/snow gltools_enthalpy in the grid cell, beg. of time step (J)
405  REAL :: &
406  eni
407 ! Sea ice/snow gltools_enthalpy in the grid cell, current (J)
408  REAL :: &
409  enn
410 ! Input heat flux on sea ice in the grid cell (bottom)
411  REAL :: &
412  bii
413 ! Input heat flux due to snow falls on sea ice
414  REAL :: &
415  nii
416 ! Input heat flux due to snow falls on leads
417  REAL :: &
418  nli
419 ! Input heat flux on sea ice in the grid cell (top)
420  REAL :: &
421  hii
422 ! Input heat flux on leads in the grid cell
423  REAL :: &
424  hli
425 ! Output heat flux on sea ice in the grid cell
426  REAL :: &
427  hio
428 ! Output heat flux on leads in the grid cell
429  REAL :: &
430  hlo
431 ! Input water flux on sea ice in the grid cell
432  REAL :: &
433  wii
434 ! Input water flux on leads in the grid cell
435  REAL :: &
436  wli
437 ! Grid cell water stored in sea ice & snow, beg. of time step (kg.m-2)
438  REAL :: &
439  fwi
440 ! Grid cell water stored in sea ice & snow, current (kg.m-2)
441  REAL :: &
442  fwn
443 ! Grid cell salt stored in sea ice, beg. of time step (kg.m-2)
444  REAL :: &
445  isi
446 ! Grid cell salt stored in sea ice, current (kg.m-2)
447  REAL :: &
448  isn
449 !
450 END TYPE t_bud
451 !
452 !
453 ! 4.5. Diagnostics
454 ! -----------------
455 !
456 ! .. Note that here the sea ice net bottom melting [kg.m-2.s-1] is deduced
457 ! from dsi-lsi-mrt-mrl
458 !
459 ! FIELDS
460 !
461 TYPE t_dia
462 !
463 ! Sea ice u and v velocity on C-grid (m.s-1)
464  REAL :: &
465  uvl,vvl
466 ! Weights for bare sea ice albedo
467  REAL :: &
468  aiw
469 ! Bare sea ice albedo, averaged over all categories * total ice fraction [0-1]
470  REAL :: &
471  asi
472 ! Melt pond albedo, averaged over all categories * total ice fraction [0-1]
473  REAL :: &
474  amp
475 ! Snow albedo, averaged over all categories * total ice fraction [0-1]
476  REAL :: &
477  asn
478 ! Atmospheric stress on sea ice (X-component) * total ice fraction [N.m-2]
479  REAL :: &
480  atx
481 ! Atmospheric stress on sea ice (Y-component) * total ice fraction [N.m-2]
482  REAL :: &
483  aty
484 ! Congelation sea ice growth rate [kg.m-2.s-1]
485  REAL :: &
486  cgl
487 ! Salt mass net change for all glt_thermo processes [kg.m-2.s-1]
488  REAL :: &
489  dsa
490 ! Salt mass net change due to dynamic processes [kg.m-2.s-1]
491  REAL :: &
492  dds
493 ! Snow mass net change for all glt_thermo processes [kg.m-2.s-1]
494  REAL :: &
495  dsn
496 ! Snow mass net change due to dynamic processes [kg.m-2.s-1]
497  REAL :: &
498  ddn
499 ! Sea ice mass change for all glt_thermo processes [kg.m-2.s-1]
500  REAL :: &
501  dsi
502 ! Sea ice mass change due to dynamic processes [kg.m-2.s-1]
503  REAL :: &
504  ddi
505 ! Sea ice mass change due to the constraint [kg.m-2.s-1]
506  REAL :: &
507  dci
508 ! Sea ice concentration for constraint [%]
509  REAL :: &
510  cst
511 ! Sea ice fresh water content change for all glt_thermo processes [kg.m-2.s-1]
512  REAL :: &
513  dwi
514 ! Liquid precipitation falling on the sea ice portion of grid cell [kg.m-2.s-1]
515  REAL :: &
516  lip
517 ! Sea ice production in leads [kg.m-2.s-1]
518  REAL :: &
519  lsi
520 ! Sea ice net bottom melting [kg.m-2.s-1]
521  REAL :: &
522  mrb
523 ! Sea ice net top melting [kg.m-2.s-1]
524  REAL :: &
525  mrt
526 ! Sea ice lateral melting [kg.m-2.s-1]
527  REAL :: &
528  mrl
529 ! Ocean stress on sea ice (X-component) * total ice fraction [N.m-2]
530  REAL :: &
531  otx
532 ! Ocean stress on sea ice (Y-component) * total ice fraction [N.m-2]
533  REAL :: &
534  oty
535 ! Enthalpy of the sea-ice part of the ice-snow slab [J.m-2]
536  REAL :: &
537  sie
538 ! Enthalpy of the snow part of the ice-snow slab [J.m-2]
539  REAL :: &
540  sne
541 ! Snow-Ice formation rate [kg.m-2.s-1]
542  REAL :: &
543  sni
544 ! Snow melt rate [kg.m-2.s-1]
545  REAL :: &
546  snm
547 ! Snow melt rate due to lateral ice erosion [kg.m-2.s-1]
548  REAL :: &
549  snml
550 ! Solid precipitation falling on the sea ice portion of grid cell [kg.m-2.s-1]
551  REAL :: &
552  sop
553 ! Net shortwave on the sea ice portion of grid cell [W.m-2]
554  REAL :: &
555  swi
556 ! Net shortwave on the leads portion of the grid cell [W.m-2]
557  REAL :: &
558  sww
559 ! Sea ice gltools_enthalpy change due to the thermodynamics [W.m-2]
560  REAL :: &
561  the
562 ! Temperature at snow-ice interface * ice fraction [K]
563  REAL :: &
564  tin
565 ! Ice concentration weights for temperature at snow-ice interface
566  REAL :: &
567  tiw
568 ! Input Fresh Water into the snow-ice-leads system [kg.m-2.s-1]
569  REAL :: &
570  ifw
571 ! Accumulated total sea ice concentation over all time steps [0-1]
572  REAL :: &
573  sic
574 ! Accumulated average sea ice thickness over all time steps [m]
575  REAL :: &
576  sit
577 ! Accumulated average snow thickness over all time steps [m]
578  REAL :: &
579  snd
580 ! Ocean-ice heat flux per m2 of marine surface [W.m-2]
581  REAL :: &
582  qoi
583 ! Transport X and Y-components [kg.s-1]
584  REAL :: &
585  xtr,ytr
586 ! Spare fields
587  REAL :: &
588  sp1,sp2
589 ! Sublimation
590  REAL :: &
591  sut
592 ! Sublimation over ice
593  REAL :: &
594  sui
595 ! Sublimation over snow
596  REAL :: &
597  sus
598 ! Sublimation over ocean
599  REAL :: &
600  suw
601 ! Sublimation over leads
602  REAL :: &
603  sul
604 ! Liquid precip that actually falls on snow
605  REAL :: &
606  s_pr
607 ! Liquid precip that actually falls on ocean
608  REAL :: &
609  o_pr
610 ! Liquid precip that actually falls on leads
611  REAL :: &
612  l_pr
613 ! Solid precip that actually falls on snow
614  REAL :: &
615  s_prsn
616 ! Solid precip that actually falls on ocean
617  REAL :: &
618  o_prsn
619 ! Solid precip that actually falls on leads
620  REAL :: &
621  l_prsn
622 ! terms of virtual salt flux (sublimation, snow2ice, updhsi, lmltsi, updsal)
623  REAL :: &
624  subcio, snicio, hsicio, lmlcio, salcio
625 ! Total damping/restoring (in equivalent energy flux)
626  REAL :: &
627  dmp
628 END TYPE t_dia
629 !
630 ! FIELD INDICES, COUNTERS
631 !
632 TYPE t_ind
633 !
634 ! Current time-step
635  INTEGER :: &
636  cur
637 ! First time-step
638  INTEGER :: &
639  beg
640 ! Last time-step
641  INTEGER :: &
642  end
643 ! Current 0d, 2d field index
644  INTEGER :: &
645  i0d,i2d
646 ! Number of time steps in accumulating
647  INTEGER :: &
648  nts
649 END TYPE t_ind
650 !
651 !
652 ! 4.6. Definitions of diagnostics
653 ! --------------------------------
654 !
655 TYPE t_def
656 !
657 ! Field definition
658  CHARACTER(80) :: &
659  def
660 ! Field long name
661  CHARACTER(80) :: &
662  lna
663 ! Field short name
664  CHARACTER(80) :: &
665  sna
666 ! Field units
667  CHARACTER(80) :: &
668  uni
669 ! Field location ('T', 'U' or 'V')
670  CHARACTER(80) :: &
671  loc
672 ! Field nature ('SCALAR' or 'VECTOR')
673  CHARACTER(80) :: &
674  typ
675 END TYPE t_def
676 !
677 !
678 !
679 ! 5. Super-structure
680 ! ===================
681 !
682 ! This structures includes all structures that are needed by the different
683 ! parts of Gelato that are invoked from parts of the host program (e.g.
684 ! surface-scheme or ocean model)
685 !
686 TYPE t_glt
687 !
688 !
689 ! 5.1. Time-step information
690 ! ---------------------------
691 !
692  TYPE(t_ind) :: &
693  ind
694 !
695 !
696 ! 5.2. Domain definition
697 ! -----------------------
698 !
699 ! Bathymetry (in meters)
700  REAL, DIMENSION(:,:), POINTER :: &
701  bat
702 ! Main domain
703  TYPE(t_dom), DIMENSION(:,:), POINTER :: &
704  dom
705 #if ! defined in_surfex
706 ! Domain for dynamics (EVP) and advection
707  REAL :: &
708  xymin
709  REAL, DIMENSION(:,:), POINTER :: &
710  dxtr,dytr,fcor,tarear
711  REAL, DIMENSION(:,:), POINTER :: &
712  HTN,HTE,HTS,HTW,dxtr4,dytr4,hm,um, &
713  dxta,dyta,dxt2r,dyt2r,HTNa,HTEa, &
714  cx,cy,cxx,cxy,cyy,cxxx,cxxy,cxyy,cyyy
715  REAL, DIMENSION(:,:,:,:), POINTER :: &
716  mne,mnw,msw,mse
717 ! Index mapping, for 1d packed vectors
718  INTEGER, DIMENSION(:,:,:), POINTER :: &
719  index2d
720 ! Ice mask, to pack arrays to 1d vectors
721  REAL, DIMENSION(:,:,:), POINTER :: &
722  mskice
723 #endif
724 !
725 !
726 ! 5.3. Input/Output of the glt_gelato main routine
727 ! ---------------------------------------------
728 !
729 ! Ocean forcing
730  TYPE(t_mxl), DIMENSION(:,:), POINTER :: &
731  oce_all
732 ! Input atmospheric forcing (independent on the nature of the surface)
733  TYPE(t_atm), DIMENSION(:,:), POINTER :: &
734  atm_all
735 ! Atmospheric forcing (potentially dependent on the nature of the surface)
736  TYPE(t_blk), DIMENSION(:,:,:), POINTER :: &
737  atm_ice, atm_mix
738 ! Atmospheric forcing on water
739  TYPE(t_blk), DIMENSION(:,:), POINTER :: &
740  atm_wat
741 ! Output to the ocean
742  TYPE(t_2oc), DIMENSION(:,:), POINTER :: &
743  all_oce
744 ! Output to the atmosphere
745  TYPE(t_2at), DIMENSION(:,:,:), POINTER :: &
746  ice_atm, mix_atm
747 ! Constraint on sea ice
748  TYPE(t_sit), DIMENSION(:,:,:), POINTER :: &
749  sit_d
750 !
751 !
752 ! 5.4. Ice and mixed layer states
753 ! --------------------------------
754 !
755  TYPE(t_evp), DIMENSION(:,:), POINTER :: &
756  evp
757 !
758  TYPE(t_jfn), DIMENSION(:,:), POINTER :: &
759  jfn
760 !
761  TYPE(t_sit), DIMENSION(:,:,:), POINTER :: &
762  sit
763 !
764  TYPE(t_vtp), DIMENSION(:,:,:,:), POINTER :: &
765  sil
766 !
767  TYPE(t_mxl), DIMENSION(:,:), POINTER :: &
768  tml
769 !
770 !
771 ! 5.5. Fluxes, budgets, diagnostics
772 ! ----------------------------------
773 !
774  REAL, DIMENSION(:,:), POINTER :: &
775  ust
776 !
777  REAL, DIMENSION(:,:,:), POINTER :: &
778  cdia0,cdia
779 !
780  TYPE(t_blk), DIMENSION(:,:), POINTER :: &
781  blkw
782 !
783  TYPE(t_blk), DIMENSION(:,:,:), POINTER :: &
784  blki
785 !
786  TYPE(t_tfl), DIMENSION(:,:), POINTER :: &
787  tfl
788 !
789  TYPE(t_bud), DIMENSION(:,:), POINTER :: &
790  bud
791 !
792  TYPE(t_dia), DIMENSION(:,:), POINTER :: &
793  dia
794 !
795 END TYPE t_glt
796 !
797 END MODULE modd_types_glt