Changeset 137 in lmdz_wrf

Timestamp:

Jul 29, 2014, 11:34:22 AM (10 years ago)

Author:

lfita

Message:

Removing checks for NaNs?, but keeping initialization of variables

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/calltherm.F90 (modified) (8 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/calltherm.F90

@@ -121 +121 @@
 !$OMP THREADPRIVATE(first)
 !********************************************************
-
-! Lluis
-      INTEGER                                            :: llp
-      CHARACTER(LEN=50)                                  :: lvarname, lfname
-      REAL                                               :: largest
-
-      llp = 734
-      lfname = 'physiq'
-      largest = 10.e5
       if (first) then
         itap=0
@@ -143 +134 @@
       Ale = 0.
 
-      PRINT *,'  Lluis SUm(Alp): ', SUM(Alp),' SUM Ale: ',SUM(Ale)
-      IF (SUM(Alp) /= 0.) STOP
-      IF (SUM(Ale) /= 0.) STOP
-
 ! Incrementer le compteur de la physique
      itap   = itap + 1
@@ -153 +140 @@
 !  ===================
 !         print*,'thermiques: WARNING on passe t au lieu de t_seri'
-
-      lfname='beginning calltherm'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
 
 ! On prend comme valeur initiale des thermiques la valeur du pas
@@ -194 +177 @@
          enddo
          if(nbptspb.GT.0) print*,'Number of points with q_seri(i,k)<=0 ',nbptspb   
-      lfname='calltherm before thermcell_main'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_the'
-      CALL check_var3D(lfname, lvarname, d_t_the, klon, klev, largest, .FALSE.)
-      lvarname = 'Alp'
-      CALL check_var(lfname, lvarname, Alp, klon, largest*10.e5, .FALSE.)
-      lvarname = 'Ale'
-      CALL check_var(lfname, lvarname, Ale, klon, largest*10.e5, .FALSE.)
 
          zdt=dtime/REAL(nsplit_thermals)
@@ -304 +278 @@
 ! Ce serait bien de changer, mai en prenant le temps de vérifier que ca
 ! fait bien ce qu'on croit.
-      lfname='calltherm after thermcell_main'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_the'
-      CALL check_var3D(lfname, lvarname, d_t_the, klon, klev, largest, .FALSE.)
-      lvarname = 'Alp'
-      CALL check_var(lfname, lvarname, Alp, klon, largest*10.e5, .FALSE.)
-      lvarname = 'Ale'
-      CALL check_var(lfname, lvarname, Ale, klon, largest*10.e5, .FALSE.)
-
        flag_bidouille_stratocu=iflag_thermals<=12.or.iflag_thermals==14.or.iflag_thermals==16.or.iflag_thermals==18
 
@@ -344 +308 @@
        fm_therm(:,klev+1)=0.
 
-      lfname='calltherm after transforming from derivate to tendency'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_the'
-      CALL check_var3D(lfname, lvarname, d_t_the, klon, klev, largest, .FALSE.)
-
-
 !  accumulation de la tendance
             d_t_ajs(:,:)=d_t_ajs(:,:)+d_t_the(:,:)
@@ -364 +321 @@
             q_seri(:,:) = q_seri(:,:) + d_q_the(:,:)
            if (prt_level.gt.10) write(lunout,*)'Apres apres thermcell_main OK'
-      lfname='calltherm after add tendencies'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_the'
-      CALL check_var3D(lfname, lvarname, d_t_the, klon, klev, largest, .FALSE.)
-      lvarname = 'Alp_bl'
-      CALL check_var(lfname, lvarname, Alp_bl, klon, largest, .FALSE.)
-      lvarname = 'Alp'
-      CALL check_var(lfname, lvarname, Alp, klon, largest*10.e5, .FALSE.)
-      lvarname = 'Ale_bl'
-      CALL check_var(lfname, lvarname, Ale_bl, klon, largest, .FALSE.)
-      lvarname = 'Ale'
-      CALL check_var(lfname, lvarname, Ale, klon, largest*10.e5, .FALSE.)
 
        DO i=1,klon
@@ -386 +330 @@
         if(prt_level.GE.10) print*,'calltherm i Alp_bl Alp Ale_bl Ale',i,Alp_bl(i),Alp(i),Ale_bl(i),Ale(i)
        ENDDO
-      lfname='calltherm after computing Al[p/e]_bl'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_the'
-      CALL check_var3D(lfname, lvarname, d_t_the, klon, klev, largest, .FALSE.)
-      lvarname = 'Alp_bl'
-      CALL check_var(lfname, lvarname, Alp_bl, klon, largest, .FALSE.)
-      lvarname = 'Alp'
-      CALL check_var(lfname, lvarname, Alp, klon, largest*10.e5, .FALSE.)
-      lvarname = 'Ale_bl'
-      CALL check_var(lfname, lvarname, Ale_bl, klon, largest, .FALSE.)
-      lvarname = 'Ale'
-      CALL check_var(lfname, lvarname, Ale, klon, largest*10.e5, .FALSE.)
 
 !IM 060508 marche pas comme cela !!!        enddo ! isplit