Changeset 149 for trunk/NEMO

Timestamp:

2004-10-19T15:32:08+02:00 (20 years ago)

Author:

opalod

Message:

CL + CT: BUGFIX094: Add "USE flx_oce" module, the computation of fr1_i0 and fr2_i0 (which was missing in couple mode) and syntax correction

File:

• trunk/NEMO/OPA_SRC/SBC/flx_coupled_ice.h90 (modified) (8 diffs)

trunk/NEMO/OPA_SRC/SBC/flx_coupled_ice.h90

@@ -33 +33 @@
       !!----------------------------------------------------------------------
       !! * Modules used
-      USE ioipsl
+      USE ioipsl               ! NetCDF IPSL library 
       USE ice_oce
-      USE cpl_oce
+      USE cpl_oce              ! coupled ocean-atmosphere variables
+      USE flx_oce              ! sea-ice/ocean forcings variables 
 
       !! * arguments
@@ -42 +43 @@
       !! * Local declarations
       INTEGER :: ji, jj, jf
-      INTEGER :: itm1,isize,iflag,icpliter
+      INTEGER :: itm1,isize,iflag
+!      INTEGER :: icpliter
       INTEGER :: info, inuread, index
       REAL(wp) ::   zfacflx,zfacwat
-      REAL(wp) ::   ztgel,zice
       REAL(wp) ::   znsolc (jpiglo,jpjglo),zqsrc (jpiglo,jpjglo)
       REAL(wp) ::   zrunoff(jpiglo,jpjglo),zec   (jpiglo,jpjglo)
@@ -55 +56 @@
       REAL(wp) ::   zrunriv (jpiglo,jpjglo),zcalving(jpiglo,jpjglo)
       REAL(wp) ::   zevap (jpiglo,jpjglo)
+      REAL(wp) ::   zcatm1 (jpiglo,jpjglo)            ! cloud fraction 
       CHARACTER (len=80) ::   clcplfnam
       REAL(wp) ::   zjulian
@@ -60 +62 @@
       ! Addition for SIPC CASE
       CHARACTER (len=3) ::   clmodinf       ! Header or not
-      CHARACTER (len=3) ::   cljobnam_r    ! Experiment name in the field brick, if any 
-      INTEGER ::   infos(3)          ! infos in the field brick, if any
+!      CHARACTER (len=3) ::   cljobnam_r    ! Experiment name in the field brick, if any 
+!      INTEGER ::   infos(3)          ! infos in the field brick, if any
       !!---------------------------------------------------------------------
 
@@ -205 +207 @@
             END DO
          END DO
+
+! Since cloud cover catm not transmitted from atmosphere, init =0. 
+ 
+         catm(:, :) =0.
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+            zcatm1(ji,jj) = 1.0    - catm  (ji,jj)  !  fractional cloud cover
+            END DO
+         END DO
+
+         !  fraction of net shortwave radiation which is not absorbed in the 
+         !  thin surface layer and penetrates inside the ice cover 
+         !  ( Maykut and Untersteiner, 1971 ; Elbert anbd Curry, 1993 )
+         !------------------------------------------------------------------
+         DO jj = 1, nlcj
+            DO ji = 1, nlci
+            fr1_i0(ji,jj) = 0.18  * zcatm1(ji,jj) + 0.35 * catm(ji,jj) 
+            fr2_i0(ji,jj) = 0.82  * zcatm1(ji,jj) + 0.65 * catm(ji,jj)
+            END DO
+         END DO
+
          ! copy in the subdomain
      
@@ -224 +247 @@
                sprecip(ji,jj) =  zpsol    ( mig(ji), mjg(jj)) * tmask(ji,jj,1) * zfacwat
                !  7: runoff      (positive upward)
-               srunoff(ji,jj) = ( zruncot ( mig(ji), mjg(jj))   &
+               rrunoff(ji,jj) = ( zruncot ( mig(ji), mjg(jj))   &
                   &              +  zrunriv ( mig(ji), mjg(jj)) ) * tmask(ji,jj,1) * zfacwat
                !  8: Derivative of non solar heat flux on sea ice
@@ -230 +253 @@
                !  13: Iceberg calving (positive upward)
                calving(ji,jj) =  zcalving ( mig(ji), mjg(jj)) * tmask(ji,jj,1) * zfacwat
+               !  1st part of the fraction of sol. rad.  which penetrate inside
+               !  the ice cover
+               fr1_i0(ji,jj)  = fr1_i0(mig(ji), mjg(jj)) * tmask(ji,jj,1)
+               ! 2nd part of the fraction of sol. rad.  which penetrate inside
+               ! the ice cover
+               fr2_i0(ji,jj)  = fr2_i0(mig(ji), mjg(jj)) * tmask(ji,jj,1)
               END DO
            END DO
@@ -240 +269 @@
          CALL lbc_lnk( tprecip , 'T', 1. )
          CALL lbc_lnk( sprecip , 'T', 1. )
-         CALL lbc_lnk( srunoff , 'T', 1. )
+         CALL lbc_lnk( rrunoff , 'T', 1. )
          CALL lbc_lnk( dqns_ice, 'T', 1. )
          CALL lbc_lnk( calving , 'T', 1. )
+         CALL lbc_lnk( fr1_i0  , 'T', 1. )
+         CALL lbc_lnk( fr2_i0  , 'T', 1. )
 
       ENDIF