Changeset 85 for trunk/SOURCES/Ant40_files/lect-anteis_mod.f90

Timestamp:

08/12/16 16:00:17 (8 years ago)

Author:

dumas

Message:

Open MP validated on Antarctique 15 km grid | New compilation option for large grid size | Call of subroutine courbure is commentated because of writing outsize the dimension of the array

File:

• trunk/SOURCES/Ant40_files/lect-anteis_mod.f90 (modified) (4 diffs)

trunk/SOURCES/Ant40_files/lect-anteis_mod.f90

@@ -114 +114 @@
      call Read_Ncdf_var('H',topo_ref,tab)
      H0(:,:)  = tab(:,:)
+!cdc correction de 3 pts qui posent probleme (pente tres forte):
+     S0(38,53)=1400.
+     Bsoc0(38,53)=-1000.
+     H0(38,53)=S0(38,53)-Bsoc0(38,53)
+     
+     S0(35,53)=1300.
+     Bsoc0(35,53)=-1000.
+     H0(35,53)=S0(35,53)-Bsoc0(35,53)
+     
+     S0(35,56)=1200.
+     Bsoc0(35,56)=-1000.
+     H0(35,56)=S0(35,56)-Bsoc0(35,56)
+
 !cdc correction point pole sud :
 !    S0(71,71)=(S0(71,70)+S0(71,72)+S0(70,71)+S0(72,71))/4.
@@ -147 +160 @@
      call Read_Ncdf_var('H',topo_dep,tab)
      H(:,:)  = tab(:,:)
+     
+!cdc correction de 3 pts qui posent probleme (pente tres forte):
+!~     do i=34,39
+!~       do j=52,57
+!~         print*,'i j', i, j
+!~         print*,'SHB',S(i,j),Bsoc(i,j),H(i,j)
+!~       enddo
+!~    enddo
+
+     S(39,54)=1400.
+     Bsoc(39,54)=-1000.
+     H(39,54)=S(39,54)-Bsoc(39,54)
+     
+     S(36,54)=1300.
+     Bsoc(36,54)=-1000.
+     H(36,54)=S(36,54)-Bsoc(36,54)
+     
+     S(36,57)=1200.
+     Bsoc(36,57)=-1000.
+     H(36,57)=S(36,57)-Bsoc(36,57)
     
 !    S(71,71)=(S(71,70)+S(71,72)+S(70,71)+S(72,71))/4.
@@ -159 +192 @@
 
 
-
 !!$
 !!$    ! socle
@@ -214 +246 @@
 ! calcul des courbures du socle
 
-     call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &
-          bidon(:,:,4),socle_cry,bidon(:,:,5))
-     socle_cry(:,:)=socle_cry(:,:)*dx*dx
+!     call courbure(nx,ny,dx,Bsoc,bidon(:,:,1),bidon(:,:,2),bidon(:,:,3), &
+!          bidon(:,:,4),socle_cry,bidon(:,:,5))
+!     socle_cry(:,:)=socle_cry(:,:)*dx*dx
 
     ! lecture des coordonnées geographiques