Changeset 116

Timestamp:

06/20/17 14:51:53 (7 years ago)

Author:

dumas

Message:

Ant-40_iterbeta added in Makefile.grisli.inc | update steps_time_loop_avec_iterbeta.f90 to be coherent with steps_time_loop.f90 | OMP DO shifted in dragging_param_beta_mod.f90

Location:

trunk/SOURCES

Files:

• Makefile.grisli.inc (modified) (2 diffs)
• dragging_param_beta_mod.f90 (modified) (1 diff)
• steps_time_loop_avec_iterbeta.f90 (modified) (6 diffs)

trunk/SOURCES/Makefile.grisli.inc

@@ -235 +235 @@
         bmelt-ant-regions_mod.o \
         bmelt-ant-regions-oce_mod.o \
-        fake-routines-ant_mod.o module_choix-antar40-0.4.o \
+        fake-routines-ant_mod.o \
+        dragging_prescr_beta_mod.o beta_iter_vitbil_mod.o \
+        module_choix-antar40-0.4.o \
         track_ant40_mod.o
         
@@ -689 +691 @@
         $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) $(Liste_BLAS)
         
+Ant-40_iterbeta : $(Dim_Ant40) $(mod_dim_communs) \
+        $(toy_recul)  \
+        $(mod_communs) \
+        $(mod_clim_tof) \
+        $(mod_no_tracers) \
+        $(mod_ell) $(Liste_Ant40) \
+        $(diagnoshelf) \
+        $(Liste_Netcdf) \
+        $(routines_communes) steps_time_loop_avec_iterbeta.o \
+        $(routine_elliptiques) \
+        $(Liste_BLAS)
+
+        $(LK) -o  ../bin/Ant-40_iterbeta \
+        $(Dim_Ant40) $(mod_dim_communs) \
+        $(toy_recul)  \
+        $(mod_communs) \
+        $(mod_clim_tof) \
+        $(mod_no_tracers) \
+        $(mod_ell) $(Liste_Ant40) \
+        $(diagnoshelf) \
+        $(Liste_Netcdf) \
+        $(routines_communes) steps_time_loop_avec_iterbeta.o \
+        $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) $(Liste_BLAS)     
+        
 Snowball : $(Dim_Snowball) $(mod_dim_communs) \
         $(toy_recul)  \

trunk/SOURCES/dragging_param_beta_mod.f90

@@ -152 +152 @@
 
 ! calcul de neff (pression effective sur noeuds majeurs)
-!$OMP DO
 if (sum(neffmx(:,:)).le.0.) neffmx(:,:) =1.e8
 if (sum(neffmy(:,:)).le.0.) neffmy(:,:) =1.e8
 
+!$OMP DO
 do j=1,ny-1
    do i=1,nx-1

trunk/SOURCES/steps_time_loop_avec_iterbeta.f90

@@ -56 +56 @@
      if (itracebug.eq.1)  call tracebug(' Dans spinup_icethick : appelle icethick')
 
-     call icethick3()
+     call icethick3
 
      call calving
@@ -140 +140 @@
   !------------------------------------------------------------------
 
-  if ((NT.eq.1) &    !   .or.(NT.eq.2)
-       .or.(mod(abs(dble(TIME)),dble(DTPROFILE)).lt.dble(dtmin)) &
-       .or.(abs(TIME-TEND).lt.dtmin)) then
+  !if ((NT.eq.1) &    !   .or.(NT.eq.2)
+  !     .or.(mod(abs(dble(TIME)),dble(DTPROFILE)).lt.dble(dtmin)) &
+  !     .or.(abs(TIME-TEND).lt.dtmin)) then
 
      !        call sortieprofile()
 
-  endif
+  !endif
 
   !  sorties netcdf Hassine (aout 2010)  (2D and 3D)
@@ -164 +164 @@
 
   ! end of outputs
+  if (isynchro.eq.1) then
+                 diff_H = 0.
+                 Bm_dtt(:,:) = 0.
+                 bmelt_dtt(:,:) = 0.
+           calv_dtt(:,:)=0.
+                 ablbord_dtt(:,:)=0.
+                 diff_H_2D(:,:)=0.
+        endif   
   !===================================================================
 
@@ -233 +241 @@
   !-------------------------------------------
   call masque()
-  call flottab()
+
   if (iter_beta.gt.0)   call init_dragging
 
@@ -345 +353 @@
         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           
            call bedrock()                                            !  bedrock adjustment
         endif
@@ -444 +452 @@
            do m=1,iter_visco  
               call diagnoshelf
-              call  mix_SIA_L1  
-              call strain_rate()
+              call mix_SIA_L1  
+              call strain_rate
            end do
 
 ! iterations pour affiner beta en fonction des vitesses de bilan
  
-           if ((time.ge.time_iter).and.(time.le.time_iter+15.)) then    !iterations seulment pendant 15 ans
+           if ((time.ge.time_iter).and.(time.le.time_iter+20.)) then    !iterations seulment pendant 15 ans
               do m=1,nb_iter_vitbil
                  call beta_iter_vitbil 
                  call diagnoshelf
-                 call  mix_SIA_L1  
-                 call strain_rate()
+                 call mix_SIA_L1  
+                 call strain_rate
               end do
            end if