Changeset 101 for branches

Timestamp:

04/13/17 10:10:00 (7 years ago)

Author:

aquiquet

Message:

few updates for GRISLI-CLIO coupling.

Location:

branches/iLoveclim/SOURCES

Files:

• 3D-physique-gen_mod.f90 (modified) (1 diff)
• bmelt_clio_coupl_mod.f90 (modified) (4 diffs)

branches/iLoveclim/SOURCES/3D-physique-gen_mod.f90

@@ -325 +325 @@
   real,dimension(nx,ny) :: NEFFMX       !< pression effective '>'
   real,dimension(nx,ny) :: NEFFMY       !< pression effective '^'
-  real(kind=8),dimension(nx,ny),target :: RUNOF_OC ! Ablation at the ice shelf that is given to CLIO
+  real(kind=8),dimension(nx,ny),target :: RUNOF_OC ! Continental runoff (grounded part)
+  real(kind=8),dimension(nx,ny),target :: BMELT_OC ! Shelf runoff (floating part)
   real,dimension(nx,ny) :: TOBMX        !< cisaillement basal '>'
   real,dimension(nx,ny) :: TOBMY        !< cisaillement basal '^'

branches/iLoveclim/SOURCES/bmelt_clio_coupl_mod.f90

@@ -103 +103 @@
 ! In CLIO, when we don't have ocean, I assumed a very high bm
 !          In this case, in GRISLI, we put constant values (0.2)
-           if (bmshelfCLIO(i,j,20).lt.999d0) then! 
+           if (bmshelfCLIO(i,j,20).lt.99d0) then! 
               bmshelf(i,j) = real(bmshelfCLIO(i,j,20))
               bmgrz(i,j) = bmshelf(i,j) * bmgrz_fact
@@ -120 +120 @@
               do noc=20,2,-1
                  if ( locdepth .gt. z_CLIO(noc)+dz_CLIO(noc)/2. ) then
-                    if (bmshelfCLIO(i,j,noc+1).lt.999d0) then
+                    if (bmshelfCLIO(i,j,noc+1).lt.99d0) then
                        bmelt(i,j) = real(bmshelfCLIO(i,j,noc+1))
                     endif
@@ -156 +156 @@
 ! In CLIO, when we don't have ocean, I assumed a very high bm
 !          In this case, in GRISLI, we put constant values (0.2)
-      if (bmshelfCLIO(i,j,20).lt.999d0) then! 
+      if (bmshelfCLIO(i,j,20).lt.99d0) then! 
          bmshelf(i,j) = real(bmshelfCLIO(i,j,20))
          bmgrz(i,j) = bmshelf(i,j) * bmgrz_fact
@@ -172 +172 @@
          do noc=20,2,-1
             if ( locdepth.gt. z_CLIO(noc)+dz_CLIO(noc)/2. ) then
-               if (bmshelfCLIO(i,j,noc+1).lt.999d0) then
+               if (bmshelfCLIO(i,j,noc+1).lt.99d0) then
                   bmelt(i,j) = real(bmshelfCLIO(i,j,noc+1))
                endif