Changeset 294

Timestamp:

04/02/20 16:53:50 (4 years ago)

Author:

aquiquet

Message:

Calving version used for ISMIP6-Greenland

File:

• trunk/SOURCES/calving_frange_ISMIP_glaciers.f90 (modified) (1 diff)

trunk/SOURCES/calving_frange_ISMIP_glaciers.f90

@@ -170 +170 @@
 
     ! selection du mask_glaciers fonction de time
-    if(time.lt.time_snap(1)) then              ! time avant le forcage
+    ! /!\ the outputs are written after a dt... so, to be consistent we have to add a dt
+    if(time.lt.(time_snap(1)+dt)) then              ! time avant le forcage
       mask_glaciers_time(:,:) = 0 ! pas de mask
-    else if (time.ge.time_snap(nb_snap)) then  ! time apres le forcage
+    else if (time.ge.(time_snap(nb_snap)+dt)) then  ! time apres le forcage
       mask_glaciers_time(:,:) = mask_glaciers(:,:,nb_snap) ! mask du dernier pas de temps
     else ! cas general
       do k = 1 , nb_snap-1
-        if((time.ge.time_snap(k)).and.(time.lt.time_snap(k+1))) then ! entre k et k+1 
+        if((time.ge.(time_snap(k)+dt)).and.(time.lt.(time_snap(k+1)+dt))) then ! entre k et k+1 
 !          print*,'time et k',time,k, time_snap(k),time_snap(k+1)
           mask_glaciers_time(:,:)=mask_glaciers(:,:,k)          ! mask de time k