Changeset 313 for LMDZ.3.3/branches/rel-LF

Timestamp:

Jan 25, 2002, 11:59:34 AM (23 years ago)

Author:

lmdzadmin

Message:

Prise en compte de la fonte des glaciers antarctiques pour le couplé
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90 (modified) (4 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -44 +44 @@
 !!$PB
   REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: tmp_rriv, tmp_rcoa
+!! pour simuler la fonte des glaciers antarctiques
+  REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: coeff_iceberg
+  real, save                              :: surf_maille 
+  real, save                              :: cte_flux_iceberg = 6.3e7
+  integer, save                           :: num_antarctic = 1
 !!$
   CONTAINS
@@ -1123 +1128 @@
 !   itime        numero du pas de temps
 !   iim, jjm     nbres de pts de grille
-!   dtime        pas de tempsde la physique
+!   dtime        pas de temps de la physique
 !   klon         nombre total de points de grille
 !   nisurf       index de la surface a traiter (1 = sol continental)
@@ -1343 +1348 @@
     CALL histsync(nidcs)
 
+! pour simuler la fonte des glaciers antarctiques
+! 
+    surf_maille = (4. * rpi * ra**2) / (iim * (jjm +1))
+    ALLOCATE(coeff_iceberg(iim,jjm+1), stat=error)
+    if (error /= 0) then
+      abort_message='Pb allocation variable coeff_iceberg'
+      call abort_gcm(modname,abort_message,1)
+    endif
+    open (12,file='flux_iceberg',form='formatted',status='old')
+    read (12,*) coeff_iceberg
+    close (12)
+    num_antarctic = max(1, count(coeff_iceberg > 0))
+    
     first_appel = .false.
   endif ! fin if (first_appel)
@@ -1544 +1562 @@
       endwhere
 !
+! pour simuler la fonte des glaciers antarctiques
+!
+        wri_rain = wri_rain      &
+      &     + coeff_iceberg * cte_flux_iceberg / (num_antarctic * surf_maille)
+
+!
 ! on passe les coordonnées de la grille
 !