Changeset 3557 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zche.F90

Timestamp:

2012-11-14T18:04:49+01:00 (11 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB minor modifications for PISCES when using AGRIF

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zche.F90 (modified) (3 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zche.F90

@@ -132 +132 @@
    REAL(wp) ::   ox10   = -3.11680e-7 
 
-   REAL(wp), DIMENSION(5)  :: devk1, devk2, devk3, devk4, devk5   ! coeff. for seawater pressure correction 
-   !                                                              ! (millero 95)
-   DATA devk1 / -25.5    , -15.82    , -29.48  , -25.60     , -48.76    /   
-   DATA devk2 / 0.1271   , -0.0219   , 0.1622  , 0.2324     , 0.5304    /   
-   DATA devk3 / 0.       , 0.        , 2.608E-3,  -3.6246E-3, 0.        /   
-   DATA devk4 / -3.08E-3 , 1.13E-3   , -2.84E-3, -5.13E-3   , -11.76E-3 /   
-   DATA devk5 / 0.0877E-3, -0.1475E-3,  0.     , 0.0794E-3  , 0.3692E-3 /
+   !                                    ! coeff. for seawater pressure correction : millero 95
+   !                                    ! AGRIF doesn't like the DATA instruction
+   REAL(wp) :: devk11  = -25.5
+   REAL(wp) :: devk12  = -15.82
+   REAL(wp) :: devk13  = -29.48
+   REAL(wp) :: devk14  = -25.60
+   REAL(wp) :: devk15  = -48.76
+   !
+   REAL(wp) :: devk21  = 0.1271
+   REAL(wp) :: devk22  = -0.0219
+   REAL(wp) :: devk23  = 0.1622
+   REAL(wp) :: devk24  = 0.2324
+   REAL(wp) :: devk25  = 0.5304
+   !
+   REAL(wp) :: devk31  = 0.
+   REAL(wp) :: devk32  = 0.
+   REAL(wp) :: devk33  = 2.608E-3
+   REAL(wp) :: devk34  = -3.6246E-3
+   REAL(wp) :: devk35  = 0.
+   !
+   REAL(wp) :: devk41  = -3.08E-3
+   REAL(wp) :: devk42  = 1.13E-3
+   REAL(wp) :: devk43  = -2.84E-3
+   REAL(wp) :: devk44  = -5.13E-3
+   REAL(wp) :: devk45  = -11.76E-3
+   !
+   REAL(wp) :: devk51  = 0.0877E-3
+   REAL(wp) :: devk52  = -0.1475E-3     
+   REAL(wp) :: devk53  = 0.
+   REAL(wp) :: devk54  = 0.0794E-3      
+   REAL(wp) :: devk55  = 0.3692E-3      
 
    !!* Substitution
@@ -310 +334 @@
                !        (CF. BROECKER ET AL., 1982)
 
-               zbuf1  = -(devk1(1)+devk2(1)*ztc+devk3(1)*ztc*ztc)
-               zbuf2  = 0.5*(devk4(1)+devk5(1)*ztc)
+               zbuf1  = -     ( devk11 + devk21 * ztc + devk31 * ztc * ztc )
+               zbuf2  = 0.5 * ( devk41 + devk51 * ztc )
                ak13(ji,jj,jk) = zak1 * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 )
 
-               zbuf1  =     - ( devk1(2) + devk2(2) * ztc + devk3(2) * ztc * ztc )
-               zbuf2  = 0.5 * ( devk4(2) + devk5(2) * ztc )
+               zbuf1  =     - ( devk12 + devk22 * ztc + devk32 * ztc * ztc )
+               zbuf2  = 0.5 * ( devk42 + devk52 * ztc )
                ak23(ji,jj,jk) = zak2 * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 )
 
-               zbuf1  =     - ( devk1(3) + devk2(3) * ztc + devk3(3) * ztc * ztc )
-               zbuf2  = 0.5 * ( devk4(3) + devk5(3) * ztc )
+               zbuf1  =     - ( devk13 + devk23 * ztc + devk33 * ztc * ztc )
+               zbuf2  = 0.5 * ( devk43 + devk53 * ztc )
                akb3(ji,jj,jk) = zakb * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 )
 
-               zbuf1  =     - ( devk1(4) + devk2(4) * ztc + devk3(4) * ztc * ztc )
-               zbuf2  = 0.5 * ( devk4(4) + devk5(4) * ztc )
+               zbuf1  =     - ( devk14 + devk24 * ztc + devk34 * ztc * ztc )
+               zbuf2  = 0.5 * ( devk44 + devk54 * ztc )
                akw3(ji,jj,jk) = zakw * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 )
 
@@ -330 +354 @@
                !        AS FUNCTION OF PRESSURE FOLLOWING MILLERO
                !        (P. 1285) AND BERNER (1976)
-               zbuf1  =     - ( devk1(5) + devk2(5) * ztc + devk3(5) * ztc * ztc )
-               zbuf2  = 0.5 * ( devk4(5) + devk5(5) * ztc )
+               zbuf1  =     - ( devk15 + devk25 * ztc + devk35 * ztc * ztc )
+               zbuf2  = 0.5 * ( devk45 + devk55 * ztc )
                aksp(ji,jj,jk) = zaksp1 * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 )
-
 
                ! TOTAL BORATE CONCENTR. [MOLES/L]