Changeset 1503 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 26, 2015, 8:41:46 PM (9 years ago)

Author:

flefevre

Message:

Nettoyage du solveur chimique ASIS

File:

• trunk/LMDZ.MARS/libf/aeronomars/photochemistry_asis.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/aeronomars/photochemistry_asis.F90

@@ -177 +177 @@
 
 phychemrat = nint(ptimestep/stimestep)
-phychemrat = 1 ! test 2014_004
+phychemrat = 1
 
 ctimestep = ptimestep/real(phychemrat)
@@ -280 +280 @@
    if (dt_corrected /= dt_guess) then  ! the timestep has been modified
 
-!  form the matrix identity + mat*dt_guess
+!  form the matrix identity + mat*dt_corrected
 
       mat(:,:) = mat1(:,:)*dt_corrected
@@ -314 +314 @@
    dt_guess = dt_corrected     ! first-guess timestep for next iteration
 
-!  if (dt_corrected == dt_min) then
-!     print *, "ilev = ", ilev, " sza = ", sza, " iter = ", iter(ilev)
-!  end if
-
-!  do iesp = 1,nesp 
-!     if (c(ilev,iesp) < 0.) then
-!        print *, "sza = ", sza, " species ", iesp, " = ", c(ilev,iesp)/dens(ilev)
-!        print *, "ilev = ", ilev, " iter = ", iter(ilev), " dt_corrected = ", dt_corrected
-!        print *, "e_o = ", e1, " e_ho2 = ", e2, " e_no2 = ", e3
-!     end if
-!  end do
-
    end do ! while (time < ptimestep)
-
-!  if (iter(ilev) >= 80) then
-!     print *, "*** iter = ", iter(ilev), " ilev = ", ilev, " sza = ", sza
-!     print *, " dt_corrected = ", dt_corrected
-!     stop
-!  end if
-
-!  do iesp = 1,nesp 
-!     if (c(ilev,iesp) < 0.) then
-!        print *, "final value is negative!!"
-!        print *, "sza = ", sza, " species ", iesp, " = ", c(ilev,iesp)/dens(ilev)
-!        print *, "ilev = ", ilev, " iter = ", iter(ilev), " dt_corrected = ", dt_corrected
-!        print *, "e_o = ", e1, " e_ho2 = ", e2, " e_no2 = ", e3
-!        stop
-!     end if
-!  end do
-
-!  save cold and dt for next timestep
-
-!  cold3d(ig,ilev,:) = cold(:)
 
 end do ! ilev