Changeset 148

Timestamp:

10/25/17 15:16:28 (7 years ago)

Author:

aquiquet

Message:

Minor bug correction in flottab + compatibility between ant40 and ant16 geometries + get rid of unused variable in eaubasale.

Location:

trunk/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) (1 diff)
• eaubasale-0.5_mod.f90 (modified) (5 diffs)

trunk/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
   
     

trunk/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')

trunk/SOURCES/Makefile.grisli.inc

@@ -249 +249 @@
         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 \

trunk/SOURCES/calving_frange.f90

@@ -412 +412 @@
       
 ! on met en calving les points detectes iceberg :
-       where (iceberg(:,:))
-                                                 calv(:,:)=-h(:,:)
+       where (iceberg(:,:).and.(H(:,:).gt.0.))
+             calv(:,:)=-h(:,:)
              ice(:,:)=0
              H(:,:)=0.

trunk/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