Changeset 830

Timestamp:

2008-03-06T19:46:43+01:00 (16 years ago)

Author:

ctlod

Message:

dev_002_LIM : add declaration of Atmospheric drag coefficients over sea-ice & ocean cai/cao

File:

• branches/dev_002_LIM/NEMO/LIM_SRC_3/iceini.F90 (modified) (3 diffs)

branches/dev_002_LIM/NEMO/LIM_SRC_3/iceini.F90

@@ -41 +41 @@
       numit                  !: iteration number
    REAL(wp), PUBLIC  ::   &  !:
-      hsndif = 0.e0 ,     &  !: computation of temp. in snow (0) or not (9999)
-      hicdif = 0.e0 ,     &  !: computation of temp. in ice (0) or not (9999)
-      tpstot                 !: time of the run in seconds
+      hsndif = 0.e0    ,  &  !: computation of temp. in snow (0) or not (9999)
+      hicdif = 0.e0    ,  &  !: computation of temp. in ice (0) or not (9999)
+      tpstot           ,  &  !: time of the run in seconds
+      cai    = 1.40e-3 ,  &  !: atmospheric drag over sea ice
+      cao    = 1.00e-3       !: atmospheric drag over ocean
    REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:
       acrit  = (/ 1.e-06 , 1.e-06 /)    !: minimum fraction for leads in 
@@ -129 +131 @@
       !!   2.0  !  03-08 (C. Ethe)  Original code
       !!-------------------------------------------------------------------
-      NAMELIST/namicerun/ ln_limdyn, acrit, hsndif, hicdif, ln_nicep
+      NAMELIST/namicerun/ ln_limdyn, acrit, hsndif, hicdif, cai, cao, ln_nicep
       !!-------------------------------------------------------------------
 
@@ -143 +145 @@
          WRITE(numout,*) '   computation of temp. in snow (=0) or not (=9999) hsndif = ', hsndif
          WRITE(numout,*) '   computation of temp. in ice  (=0) or not (=9999) hicdif = ', hicdif
+         WRITE(numout,*) '   atmospheric drag over sea ice                           = ', cai
+         WRITE(numout,*) '   atmospheric drag over ocean                             = ', cao
          WRITE(numout,*) '   Several ice points in the ice or not in ocean.output = ', ln_nicep
       ENDIF