Changeset 197

Timestamp:

06/21/18 18:36:23 (6 years ago)

Author:

aquiquet

Message:

Bug correction in bmelt_clio_coupl for iloveclim branch

File:

• branches/iLoveclim/SOURCES/bmelt_clio_coupl_mod.f90 (modified) (4 diffs)

branches/iLoveclim/SOURCES/bmelt_clio_coupl_mod.f90

@@ -50 +50 @@
      write(42,*)'fusion basale sous les ice shelves : bmelt from CLIO   '
      write(42,*)'-------------------------------------------------------'
-
 
 !the depths of the centers of CLIO vertical layers:
@@ -110 +109 @@
               bmgrz(i,j) = bmshelf(i,j) * bmgrz_fact
            endif
-
 ! Now, we look at the depth of the ice shelves to use the right bm
            !if (flot(i,j).and.(H(i,j).gt.1.)) then
@@ -120 +118 @@
               do noc=20,2,-1
                  if ( locdepth .gt. z_CLIO(noc)+dz_CLIO(noc)/2. ) then
-                    if (bmshelfCLIO(i,j,noc+1).lt.99d0) then
-                       bmelt(i,j) = real(bmshelfCLIO(i,j,noc+1))
+                    if (bmshelfCLIO(i,j,noc-1).lt.99d0) then
+                       bmelt(i,j) = real(bmshelfCLIO(i,j,noc-1))
                     endif
                  else
@@ -172 +170 @@
          do noc=20,2,-1
             if ( locdepth.gt. z_CLIO(noc)+dz_CLIO(noc)/2. ) then
-               if (bmshelfCLIO(i,j,noc+1).lt.99d0) then
-                  bmelt(i,j) = real(bmshelfCLIO(i,j,noc+1))
+               if (bmshelfCLIO(i,j,noc-1).lt.99d0) then
+                  bmelt(i,j) = real(bmshelfCLIO(i,j,noc-1))
                endif
             else