Changeset 161 for trunk/SOURCES/steps_time_loop_avec_iterbeta.f90

Timestamp:

11/21/17 17:00:25 (6 years ago)

Author:

dumas

Message:

Attempt to debug the classical crash in diagnoshelf : determin_tache now called by steps_time_loop | isostasy after diagnoshelf | flgzmx & flgzmy .false. if there is no ice

File:

• trunk/SOURCES/steps_time_loop_avec_iterbeta.f90 (modified) (3 diffs)

trunk/SOURCES/steps_time_loop_avec_iterbeta.f90

@@ -58 +58 @@
      call icethick3
      call flottab
+     call determin_tache
      call calving
      call ablation_bord
@@ -305 +306 @@
      end if                                      ! Dans le shelf flgz = Ubar = Vbil et udef=0
 
-
-     ! isostasy
-     !=================
-
-
-     spinup_3_bed: if (ispinup.eq.0.or.nt.lt.nt_init_tm) then
-        !---------------------------------------------------------------------------------
-        ! spinup_3_bed
-        !----------------
-        ! Ne passe dans ce bloc que quand on veut l'isostasie, donc si variations epaisseur
-        ! run standard  (ispinup=0) et pour les initialisations                  
-        !---------------------------------------------------------------------------------
-
-        if (itracebug.eq.1)  call tracebug(' Dans spinup_3_bed')
-
-        if ((nbed.eq.1).and.nt.ne.1.and.isynchro.eq.1.and.(mod(abs(TIME),50.).lt.dtmin)) then           
-           call bedrock()                                            !  bedrock adjustment
-        endif
-
-        call lake_to_slv                                             ! pour traiter dynamiquement les shelves sur lac
-
-     end if spinup_3_bed
-
      call tracer      ! aurel neem (attention, position a verifier)   
 
@@ -441 +419 @@
 
   end if spinup_4_vitdyn
+  
+  ! isostasy
+  !=================
+
+  spinup_3_bed: if ((ISYNCHRO.eq.1).and.(ispinup.eq.0.or.nt.lt.nt_init_tm)) then
+  !---------------------------------------------------------------------------------
+  ! spinup_3_bed
+  !----------------
+  ! Ne passe dans ce bloc que quand on veut l'isostasie, donc si variations epaisseur
+  ! run standard  (ispinup=0) et pour les initialisations                  
+  !---------------------------------------------------------------------------------
+
+  if (itracebug.eq.1)  call tracebug(' Dans spinup_3_bed')
+
+  if ((nbed.eq.1).and.nt.ne.1.and.isynchro.eq.1.and.(mod(abs(TIME),50.).lt.dtmin)) then
+    call bedrock                 !  bedrock adjustment
+  endif
+
+  call lake_to_slv               ! pour traiter dynamiquement les shelves sur lac
+
+  end if spinup_3_bed
+    
 
 end subroutine step_thermomeca