Changeset 2708 for trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90

Timestamp:

Jun 21, 2022, 11:05:39 AM (2 years ago)

Author:

aslmd

Message:

error messages when setting a traceur.def with wrong names , or with vap/ice tracer missing for the same specie , or if the specie is not in table_tracers_condensable

File:

• trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90

@@ -82 +82 @@
         qevap(1:ngrid,1:nlayer,1:nq)=0.0
         tevap(1:ngrid,1:nlayer)=0.0
-        igcm_generic_vap = -1
-        igcm_generic_ice = -1
         ! Let's loop on tracers 
         do iq=1,nq
@@ -95 +93 @@
 
                         igcm_generic_vap=iq
+                        igcm_generic_ice = -1
 
                         ! We look for the corresponding ice traceur (before or after in the list of traceurs, maybe could be generalised to the whole list)
-                        if ((noms(iq)(9:len(trim(noms(iq)))-4) .eq. noms(iq+1)(9:len(trim(noms(iq+1)))-4)) .and. (index(noms(iq+1),"ice") .ne. 0)) then
-                                igcm_generic_ice = iq+1
-                        else if ((iq .gt. 1) .and. (noms(iq)(9:len(trim(noms(iq)))-4) .eq. noms(iq-1)(9:len(trim(noms(iq-1)))-4)) .and. (index(noms(iq-1),"ice") .ne. 0)) then
-                                igcm_generic_ice = iq-1
-                        else
+                        if (iq .ne. nq) then
+                                if ((noms(iq)(1:len(trim(noms(iq)))-4) .eq. noms(iq+1)(1:len(trim(noms(iq+1)))-4)) .and. (index(noms(iq+1),"ice") .ne. 0)) then
+                                        igcm_generic_ice = iq+1
+                                end if
+                        end if
+                        if ((iq .gt. 1)) then
+                                if ((noms(iq)(1:len(trim(noms(iq)))-4) .eq. noms(iq-1)(1:len(trim(noms(iq-1)))-4)) .and. (index(noms(iq-1),"ice") .ne. 0)) then
+                                        igcm_generic_ice = iq-1
+                                end if
+                        end if
+                        if (igcm_generic_ice .eq. -1) then
                                 write(*,*) "ERROR : You set a vap traceur but you forgot to set the corresponding ice traceur"
                         endif
 
                         ! Need to call specie_parameters of the considered specie 
-                        call specie_parameters(noms(iq)(9:len(trim(noms(iq)))-4))
-                        !call specie_parameters_table(noms(iq)(9:len(trim(noms(iq)))-4))
+                        !call specie_parameters(noms(iq)(9:len(trim(noms(iq)))-4))
+                        write(*,*) 'looking specie parameters for : ',noms(iq)(1:len(trim(noms(iq)))-4)
+                        call specie_parameters_table(noms(iq)(1:len(trim(noms(iq)))-4))
                         
                         Lcp=RLVTT/cpp ! need to be init here