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.
trcice_pisces.F90 in branches/UKMO/dev_r5107_mld_zint/NEMOGCM/NEMO/TOP_SRC/PISCES – NEMO

source: branches/UKMO/dev_r5107_mld_zint/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90 @ 5447

Last change on this file since 5447 was 5447, checked in by davestorkey, 9 years ago

Update UKMO/dev_r5107_mld_zint branch to revision 5442 of the trunk.

File size: 16.4 KB
Line 
1MODULE trcice_pisces
2   !!======================================================================
3   !!                         ***  MODULE trcice_pisces  ***
4   !! TOP :   initialisation of the PISCES biochemical model
5   !!======================================================================
6   !! History :  3.5  ! 2013    (M. Vancoppenolle, O. Aumont, G. Madec), original code
7   !! Comment ! probably not properly done when the second particle export
8   !! scheme (kriest) is used
9   !!----------------------------------------------------------------------
10#if defined key_pisces || defined key_pisces_reduced
11   !!----------------------------------------------------------------------
12   !!   'key_pisces'                                       PISCES bio-model
13   !!----------------------------------------------------------------------
14   !! trc_ice_pisces   : PISCES fake sea ice model setting
15   !!----------------------------------------------------------------------
16   USE par_trc         ! TOP parameters
17   USE par_pisces      ! PISCES parameters
18   USE oce_trc         ! Shared variables between ocean and passive tracers
19   USE trc             ! Passive tracers common variables
20   USE phycst          ! Ocean physics parameters
21   USE sms_pisces      ! PISCES Source Minus Sink variables
22   USE in_out_manager
23
24   IMPLICIT NONE
25   PRIVATE
26
27   PUBLIC   trc_ice_ini_pisces ! called by trcini.F90 module
28
29CONTAINS
30
31   SUBROUTINE trc_ice_ini_pisces
32      !!----------------------------------------------------------------------
33      !!                   ***  ROUTINE trc_ice_ini_pisces ***
34      !!
35      !! ** Purpose :   PISCES fake sea ice model setting
36      !!    Method  :   Assign prescribe values to tracer concentrations in sea ice
37      !!
38      !! For levitating sea ice, constant ocean tracer concentrations also have to be defined.
39      !! This is done specifically for Global, Arctic, Antarctic and Baltic regions
40      !!
41      !! Sea ice concentrations are by default prescribed as follows
42      !!  trc_i = zratio * trc_o
43      !!
44      !! This formulation is modulated by the namelist parameter trc_ice_ratio
45      !!
46      !! trc_ice_ratio  * betw 0 and 1: prescribed ice/ocean tracer concentration ratio
47      !!                * -1 => the ice-ocean tracer concentration ratio follows the
48      !!                         ice-ocean salinity ratio
49      !!                * -2 => no ice-ocean tracer concentration is used
50      !!                        instead, the tracer concentration in sea ice
51      !!                        is prescribed to trc_ice_prescr
52      !!
53      !! cn_trc_o  specifies which disinctions are made for prescribed tracer concentration
54      !!                * 'GL' use global ocean values making distinction for Baltic Sea only
55      !!                * 'AA' use Arctic/Antarctic contrasted values, + Baltic
56      !!
57      !!----------------------------------------------------------------------
58
59                                        !--- Dummy variables
60      REAL(wp), DIMENSION(jptra,2) &
61               ::  zratio            ! effective ice-ocean tracer cc ratio
62      REAL(wp), DIMENSION(2) :: zrs  ! ice-ocean salinity ratio, 1 - global, 2- Baltic
63      REAL(wp) :: zsice_bal          ! prescribed ice salinity in the Baltic
64      REAL(wp) :: zsoce_bal          ! prescribed ocean salinity in the Baltic
65      REAL(wp) :: zfeoce_glo         ! prescribed iron concentration in the global ocean
66      REAL(wp) :: zfeoce_bal         ! prescribed iron concentration in the global ocean
67      INTEGER  :: jn                 ! dummy loop index
68
69      !!----------------------------------------------------------------------
70
71      IF(lwp) WRITE(numout,*)
72      IF(lwp) WRITE(numout,*) ' trc_ice_ini_pisces: Prescribed sea ice biogeochemistry '
73      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~'
74
75      !--------------------------------------------
76      ! Initialize ocean prescribed concentrations
77      !--------------------------------------------
78      ! values taken from a 500 yr equilibrium run
79      ! used only in the levitating sea ice case with virtual salt / tracer
80      ! fluxes
81
82      !--- Global case
83      IF ( cn_trc_o(jpdic) == 'GL ' ) trc_o(:,:,jpdic) =  1.99e-3_wp 
84      IF ( cn_trc_o(jpdoc) == 'GL ' ) trc_o(:,:,jpdoc) =  2.04e-5_wp 
85      IF ( cn_trc_o(jptal) == 'GL ' ) trc_o(:,:,jptal) =  2.31e-3_wp 
86      IF ( cn_trc_o(jpoxy) == 'GL ' ) trc_o(:,:,jpoxy) =  2.47e-4_wp
87      IF ( cn_trc_o(jpcal) == 'GL ' ) trc_o(:,:,jpcal) =  1.04e-8_wp
88      IF ( cn_trc_o(jppo4) == 'GL ' ) trc_o(:,:,jppo4) =  5.77e-7_wp / po4r 
89      IF ( cn_trc_o(jppoc) == 'GL ' ) trc_o(:,:,jppoc) =  1.27e-6_wp 
90#  if ! defined key_kriest
91      IF ( cn_trc_o(jpgoc) == 'GL ' ) trc_o(:,:,jpgoc) =  5.23e-8_wp 
92      IF ( cn_trc_o(jpbfe) == 'GL ' ) trc_o(:,:,jpbfe) =  9.84e-13_wp 
93#  else
94      IF ( cn_trc_o(jpnum) == 'GL ' ) trc_o(:,:,jpnum) = 0. ! could not get this value since did not use it
95#  endif
96      IF ( cn_trc_o(jpsil) == 'GL ' ) trc_o(:,:,jpsil) =  7.36e-6_wp 
97      IF ( cn_trc_o(jpdsi) == 'GL ' ) trc_o(:,:,jpdsi) =  1.07e-7_wp 
98      IF ( cn_trc_o(jpgsi) == 'GL ' ) trc_o(:,:,jpgsi) =  1.53e-8_wp
99      IF ( cn_trc_o(jpphy) == 'GL ' ) trc_o(:,:,jpphy) =  9.57e-8_wp
100      IF ( cn_trc_o(jpdia) == 'GL ' ) trc_o(:,:,jpdia) =  4.24e-7_wp
101      IF ( cn_trc_o(jpzoo) == 'GL ' ) trc_o(:,:,jpzoo) =  6.07e-7_wp
102      IF ( cn_trc_o(jpmes) == 'GL ' ) trc_o(:,:,jpmes) =  3.44e-7_wp
103      IF ( cn_trc_o(jpfer) == 'GL ' ) trc_o(:,:,jpfer) =  4.06e-10_wp
104      IF ( cn_trc_o(jpsfe) == 'GL ' ) trc_o(:,:,jpsfe) =  2.51e-11_wp
105      IF ( cn_trc_o(jpdfe) == 'GL ' ) trc_o(:,:,jpdfe) =  6.57e-12_wp
106      IF ( cn_trc_o(jpnfe) == 'GL ' ) trc_o(:,:,jpnfe) =  1.76e-11_wp
107      IF ( cn_trc_o(jpnch) == 'GL ' ) trc_o(:,:,jpnch) =  1.67e-7_wp
108      IF ( cn_trc_o(jpdch) == 'GL ' ) trc_o(:,:,jpdch) =  1.02e-7_wp
109      IF ( cn_trc_o(jpno3) == 'GL ' ) trc_o(:,:,jpno3) =  5.79e-6_wp / rno3 
110      IF ( cn_trc_o(jpnh4) == 'GL ' ) trc_o(:,:,jpnh4) =  3.22e-7_wp / rno3
111
112      !--- Arctic specificities (dissolved inorganic & DOM)
113      IF ( cn_trc_o(jpdic) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdic) =  1.98e-3_wp  ; END WHERE ; ENDIF
114      IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdoc) =  6.00e-6_wp  ; END WHERE ; ENDIF
115      IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jptal) =  2.13e-3_wp  ; END WHERE ; ENDIF
116      IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpoxy) =  3.65e-4_wp  ; END WHERE ; ENDIF
117      IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpcal) =  1.50e-9_wp  ; END WHERE ; ENDIF
118      IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppo4) =  4.09e-7_wp / po4r ; END WHERE ; ENDIF
119      IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppoc) =  4.05e-7_wp  ; END WHERE ; ENDIF
120#  if ! defined key_kriest
121      IF ( cn_trc_o(jpgoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgoc) =  2.84e-8_wp  ; END WHERE ; ENDIF
122      IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpbfe) =  7.03e-13_wp ; END WHERE ; ENDIF
123#  else
124      IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnum) =  0.00e-00_wp ; END WHERE ; ENDIF
125#  endif
126      IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsil) =  6.87e-6_wp  ; END WHERE ; ENDIF
127      IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdsi) =  1.73e-7_wp  ; END WHERE ; ENDIF
128      IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgsi) =  7.93e-9_wp  ; END WHERE ; ENDIF
129      IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpphy) =  5.25e-7_wp  ; END WHERE ; ENDIF
130      IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdia) =  7.75e-7_wp  ; END WHERE ; ENDIF
131      IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpzoo) =  3.34e-7_wp  ; END WHERE ; ENDIF
132      IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpmes) =  2.49e-7_wp  ; END WHERE ; ENDIF
133      IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpfer) =  1.43e-9_wp  ; END WHERE ; ENDIF
134      IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsfe) =  2.21e-11_wp ; END WHERE ; ENDIF
135      IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdfe) =  2.04e-11_wp ; END WHERE ; ENDIF
136      IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnfe) =  1.75e-11_wp ; END WHERE ; ENDIF
137      IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnch) =  1.46e-07_wp ; END WHERE ; ENDIF
138      IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdch) =  2.36e-07_wp ; END WHERE ; ENDIF
139      IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpno3) =  3.51e-06_wp / rno3 ; END WHERE ; ENDIF
140      IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnh4) =  6.15e-08_wp / rno3 ; END WHERE ; ENDIF
141
142      !--- Antarctic specificities (dissolved inorganic & DOM)
143      IF ( cn_trc_o(jpdic) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdic) =  2.20e-3_wp  ; END WHERE ; ENDIF
144      IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdoc) =  7.02e-6_wp  ; END WHERE ; ENDIF
145      IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jptal) =  2.37e-3_wp  ; END WHERE ; ENDIF
146      IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpoxy) =  3.42e-4_wp  ; END WHERE ; ENDIF
147      IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpcal) =  3.17e-9_wp  ; END WHERE ; ENDIF
148      IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jppo4) =  1.88e-6_wp / po4r  ; END WHERE ; ENDIF
149      IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jppoc) =  1.13e-6_wp  ; END WHERE ; ENDIF
150#  if ! defined key_kriest
151      IF ( cn_trc_o(jpgoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpgoc) =  2.89e-8_wp  ; END WHERE ; ENDIF
152      IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpbfe) =  5.63e-13_wp ; END WHERE ; ENDIF
153#  else
154      IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnum) =  0.00e-00_wp ; END WHERE ; ENDIF
155#  endif
156      IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpsil) =  4.96e-5_wp  ; END WHERE ; ENDIF
157      IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdsi) =  5.63e-7_wp  ; END WHERE ; ENDIF
158      IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpgsi) =  5.35e-8_wp  ; END WHERE ; ENDIF
159      IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpphy) =  8.10e-7_wp  ; END WHERE ; ENDIF
160      IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdia) =  5.77e-7_wp  ; END WHERE ; ENDIF
161      IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpzoo) =  6.68e-7_wp  ; END WHERE ; ENDIF
162      IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpmes) =  3.55e-7_wp  ; END WHERE ; ENDIF
163      IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpfer) =  1.62e-10_wp ; END WHERE ; ENDIF
164      IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpsfe) =  2.29e-11_wp ; END WHERE ; ENDIF
165      IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdfe) =  8.75e-12_wp ; END WHERE ; ENDIF
166      IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnfe) =  1.48e-11_wp ; END WHERE ; ENDIF
167      IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnch) =  2.02e-7_wp  ; END WHERE ; ENDIF
168      IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdch) =  1.60e-7_wp  ; END WHERE ; ENDIF
169      IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpno3) =  2.64e-5_wp / rno3  ; END WHERE ; ENDIF
170      IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnh4) =  3.39e-7_wp / rno3  ; END WHERE ; ENDIF
171
172      !--- Baltic Sea particular case for ORCA configurations
173      IF( cp_cfg == "orca" ) THEN            ! Baltic mask
174         WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    &
175                54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp )
176         trc_o(:,:,jpdic) = 1.14e-3_wp
177         trc_o(:,:,jpdoc) = 1.06e-5_wp
178         trc_o(:,:,jptal) = 1.16e-3_wp
179         trc_o(:,:,jpoxy) = 3.71e-4_wp
180         trc_o(:,:,jpcal) = 1.51e-9_wp
181         trc_o(:,:,jppo4) = 2.85e-9_wp / po4r
182         trc_o(:,:,jppoc) = 4.84e-7_wp
183#  if ! defined key_kriest
184         trc_o(:,:,jpgoc) = 1.05e-8_wp
185         trc_o(:,:,jpbfe) = 4.97e-13_wp
186#  else
187         trc_o(:,:,jpnum) = 0. ! could not get this value
188#  endif
189         trc_o(:,:,jpsil) = 4.91e-5_wp
190         trc_o(:,:,jpdsi) = 3.25e-7_wp
191         trc_o(:,:,jpgsi) = 1.93e-8_wp
192         trc_o(:,:,jpphy) = 6.64e-7_wp
193         trc_o(:,:,jpdia) = 3.41e-7_wp
194         trc_o(:,:,jpzoo) = 3.83e-7_wp
195         trc_o(:,:,jpmes) = 0.225e-6_wp
196         trc_o(:,:,jpfer) = 2.45e-9_wp
197         trc_o(:,:,jpsfe) = 3.89e-11_wp
198         trc_o(:,:,jpdfe) = 1.33e-11_wp
199         trc_o(:,:,jpnfe) = 2.62e-11_wp
200         trc_o(:,:,jpnch) = 1.17e-7_wp
201         trc_o(:,:,jpdch) = 9.69e-8_wp
202         trc_o(:,:,jpno3) = 5.36e-5_wp / rno3
203         trc_o(:,:,jpnh4) = 7.18e-7_wp / rno3
204         END WHERE
205      ENDIF ! cfg
206
207      !-----------------------------
208      ! Assign ice-ocean cc ratios
209      !-----------------------------
210      ! 0 means zero concentration in sea ice
211      ! 1 means same concentration in the sea ice as in the ocean
212
213      ! Ice ocean salinity ratio
214      zsoce_bal   = 4. ; zsice_bal   = 2. !! Baltic ocean and sea ice salinities
215      zrs(1) = sice / soce                !! ice-ocean salinity ratio, global case
216      zrs(2) = zsice_bal / zsoce_bal      !! ice-ocean salinity ratio, Baltic case
217
218      DO jn = jp_pcs0, jp_pcs1
219         IF ( trc_ice_ratio(jn) >= 0._wp )  zratio(jn,:) = trc_ice_ratio(jn)
220         IF ( trc_ice_ratio(jn) == -1._wp ) zratio(jn,:) = zrs(:)
221         IF ( trc_ice_ratio(jn) == -2._wp ) zratio(jn,:) = -9999.99_wp
222      END DO
223
224      !-------------------------------
225      ! Sea ice tracer concentrations
226      !-------------------------------
227      DO jn = jp_pcs0, jp_pcs1
228         !-- Everywhere but in the Baltic
229         IF ( trc_ice_ratio(jn) >= -1._wp ) THEN !! no prescribed concentration
230                                              !! (typically everything but iron)
231            trc_i(:,:,jn) = zratio(jn,1) * trc_o(:,:,jn) 
232         ELSE                                 !! prescribed concentration
233            trc_i(:,:,jn) = trc_ice_prescr(jn)
234         ENDIF
235       
236         !-- Baltic
237         IF( cp_cfg == "orca" ) THEN !! Baltic treated seperately for ORCA configs
238            IF ( trc_ice_ratio(jn) >= - 1._wp ) THEN !! no prescribed concentration
239                                                 !! (typically everything but iron)
240               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    &
241                      54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp )
242                     trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn) 
243               END WHERE
244            ELSE                                 !! prescribed tracer concentration in ice
245               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    &
246                   54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp )
247                     trc_i(:,:,jn) = trc_ice_prescr(jn)
248               END WHERE
249            ENDIF ! trc_ice_ratio
250         ENDIF
251      !
252      END DO ! jn
253
254   END SUBROUTINE trc_ice_ini_pisces
255
256#else
257   !!----------------------------------------------------------------------
258   !!   Dummy module                            No PISCES biochemical model
259   !!----------------------------------------------------------------------
260CONTAINS
261   SUBROUTINE trc_ice_ini_pisces         ! Empty routine
262   END SUBROUTINE trc_ice_ini_pisces
263#endif
264
265   !!======================================================================
266END MODULE trcice_pisces
Note: See TracBrowser for help on using the repository browser.