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

Timestamp:

09/15/16 17:52:05 (8 years ago)

Author:

dumas

Message:

Bug correction in flottab : call to determin_tache suppressed | sealevel added in Netcdf output class 1

File:

• trunk/SOURCES/Ant40_files/lect-anteis_mod.f90 (modified) (3 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 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 :
@@ -161 +161 @@
      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)
+!~ !cdc correction de 3 pts qui posent probleme (pente tres forte):
+!~      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.
@@ -235 +228 @@
     do J=1,NY
        do I=1,NX
-          if ((BSOC(I,J)+H(I,J)*RO/ROW -SEALEVEL).LT.0.) then
+          if (((BSOC(I,J)+H(I,J)*RO/ROW -SEALEVEL).LT.0.).and.(H(I,J).gt.1.E-3)) then
              FLOT(I,J)=.TRUE.
           else