Changeset 146 for branches/iLoveclim/SOURCES/calving_frange.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/calving_frange.f90 (modified) (3 diffs)

branches/iLoveclim/SOURCES/calving_frange.f90

@@ -122 +122 @@
   mass_total=0d0
 !afq for CONSEAU  calvtest=0.
-  mass_total=0d0
 
 ! calcul du dhdt lagrangien
@@ -159 +158 @@
 
 !ifint:    if((front(i,j).gt.0).and.(front(i,j).lt.4)) then
-        ifint:    if(front(i,j).lt.4) then
+!cdc pb avec front        ifint:    if(front(i,j).lt.4) then
+        ifint:    if((H(i-1,j).lt.2.).or.(H(i+1,j).lt.2.).or.(H(i,j-1).lt.2.).or.(H(i,j+1).lt.2)) then ! si on est au bord avec test sur H
 ! ifint: le point doit avoir au - 1 voisin mais ne pas etre entouré de glace
 ! ce qui evite la formation des polynies dans les shelfs
@@ -491 +491 @@
                      S(i,j)=H(i,j)*(1.-ro/row) + sealevel
                      B(i,j)=S(i,j) - H(i,j)
+          ! ATTENTION ne pas mettre ice=0 sinon degradation bilan d'eau (bm et bmelt non comptabilises dans ce cas)
                    endif
                    !H(i,j)=1