Changeset 149

Timestamp:

10/25/17 16:53:45 (7 years ago)

Author:

aquiquet

Message:

Grisli-iLoveclim branch: merged to trunk at revision 148 for calving_frange bug correction

Location:

branches/iLoveclim/SOURCES

Files:

• Ant16_files/bmelt-ant-regions-initmip_mod.f90 (modified) (2 diffs)
• Ant40_files/output_anta40_mod-0.4.f90 (modified) (1 diff)
• Makefile.grisli.inc (modified) (1 diff)
• calving_frange.f90 (modified) (5 diffs)
• eaubasale-0.5_mod.f90 (modified) (5 diffs)

branches/iLoveclim/SOURCES/Ant16_files/bmelt-ant-regions-initmip_mod.f90

@@ -46 +46 @@
     ! Cette routine fait l'initialisation pour la fusion basale.
     ! Elle est appelée par inputfile-vec-0.5.f90
-        
+
     real*8, dimension(:,:),   pointer    :: tab               !< tableau 2d real pointer
     character(len=100) :: file_number_shelves ! fichier avec les zones ice-shelves
@@ -125 +125 @@
   write(num_rep_42,428)'!_______________________________________________________________________' 
 
-  
-  file_bmelt_anom  = trim(dirnameinp)//trim(file_bmelt_anom)
-  call Read_Ncdf_var('abmb',file_bmelt_anom,tab)
-  bmelt_anom (:,:) = Tab(:,:)
+  if ( bmelt_time .eq. 1 ) then  
+     file_bmelt_anom  = trim(dirnameinp)//trim(file_bmelt_anom)
+     call Read_Ncdf_var('abmb',file_bmelt_anom,tab)
+     bmelt_anom (:,:) = Tab(:,:)
+  end if
   
     

branches/iLoveclim/SOURCES/Ant40_files/output_anta40_mod-0.4.f90

@@ -82 +82 @@
 
 ! lecture du fichier avec les corrections de surface
-  call Read_Ncdf_var('z',trim(DIRNAMEINP)//'/corrsurf-initMIP-16km.grd',tab)
-  corrsurf(:,:)=tab(:,:)
+  if (geoplace.eq.'ant16km') then
+      call Read_Ncdf_var('z',trim(DIRNAMEINP)//'/corrsurf-initMIP-16km.grd',tab)
+      corrsurf(:,:)=tab(:,:)
+  else
+      corrsurf(:,:)= 1.
+  end if
 
   open(568,file=trim(dirsource)//'/Fichiers-parametres/short-initMIPnc.dat',status='old')

branches/iLoveclim/SOURCES/Makefile.grisli.inc

@@ -247 +247 @@
         bmelt-ant-regions_mod.o \
         bmelt-ant-regions-oce_mod.o \
+        bmelt-ant-regions-initmip_mod.o \
         fake-routines-ant_mod.o \
         dragging_prescr_beta_mod.o beta_iter_vitbil_mod.o \

branches/iLoveclim/SOURCES/calving_frange.f90

@@ -158 +158 @@
 
 !ifint:    if((front(i,j).gt.0).and.(front(i,j).lt.4)) then
-!cdc pb avec front        ifint:    if(front(i,j).lt.4) then
+!cdc pb avec front    ifint:    if((front(i,j).lt.4).or.((front(i-1,j)+front(i+1,j)+front(i,j-1)+front(i,j+1)).lt.16)) then
         ifint:    if((H(i-1,j).lt.2.).or.(H(i+1,j).lt.2.).or.(H(i,j-1).lt.2.).or.(H(i,j+1).lt.2)) then ! si on est au bord avec test sur H
 ! ifint: le point doit avoir au - 1 voisin mais ne pas etre entouré de glace
@@ -471 +471 @@
 
 
-if ((.not.(testmij.or.testpij.or.testimj.or.testipj))  &  ! pas suffisament alimente
+ if ((.not.(testmij.or.testpij.or.testimj.or.testipj))  &  ! pas suffisament alimente
                 .and.(.not.interieur)) then                ! et pas interieur
 
@@ -493 +493 @@
           ! ATTENTION ne pas mettre ice=0 sinon degradation bilan d'eau (bm et bmelt non comptabilises dans ce cas)
                    endif
-                   !H(i,j)=1             
-                !endif
 
                   end if ifint
@@ -507 +505 @@
       
 ! on met en calving les points detectes iceberg :
-       where (iceberg(:,:))
+       where (iceberg(:,:).and.(H(:,:).gt.0.))
              calv(:,:)=-h(:,:)
              ice(:,:)=0
@@ -515 +513 @@
        endwhere   
       
-       calv_dtt(:,:) = calv_dtt(:,:) + calv(:,:) ! somme du calving sur dtt   
-       ! calv_dtt est remis à 0 dans bilan_eau (tous les dtt)
+       calv_dtt(:,:) = calv_dtt(:,:) + calv(:,:) ! somme du calving sur dtt  
+       ! calv_dtt est remis a 0 dans steps_time_loop (tous les dtt)
        
     end subroutine calving

branches/iLoveclim/SOURCES/eaubasale-0.5_mod.f90

@@ -30 +30 @@
     real :: kond0
     REAL :: INFILTR
-    REAL :: DTWAT               !< pas de temps pour l'eau
-    REAL :: compress_w
     REAL :: hmax_till    !< épaisseur de la couche de till
     REAL :: hmax_wat     !< épaisseur de la couche d'eau dans le till
@@ -51 +49 @@
     REAL,dimension(NX,NY) :: bmelt_w          !< fusion (terme source) exprimé en m d'eau
     REAL,dimension(NX,NY) :: vieuxhwater      !< valeur de hwater au debut de l'appel
-    REAL,dimension(NX,NY) :: tetar            !< pour le routage de l'eau dans l'ocean
 
 
@@ -111 +108 @@
 
 
-
-kond(:,:)=kond0
-
-
 ! Conductivite hydraulique : cond passée en m/an ( car le dt est en années)
 kond(:,:)=kond0
@@ -197 +190 @@
      !$OMP DO
      do j=1,ny
-        do i=1,nx
-           tetar(i,j)=(xlong(i,j))*PI/180. ! pourrait etre fait une fois pour toute
-           PGX(I,J)=101*sin(tetar(i,j))*1.e-2                 
-           PGY(I,J)=101*cos(tetar(i,j))*1.e-2                 
+        do i=1,nx              
 
            !   calcul des potentiels
@@ -215 +205 @@
 
      ! sorties debug 17 juillet 2007
-!     debug_3D(:,:,5)=pot_w(:,:)
-!     debug_3D(:,:,6)=pot_f(:,:)
-!     debug_3D(:,:,7)=pot_w(:,:)+rofreshg*hwater(:,:)
-!     debug_3D(:,:,8)=hwater(:,:)
+     debug_3D(:,:,5)=pot_w(:,:)
+     debug_3D(:,:,6)=pot_f(:,:)
+     debug_3D(:,:,7)=pot_w(:,:)+rofreshg*hwater(:,:)
+     debug_3D(:,:,8)=hwater(:,:)
 !$OMP DO
      do j=2,ny