Changeset 165

Timestamp:

12/13/17 18:03:43 (7 years ago)

Author:

dumas

Message:

Spinup with iterbeta is now compatible with steps_time_loop.f90

File:

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

trunk/SOURCES/steps_time_loop.f90

@@ -31 +31 @@
 
   if (itracebug.eq.1)  call tracebug('Entree dans step_time_loop ')
-  if (itracebug.eq.1) write(num_tracebug,*) 'nt = ',nt,'  iter_beta = ',iter_beta
+  if (itracebug.eq.1) write(num_tracebug,*) 'nt = ',nt
 
 
@@ -38 +38 @@
 
 
-  spinup_icethick: if ((ispinup.eq.0.or.ispinup.eq.2.or.nt.lt.nt_init) &
-       .and.(iter_beta.eq.0)) then 
+  spinup_icethick: if (ispinup.eq.0.or.ispinup.eq.2.or.nt.lt.nt_init) then 
 
      !-------------------------------------------------------------------
@@ -45 +44 @@
      !----------------
      !   - est fait dans tous les cas pour les premieres boucles temporelles
-     !   - pour ispinup = 0 (run standard) sauf si iter_beta = 0
+     !   - pour ispinup = 0 (run standard)
      !   - pour ispinup = 2 (force avec vitesses de bilan)
      !
@@ -121 +120 @@
   call testsort_time_ncdf(time)
   if (iglob_ncdf .EQ. 1) call sortie_ncdf_cat
-  if (isynchro.eq.1) then
+  if ((isynchro.eq.1).or.(nt.eq.1)) then
     call shortoutput
     diff_H = 0.
@@ -189 +188 @@
   timemax=time
   if (itracebug.eq.1) write(num_tracebug,*) 'entree dans step_thermomeca', &
-       '   nt=',nt,'iter_beta=',iter_beta
+       '   nt=',nt
 
   !--------------------------------------------------------------------
@@ -208 +207 @@
   !-------------------------------------------
   call masque()
-
-  if (iter_beta.gt.0)   call init_dragging
 
 
@@ -321 +318 @@
      ! Passe dans ce bloc que quand on veut calculer les vitesses dynamiques :
      ! - initialisations (nt.lt.nt_init_tm)
-     ! - run standard  (ispinup=0) (y compris iter_beta.ne.0)
+     ! - run standard  (ispinup=0)
      ! - spinup temperature avec vitesses dynamiques  (ispinup = 1)                
      !---------------------------------------------------------------------------------
@@ -331 +328 @@
      call diffusiv()
 
-
-     iterbeta:     if (iter_beta.ne.0) then          ! only for iterations on beta
-
-        if (.not.shelf_vitbil) then    !  bloc si les vitesses shelves viennent de diagnoshelf
-
-           !  Vcol declare dans le module spinup_mod 
-           !  ou dans le dragging si no_spinup
-
-           where (flotmx(:,:))
-              uxbar(:,:)=uxflgz(:,:)
-           elsewhere
-              uxbar(:,:)=Vcol_x(:,:)
-           end where
-
-           where (flotmy(:,:))
-              uybar(:,:)=uyflgz(:,:)
-           elsewhere
-              uybar(:,:)=Vcol_y(:,:)
-           end where
-        else
-           uxbar(:,:)=Vcol_x(:,:)
-           uybar(:,:)=Vcol_y(:,:)
-        end if
-     end if iterbeta
-
-
      ! strain rate calculations
      !===========================
 
      call strain_rate()
-
-
-     if (iter_beta.ne.0) then               ! on n'a utilise les vitesses de bilan que pour strain_rate "Dirichlet ?"
-        uxbar(:,:)=uxflgz(:,:)+uxdef(:,:)
-        uybar(:,:)=uyflgz(:,:)+uydef(:,:)
-     end if
-
-
 
      ! ice shelves and ice streams velocities : SSA
@@ -393 +356 @@
         end if
  
-
+           test_iter_diag=0.
            do m=1,iter_visco              
               call diagnoshelf
               call mix_SIA_L1  
               call strain_rate
-              if (iter_beta.eq.0) then
-                call dragging          
-              endif
               if (test_iter_diag.lt.1.e-2) exit
            end do
+           
+           ! iterations pour inversion du beta en fonction des vitesses de bilan
+           if ((inv_beta.eq.1).and.(time.ge.time_iter).and.(time.le.time_iter_end)) then  ! iterations pendant 20 ans
+              do m=1,nb_iter_vitbil
+                 call beta_iter_vitbil(m)
+                 call diagnoshelf
+                 call mix_SIA_L1
+                 call strain_rate
+              end do
+           end if
+           ! fin des iterations beta_iter_vitbil
+
         endif isync_2