Changeset 1124 for LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/cpl_mod.F90

Timestamp:

Mar 17, 2009, 10:53:39 AM (16 years ago)

Author:

jghattas

Message:

Modification pour couplage conservatif avec NEMO : passage des champs totaux et champs sur la glace au lieu de passage des champs sur l'ocean et champs sur la glace.

Arnaud Caubel

File:

• LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/cpl_mod.F90 (modified) (4 diffs)

LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/cpl_mod.F90

@@ -954 +954 @@
     tab_flds(:,:,7)  = cpl_windsp2D(:,:)
     tab_flds(:,:,8)  = cpl_sols2D(:,:,2)
-    tab_flds(:,:,9)  = cpl_sols2D(:,:,1)
     tab_flds(:,:,10) = cpl_nsol2D(:,:,2)
-    tab_flds(:,:,11) = cpl_nsol2D(:,:,1)
     tab_flds(:,:,12) = cpl_fder2D(:,:,2)
     
@@ -962 +960 @@
        tab_flds(:,:,18) = cpl_rriv2D(:,:) + cpl_rcoa2D(:,:)
     ELSE IF (version_ocean=='opa8') THEN
+       tab_flds(:,:,9)  = cpl_sols2D(:,:,1)
+       tab_flds(:,:,11) = cpl_nsol2D(:,:,1)
+
        tab_flds(:,:,13) = cpl_evap2D(:,:,2)
        tab_flds(:,:,14) = cpl_evap2D(:,:,1)
@@ -1030 +1031 @@
 
     IF (version_ocean=='nemo') THEN
-       tab_flds(:,:,13) = cpl_evap2D(:,:,2) - (cpl_rain2D(:,:,2) + cpl_snow2D(:,:,2))
-
+
+       tab_flds(:,:,9)  = 0.0
+       tab_flds(:,:,11) = 0.0
+       tab_flds(:,:,13) = 0.0
        tab_flds(:,:,14) = 0.0
+       tab_flds(:,:,15) = 0.0
        tmp_taux(:,:)    = 0.0
        tmp_tauy(:,:)    = 0.0
+       
        ! For all valid grid cells containing some fraction of ocean or sea-ice
        WHERE ( deno(:,:) /= 0 )
@@ -1041 +1046 @@
            tmp_tauy = cpl_tauy2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
               cpl_tauy2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
+           tab_flds(:,:,9) = cpl_sols2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+              cpl_sols2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
+           tab_flds(:,:,11) = cpl_nsol2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+              cpl_nsol2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
+           tab_flds(:,:,13) = cpl_rain2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+               cpl_rain2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
            tab_flds(:,:,14) = cpl_snow2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
               cpl_snow2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
+           tab_flds(:,:,15) = cpl_evap2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+              cpl_evap2D(:,:,2)  * pctsrf2D(:,:,is_sic) / deno(:,:)
        ENDWHERE
 
-       WHERE (pctsrf2D(:,:,is_sic) == 0)
-          cpl_evap2D(:,:,2) = 0.0
-       END WHERE
-
-       tab_flds(:,:,14) = cpl_evap2D(:,:,2) - tab_flds(:,:,14)
-       tab_flds(:,:,15) = cpl_evap2D(:,:,1) - (cpl_rain2D(:,:,1) + cpl_snow2D(:,:,1))
        tab_flds(:,:,16) = cpl_evap2D(:,:,2)