Changeset 157 for trunk/SOURCES/steps_time_loop.f90

Timestamp:

11/10/17 17:04:49 (7 years ago)

Author:

dumas

Message:

determin_marais is now called only when isynchro=1 | sortie_hz_multi and hz_output are suppressed | step_time_loop_avec_iterbeta updated as step_time_loop

File:

• trunk/SOURCES/steps_time_loop.f90 (modified) (5 diffs)

trunk/SOURCES/steps_time_loop.f90

@@ -61 +61 @@
      call ablation_bord
      call bilan_eau
-     call bilan_flux_output
-
-!~      if (isynchro.eq.1) then
-!~          call shortoutput
-!~          diff_H = 0.
-!~          Bm_dtt(:,:) = 0.
-!~          bmelt_dtt(:,:) = 0.
-!~          calv_dtt(:,:)=0.
-!~          ablbord_dtt(:,:)=0.
-!~          diff_H_2D(:,:)=0.
-!~          grline_dtt(:,:)=0.
-!~       endif
-     
+     call bilan_flux_output     
      call flottab
-
-     if (itracebug.eq.1)  call tracebug('apres calving')
+     if (isynchro.eq.1.or.nt.eq.1) call determin_marais
+     if (itracebug.eq.1)  call tracebug('apres marais')
 
      ! new calculation of ice surface elevation
@@ -118 +106 @@
   !====================================================================
 
-  ! horizontal plan snapshots
-  !------------------------------------------------------------------
-  !  Premiere partie des sorties horizontales
-
-  !     if (mod(abs(dble(TIME)),dble(DTSORTIE)).lt.dble(dtmin)) then
-  !
-
-
-  if ((mod(abs(TIME),dtt).lt.dtmin).or.(isynchro.eq.1)) then
-     isynchro=1
-
-     !  sorties hz 
-     call testsort_time(real(time))
-     if (iglob_hz.eq.1) then
-        call sortie_hz_multi   ! pour les variables déclarées dans 3D
-        !              call hz_output(real(time))
-     endif
-
-  else
-     iglob_hz=0
-
-  endif
-
-
   !   vertical plan snapshots (profiles)
   !------------------------------------------------------------------
@@ -148 +112 @@
   !     .or.(mod(abs(dble(TIME)),dble(DTPROFILE)).lt.dble(dtmin)) &
   !     .or.(abs(TIME-TEND).lt.dtmin)) then
-
-     !        call sortieprofile()
-
+  !  call sortieprofile()
   !endif
 
@@ -357 +319 @@
         if (itracebug.eq.1)  call tracebug(' Dans spinup_3_bed')
 
-!        if ((nbed.eq.1).and.nt.ne.1.and.isynchro.eq.1) then
-         if ((nbed.eq.1).and.nt.ne.1.and.isynchro.eq.1.and.(mod(abs(TIME),50.).lt.dtmin)) then
+        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
@@ -449 +410 @@
         else if ((icompteur.ne.1).and.(nt.le.1)) then  ! initialisation avec reprise partielle 
                                                        ! ou nulle du vecteur d' etat
-          iter_visco =  10
+          iter_visco = 10
         else
-           iter_visco = 2
+          iter_visco = 2
         end if
  
-!cdc        iter_visco= 10                                 ! warning, test sur l'impact iteration dragging et visco 
-                                                       ! Cat 18 janv 2013
 
            do m=1,iter_visco              
               call diagnoshelf
-              call  mix_SIA_L1  
-              call strain_rate()
+              call mix_SIA_L1  
+              call strain_rate
               if (iter_beta.eq.0) then
                 call dragging