Changeset 3799 for palm/trunk/UTIL/chemistry/gasphase_preproc/kpp4palm/bin/kpp4palm.ksh

Timestamp:

Mar 15, 2019 3:36:25 PM (5 years ago)

Author:

forkel

Message:

editing in kpp4palm: add statements for avoiding unused variables, remove $Id

File:

• palm/trunk/UTIL/chemistry/gasphase_preproc/kpp4palm/bin/kpp4palm.ksh (modified) (2 diffs)

palm/trunk/UTIL/chemistry/gasphase_preproc/kpp4palm/bin/kpp4palm.ksh

@@ -40 +40 @@
 # -----------------
 # $Id$ 
+# Removed edit of phot(nphot) for version edit phot(nphot), now done in fortran_file.C (15.03.2019, forkel)
+# Editing with sed to add dummy statements toavoud unused variables
+# 
+# 3780 2019-03-05 11:19:45Z forkel
 # forkel   05. March 2019: cs_mech and set_cm
 # forkel   30. Oktober 2018: Integrating contents of kp4_compress into chem_gasphase_mod.f90
@@ -337 +341 @@
 $BASE/bin/kpp4palm.exe $PREFIX $MODE $VLEN $DE_INDEX $DE_INDEX_FAST
 
-#Prelimanary, substitution has to be moved into kpp4palm.exe
+# Add dummy statements in order to prevent warnings due to unused variables
+#
+sed -i -e '/cfactor =/a !  ' kk_kpp.f90
+sed -i -e '/cfactor =/a    IF ( lu_crow(1) == 1  .OR.  lu_icol(1) == 1  .OR.  lu_irow(1) == 1 )  CONTINUE ' kk_kpp.f90
+sed -i -e '/cfactor =/a ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/cfactor =/a !  ' kk_kpp.f90
+
 if [[ $MODE = "vector" ]]
 then
-  sed -i -e 's/phot(nphot/phot(vl_dim,nphot/g' kk_kpp.f90
-fi
+sed -i -e '/! Computation of equation rates/i ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/! Computation of equation rates/i IF ( f(vl,nfix) > 0.0_dp )  CONTINUE' kk_kpp.f90
+sed -i -e '/! Computation of equation rates/i !  ' kk_kpp.f90
+else
+sed -i -e '/! Computation of equation rates/i ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/! Computation of equation rates/i IF ( f(nfix) > 0.0_dp )  CONTINUE' kk_kpp.f90
+sed -i -e '/! Computation of equation rates/i !  ' kk_kpp.f90
+fi
+
+if [[ $MODE = "vector" ]]
+then
+sed -i -e '/REAL(kind=dp) :: b/a   IF ( f(vl,nfix) > 0.0_dp )  CONTINUE' kk_kpp.f90
+sed -i -e '/REAL(kind=dp) :: b/a ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/REAL(kind=dp):: b/a !' kk_kpp.f90
+else
+sed -i -e '/REAL(kind=dp):: b/a   IF ( f(nfix) > 0.0_dp )  CONTINUE' kk_kpp.f90
+sed -i -e '/REAL(kind=dp):: b/a ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/REAL(kind=dp):: b/a !' kk_kpp.f90
+fi
+
+sed -i -e '/one=1.0_dp/a       IF ( incx == 0 )  CONTINUE' kk_kpp.f90
+sed -i -e '/one=1.0_dp/a ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/one=1.0_dp/a !  ' kk_kpp.f90
+
+sed -i -e '/IF (alpha .eq. zero)RETURN/i !  ' kk_kpp.f90
+sed -i -e '/IF (alpha .eq. zero)RETURN/i ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/IF (alpha .eq. zero)RETURN/i     IF ( incx == 0  .OR.  incy == 0 )  CONTINUE' kk_kpp.f90
+
+sed -i -e '/INTENT(INOUT):: b(n)/a       IF ( pivot(1) == 0 )  CONTINUE' kk_kpp.f90
+sed -i -e '/INTENT(INOUT):: b(n)/a ! Following line is just to avoid compiler message about unused variables' kk_kpp.f90
+sed -i -e '/INTENT(INOUT):: b(n)/a !  ' kk_kpp.f90
 
 if [[ -e $OUTDIR/${OUTFILE}.f90 ]]