Changeset 146 for branches/iLoveclim/SOURCES/spinup_mod.f90

Timestamp:

10/20/17 09:31:39 (7 years ago)

Author:

aquiquet

Message:

Grisli-iLoveclim branch: merged to trunk at revision 145

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/spinup_mod.f90 (modified) (2 diffs)

branches/iLoveclim/SOURCES/spinup_mod.f90

@@ -102 +102 @@
 case default
   print*,'type_vitbil valeur invalide dans spinup_vitbil'
+  stop
 end select
 
@@ -257 +258 @@
     ddbx(:,:)=0.
     ddby(:,:)=0.
-
+    gzmx(:,:)=.true.
+    gzmy(:,:)=.true.
+    flgzmx(:,:)=.false.
+    flgzmy(:,:)=.false.
+    fleuvemx(:,:)=.false.
+    fleuvemy(:,:)=.false.
+    
+    do j=2,ny
+      do i=2,nx
+          if (flot(i,j).and.(flot(i-1,j))) then
+            flotmx(i,j)=.true.
+            flgzmx(i,j)=.true.
+            gzmx(i,j)=.false.
+          end if
+          if (flot(i,j).and.(flot(i,j-1))) then
+            flotmy(i,j)=.true.
+            flgzmy(i,j)=.true.
+            gzmy(i,j)=.false.
+          end if
+      end do
+    end do   
+    
+    where (coefmxbmelt(:,:).le.0.)
+      gzmx(:,:)=.false.
+    endwhere
+    where (coefmybmelt(:,:).le.0.)
+      gzmy(:,:)=.false.
+    endwhere
+    flgzmx(:,:) = flgzmx(:,:) .or. gzmx(:,:)
+    flgzmy(:,:) = flgzmy(:,:) .or. gzmy(:,:)
+    fleuvemx(:,:)=gzmx(:,:)
+    fleuvemy(:,:)=gzmy(:,:)
+    
     if (itracebug.eq.1)  call tracebug(' Subroutine  calc_coef_vitbil')
     call slope_surf