Changeset 98 for LMDZ.3.3/branches/rel-LF/libf/phylmd/oasis.F

Timestamp:

Jul 5, 2000, 4:58:04 PM (24 years ago)

Author:

lmdzadmin

Message:

Interface avec les differentes surface, version de travail.LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/oasis.F (modified) (5 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/oasis.F

@@ -386 +386 @@
       END
 
-      SUBROUTINE fromcpl(kt, imjm, sst, gla)
+      SUBROUTINE fromcpl(kt, imjm, sst,sic, alb_sst, alb_sic )
       IMPLICIT none
 c
@@ -399 +399 @@
       INTEGER imjm, kt
       REAL sst(imjm)          ! sea-surface-temperature
-      REAL gla(imjm)          ! sea-ice
+      REAL alb_sst(imjm)  ! open sea albedo
+      REAL sic(imjm)      ! sea ice cover
+      REAL alb_sic(imjm)  ! sea ice albedo
+
 c
       INTEGER nuout             ! listing output unit
@@ -447 +450 @@
      $          nuout)
             IF (jf.eq.2)
-     $          CALL locread(cl_read(jf), gla, imjm, nuread, iflag,
+     $          CALL locread(cl_read(jf), sic, imjm, nuread, iflag,
+     $          nuout)
+            IF (jf.eq.3)
+     $          CALL locread(cl_read(jf), alb_sst, imjm, nuread, iflag,
+     $          nuout)
+            IF (jf.eq.4)
+     $          CALL locread(cl_read(jf), alb_sic, imjm, nuread, iflag,
      $          nuout)
             CLOSE (nuread)
@@ -474 +483 @@
 c
           CALL SIPC_Read_Model(index, imjm, cmodinf, 
-     $              cljobnam_r,infos, gla)
+     $              cljobnam_r,infos, sic)
+c         Index of open sea albedo in total number of fields jpfldo2a: 
+          index = 3
+c
+          CALL SIPC_Read_Model(index, imjm, cmodinf, 
+     $              cljobnam_r,infos, alb_sst)
+c         Index of sea-ice albedo in total number of fields jpfldo2a: 
+          index = 4
+c
+          CALL SIPC_Read_Model(index, imjm, cmodinf, 
+     $              cljobnam_r,infos, alb_sic)
 c
 c
@@ -482 +501 @@
 c exchanges from ocean=CPL to atmosphere
 c
-          DO jf=1,jpfldo2a
-            IF (jf.eq.1) CALL CLIM_Import (cl_read(jf) , kt, sst, info)
-            IF (jf.eq.2) CALL CLIM_Import (cl_read(jf) , kt, gla, info)
-            IF ( info .NE. CLIM_Ok) THEN
+        DO jf=1,jpfldo2a
+          IF (jf.eq.1) CALL CLIM_Import (cl_read(jf) , kt, sst, info)
+          IF (jf.eq.2) CALL CLIM_Import (cl_read(jf) , kt, sic, info)
+         IF (jf.eq.3) CALL CLIM_Import (cl_read(jf) , kt, alb_sst, info)
+         IF (jf.eq.4) CALL CLIM_Import (cl_read(jf) , kt, alb_sic, info)
+          IF ( info .NE. CLIM_Ok) THEN
                 WRITE(nuout,*)'Pb in reading ', cl_read(jf), jf
                 WRITE(nuout,*)'Couplage kt is = ',kt