New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
limdia_2.F90 in branches/dev_002_LIM/NEMO/LIM_SRC_2 – NEMO

source: branches/dev_002_LIM/NEMO/LIM_SRC_2/limdia_2.F90 @ 823

Last change on this file since 823 was 823, checked in by rblod, 16 years ago

Final step to rename LIM_SRC in LIM_SRC_2

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 12.8 KB
Line 
1MODULE limdia_2
2   !!======================================================================
3   !!                       ***  MODULE limdia_2   ***
4   !!                      diagnostics of ice model
5   !!======================================================================
6   !! History :   8.0  !  97-06  (Louvain-La-Neuve)  Original code
7   !!             8.5  !  02-09  (C. Ethe , G. Madec )  F90: Free form and module
8   !!             9.0  !  06-08  (S. Masson)  change frequency output control
9   !!-------------------------------------------------------------------
10#if defined key_lim2
11   !!----------------------------------------------------------------------
12   !!   'key_lim2' :                                  LIM 2.0 sea-ice model
13   !!----------------------------------------------------------------------
14   !!----------------------------------------------------------------------
15   !!   lim_dia_2      : computation of the time evolution of keys var.
16   !!   lim_dia_init_2 : initialization and namelist read
17   !!----------------------------------------------------------------------
18   USE phycst          !
19   USE par_ice_2       ! ice parameters
20   USE ice_oce         ! ice variables
21   USE daymod          !
22   USE dom_ice_2       !
23   USE ice_2           !
24   USE limistate_2     !
25   USE in_out_manager  ! I/O manager
26
27   IMPLICIT NONE
28   PRIVATE
29
30   PUBLIC               lim_dia_2          ! called by ice_step
31   INTEGER, PUBLIC ::   ntmoy   = 1 ,   &  !: instantaneous values of ice evolution or averaging ntmoy
32      &                 ninfo   = 1        !: frequency of ouputs on file ice_evolu in case of averaging
33
34   INTEGER, PARAMETER ::   &  ! Parameters for outputs to files "evolu"
35      jpinfmx = 100         ,    &  ! maximum number of key variables
36      jpchinf = 5           ,    &  ! ???
37      jpchsep = jpchinf + 2         ! ???
38
39   INTEGER ::   &
40      nfrinf  = 4 ,     &  ! number of variables written in one line
41      nferme ,          &  ! last time step at which the var. are written on file
42      nvinfo ,          &  ! number of total variables
43      nbvt   ,          &  ! number of time variables
44      naveg                ! number of step for accumulation before averaging
45
46   CHARACTER(len= 8) ::   fmtinf  = '1PE13.5 ' ! format of the output values 
47   CHARACTER(len=30) ::   fmtw  ,           &  ! formats
48      &                   fmtr  ,           &  ! ???
49      &                   fmtitr               ! ???
50   CHARACTER(len=jpchsep), DIMENSION(jpinfmx) ::   titvar               ! title of key variables
51 
52   REAL(wp)                     ::   epsi06 = 1.e-06      ! ???
53   REAL(wp), DIMENSION(jpinfmx) ::   vinfom               ! temporary working space
54   REAL(wp), DIMENSION(jpi,jpj) ::   aire                 ! masked grid cell area
55
56   !! * Substitutions
57#  include "vectopt_loop_substitute.h90"
58   !!----------------------------------------------------------------------
59   !!   LIM 2.0,  UCL-LOCEAN-IPSL (2005)
60   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/LIM_SRC/limdia.F90,v 1.9 2007/06/29 17:03:12 opalod Exp $
61   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
62   !!----------------------------------------------------------------------
63
64CONTAINS
65
66   SUBROUTINE lim_dia_2( kt )
67      !!--------------------------------------------------------------------
68      !!                  ***  ROUTINE lim_dia_2  ***
69      !!   
70      !! ** Purpose : Computation and outputs on file ice.evolu
71      !!      the temporal evolution of some key variables
72      !!-------------------------------------------------------------------
73      INTEGER, INTENT(in) ::   kt     ! number of iteration
74      !!
75      INTEGER  ::   jv,ji, jj   ! dummy loop indices
76      INTEGER  ::   nv          ! indice of variable
77      REAL(wp) ::   zarea    , zldarea  ,    &  ! sea-ice and leads area
78         &          zextent15, zextent85,    &  ! sea-ice extent (15% and 85%)
79         &          zicevol  , zsnwvol  ,    &  ! sea-ice and snow volume volume
80         &          zicespd                     ! sea-ice velocity
81      REAL(wp), DIMENSION(jpinfmx) ::   vinfor  ! temporary working space
82      !!-------------------------------------------------------------------
83
84      IF( kt == nit000 )   CALL lim_dia_init_2  ! initialisation of ice_evolu file     
85
86      ! computation of key variables at each time step   
87
88      nv = 1 
89      vinfor(nv) = REAL( kt + nfice - 1 )
90      nv = nv + 1
91      vinfor(nv) = nyear
92 
93      DO jv = nbvt + 1, nvinfo
94         vinfor(jv) = 0.e0
95      END DO
96
97      zextent15 = 0.e0
98      zextent85 = 0.e0
99      ! variables in northern Hemis
100      DO jj = njeq, jpjm1
101         DO ji = fs_2, fs_jpim1   ! vector opt.
102            IF( tms(ji,jj) == 1 ) THEN
103               zarea = ( 1.0 - frld(ji,jj) ) * aire(ji,jj)
104               IF (frld(ji,jj) <= 0.15 ) zextent15 = aire(ji,jj)   
105               IF (frld(ji,jj) <= 0.85 ) zextent85 = aire(ji,jj)   
106               zldarea = zarea   / MAX( ( 1 - frld(ji,jj) ) , epsi06 )
107               zicevol = zarea   * hicif(ji,jj)
108               zsnwvol = zarea   * hsnif(ji,jj)
109               zicespd = zicevol * ( u_ice(ji,jj) * u_ice(ji,jj)   &
110                  &                + v_ice(ji,jj) * v_ice(ji,jj) )
111               vinfor(nv+ 1) = vinfor(nv+ 1) + zarea
112               vinfor(nv+ 3) = vinfor(nv+ 3) + zextent15
113               vinfor(nv+ 5) = vinfor(nv+ 5) + zextent85
114               vinfor(nv+ 7) = vinfor(nv+ 7) + zldarea
115               vinfor(nv+ 9) = vinfor(nv+ 9) + zicevol
116               vinfor(nv+11) = vinfor(nv+11) + zsnwvol
117               vinfor(nv+13) = vinfor(nv+13) + zicespd
118            ENDIF
119         END DO
120      END DO
121      vinfor(nv+13) = SQRT( vinfor(nv+13) / MAX( vinfor(nv+9) , epsi06 ) )
122
123
124      ! variables in southern Hemis
125       nv = nv + 1
126       DO jj = 2, njeqm1
127          DO ji = fs_2, fs_jpim1   ! vector opt.
128             IF( tms(ji,jj) == 1 ) THEN
129                zarea = ( 1.0 - frld(ji,jj) ) * aire(ji,jj)
130                IF (frld(ji,jj) <= 0.15 ) zextent15 = aire(ji,jj)   
131                IF (frld(ji,jj) <= 0.85 ) zextent85 = aire(ji,jj)   
132                zldarea = zarea   / MAX( ( 1 - frld(ji,jj) ) , epsi06 )
133                zicevol = zarea   * hicif(ji,jj)
134                zsnwvol = zarea   * hsnif(ji,jj)
135                zicespd = zicevol * ( u_ice(ji,jj) * u_ice(ji,jj)   &
136                   &                + v_ice(ji,jj) * v_ice(ji,jj) )
137                vinfor(nv+ 1) = vinfor(nv+ 1) + zarea
138                vinfor(nv+ 3) = vinfor(nv+ 3) + zextent15
139                vinfor(nv+ 5) = vinfor(nv+ 5) + zextent85
140                vinfor(nv+ 7) = vinfor(nv+ 7) + zldarea
141                vinfor(nv+ 9) = vinfor(nv+ 9) + zicevol
142                vinfor(nv+11) = vinfor(nv+11) + zsnwvol
143                vinfor(nv+13) = vinfor(nv+13) + zicespd
144             ENDIF
145          END DO
146       END DO
147       vinfor(nv+13) = SQRT( vinfor(nv+13) / MAX( vinfor(nv+9) , epsi06 ) )   
148
149       !  Accumulation before averaging
150       DO jv = 1, nvinfo
151          vinfom(jv) = vinfom(jv) + vinfor(jv)
152       END DO
153       naveg = naveg + 1 
154   
155       ! oututs on file ice_evolu   
156       IF( MOD( kt + nfice - 1, ninfo ) == 0 ) THEN
157          WRITE(numevo_ice,fmtw) ( titvar(jv), vinfom(jv)/naveg, jv = 1, nvinfo )
158          naveg = 0
159          DO jv = 1, nvinfo
160             vinfom(jv) = 0.e0
161          END DO
162       ENDIF
163       !
164    END SUBROUTINE lim_dia_2
165 
166
167    SUBROUTINE lim_dia_init_2
168       !!-------------------------------------------------------------------
169       !!                  ***  ROUTINE lim_dia_init_2  ***
170       !!             
171       !! ** Purpose : Preparation of the file ice_evolu for the output of
172       !!      the temporal evolution of key variables
173       !!
174       !! ** input   : Namelist namicedia
175       !!-------------------------------------------------------------------
176       CHARACTER(len=jpchinf) ::   titinf
177       CHARACTER(len=50)      ::   clname 
178       INTEGER  ::   jv            ! dummy loop indice
179       INTEGER  ::   ntot , ndeb , irecl
180       INTEGER  ::   nv            ! indice of variable
181       REAL(wp) ::   zxx0, zxx1    ! temporary scalars
182
183       NAMELIST/namicedia/fmtinf, nfrinf, ninfo, ntmoy
184       !!-------------------------------------------------------------------
185
186       ! Read Namelist namicedia
187       REWIND ( numnam_ice )
188       READ   ( numnam_ice  , namicedia )
189       
190       IF(lwp) THEN
191          WRITE(numout,*)
192          WRITE(numout,*) 'lim_dia_init_2 : ice parameters for ice diagnostics '
193          WRITE(numout,*) '~~~~~~~~~~~~~~'
194          WRITE(numout,*) '   format of the output values                                 fmtinf = ', fmtinf
195          WRITE(numout,*) '   number of variables written in one line                     nfrinf = ', nfrinf 
196          WRITE(numout,*) '   Instantaneous values of ice evolution or averaging          ntmoy  = ', ntmoy
197          WRITE(numout,*) '   frequency of ouputs on file ice_evolu in case of averaging  ninfo  = ', ninfo
198       ENDIF
199
200       ! masked grid cell area
201       aire(:,:) = area(:,:) * tms(:,:)
202
203       ! Titles of ice key variables :
204       nv = 1
205       titvar(nv) = 'NoIt'  ! iteration number
206       nv = nv + 1
207       titvar(nv) = 'T yr'  ! time step in years
208       
209       nbvt = nv - 1
210
211       nv = nv + 1   ;   titvar(nv) = 'AEFN' ! sea ice area in the northern Hemisp.(10^12 km2)
212       nv = nv + 1   ;   titvar(nv) = 'AEFS' ! sea ice area in the southern Hemisp.(10^12 km2)
213       nv = nv + 1   ;   titvar(nv) = 'A15N'  ! sea ice extent (15%) in the northern Hemisp.(10^12 km2)
214       nv = nv + 1   ;   titvar(nv) = 'A15S'  ! sea ice extent (15%) in the southern Hemisp.(10^12 km2)
215       nv = nv + 1   ;   titvar(nv) = 'A85N'  ! sea ice extent (85%) in the northern Hemisp.(10^12 km2)
216       nv = nv + 1   ;   titvar(nv) = 'A85S'  ! sea ice extent (85%) in the southern Hemisp.(10^12 km2)
217       nv = nv + 1   ;   titvar(nv) = 'ALEN'  ! leads area in the northern Hemisp.(10^12 km2)
218       nv = nv + 1   ;   titvar(nv) = 'ALES'  ! leads area in the southern Hemisp.(10^12 km2)
219       nv = nv + 1   ;   titvar(nv) = 'VOLN'  ! sea ice volume in the northern Hemisp.(10^3 km3)
220       nv = nv + 1   ;   titvar(nv) = 'VOLS'  ! sea ice volume in the southern Hemisp.(10^3 km3)
221       nv = nv + 1   ;   titvar(nv) = 'VONN'  ! snow volume over sea ice in the northern Hemisp.(10^3 km3)
222       nv = nv + 1   ;   titvar(nv) = 'VONS'  ! snow volume over sea ice in the southern Hemisp.(10^3 km3)
223       nv = nv + 1   ;   titvar(nv) = 'ECGN'  ! mean sea ice velocity in the northern Hemisp.(m/s)
224       nv = nv + 1   ;   titvar(nv) = 'ECGS'  ! mean sea ice velocity in the southern Hemisp.(m/s)
225
226       nvinfo = nv
227
228       ! Definition et Ecriture de l'entete : nombre d'enregistrements
229       ndeb   = ( nit000 - 1 + nfice - 1 ) / ninfo
230       IF( nit000 - 1 + nfice == 1 ) ndeb = -1
231
232       nferme = ( nitend + nfice - 1 ) / ninfo ! nit000 - 1 + nfice - 1 + nitend - nit000 + 1
233       ntot   = nferme - ndeb
234       ndeb   = ninfo * ( 1 + ndeb )
235       nferme = ninfo * nferme
236
237       ! definition of formats
238       WRITE( fmtw  , '(A,I3,A2,I1,A)' )  '(', nfrinf, '(A', jpchsep, ','//fmtinf//'))'
239       WRITE( fmtr  , '(A,I3,A,I1,A)'  )  '(', nfrinf, '(', jpchsep, 'X,'//fmtinf//'))'
240       WRITE( fmtitr, '(A,I3,A,I1,A)'  )  '(', nvinfo, 'A', jpchinf, ')'
241
242       ! opening  "ice_evolu" file
243       clname = 'ice_evolu' 
244       irecl = ( jpchinf + 1 ) * nvinfo 
245       CALL ctlopn( numevo_ice, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL',   &
246          &           irecl, numout, lwp, 1 )
247
248       !- ecriture de 2 lignes d''entete :
249       WRITE(numevo_ice,1000) fmtr, fmtw, fmtitr, nvinfo, ntot, 0, nfrinf
250       zxx0 = 0.001 * REAL( ninfo )
251       zxx1 = 0.001 * REAL( ndeb  )
252       WRITE(numevo_ice,1111) REAL(jpchinf), 0., zxx1, zxx0, 0., 0., 0
253
254       !- ecriture de 2 lignes de titre :
255       WRITE(numevo_ice,'(A,I8,A,I8,A,I5)')                 &
256          'Evolution chronologique - Experience '//cexper   &
257          //'   de', ndeb, ' a', nferme, ' pas', ninfo
258       WRITE(numevo_ice,fmtitr) ( titvar(jv), jv = 1, nvinfo )
259
260
261       !--preparation de "titvar" pour l''ecriture parmi les valeurs numeriques :
262       DO  jv = 2 , nvinfo
263          titinf     = titvar(jv)(:jpchinf)
264          titvar(jv) = '  '//titinf
265       END DO
266
267       !--Initialisation of the arrays for the accumulation
268       DO  jv = 1, nvinfo
269          vinfom(jv) = 0.
270       END DO
271       naveg = 0
272
2731000   FORMAT( 3(A20),4(1x,I6) )
2741111   FORMAT( 3(F7.1,1X,F7.3,1X),I3,A ) 
275      !
276    END SUBROUTINE lim_dia_init_2
277
278#else
279   !!----------------------------------------------------------------------
280   !!   Default option :                           NO LIM 2.0 sea-ice model
281   !!----------------------------------------------------------------------
282CONTAINS
283   SUBROUTINE lim_dia_2         ! Empty routine
284   END SUBROUTINE lim_dia_2
285#endif
286
287   !!======================================================================
288END MODULE limdia_2
Note: See TracBrowser for help on using the repository browser.