Changeset 271 for trunk/SOURCES/Temperature-routines

Timestamp:

07/08/19 17:25:37 (5 years ago)

Author:

aquiquet

Message:

Further computations of mass conservation in double precision for small timesteps

File:

• trunk/SOURCES/Temperature-routines/icetemp_mod.f90 (modified) (5 diffs)

trunk/SOURCES/Temperature-routines/icetemp_mod.f90

@@ -135 +135 @@
     Call Qprod_ice(Iq)
 
-  
+
     ! Dans Le Socle : Les elements de La Matrice Tridiag Sont Invariants Dans L'Espace
     !$OMP PARALLEL
@@ -179 +179 @@
                Advecx(I,J,:),Advecy(I,J,:),Advec(I,J,:),Ct(I,J,:),Iadvec_w(I,J),Iadvec_e(I+1,J),Iadvec_s(I,J),Iadvec_n(I,J+1), &
                                         Ux(i,j,:),Ux(i+1,j,:),Uy(i,j,:),Uy(i,j+1,:),T(i,j,:),T(i-1,j,:),T(i+1,j,:),T(i,j-1,:),T(i,j+1,:), &
-                                        H(i,j),Ts(i,j),Aa,Bb,Cc,Rr,Dx11,Dou,DTT,Dee,Uzr(i,j,:),Dttdx)
+                                        real(H(i,j)),Ts(i,j),Aa,Bb,Cc,Rr,Dx11,Dou,DTT,Dee,Uzr(i,j,:),Dttdx)
 
           ! variables de calcul dans la glace
@@ -188 +188 @@
 !   If (Itracebug.Eq.1) Write(num_tracebug,*)' avant appel  temp_col  i,j', i,j
           Call temp_col(Nz,Nzm,Nn,T3d_new(i,j,:),Ct(i,j,:),&
-               flot(i,j),H(i,j),Tbmer(i,j),Ibase(i,j),T(i,j,:),Tpmp(i,j,:),Aa,Bb,Cc,Rr,Hh,Ncond,Dee,Cm,&
+               flot(i,j),real(H(i,j)),Tbmer(i,j),Ibase(i,j),T(i,j,:),Tpmp(i,j,:),Aa,Bb,Cc,Rr,Hh,Ncond,Dee,Cm,&
                Dzm,Phid(i,j),Ghf(i,j),Ifail1,Tbdot(i,j),DTT,Dou,Ts(i,j))
 
@@ -207 +207 @@
          !$OMP END DO
 
+
          !$OMP DO
     Do J=1,Ny
@@ -254 +255 @@
     End Do
     !$OMP END DO
+
 
     !     Temperature Sur Les Bords :