Remove debug output added while fixing synchdmlc

ftessier · ftessier · commit 91b6574c9d47 · 2022-06-30T17:52:07.000-04:00
diff --git a/HEN_HOUSE/omega/beamnrc/CMs/SYNCHDMLC_cm.mortran b/HEN_HOUSE/omega/beamnrc/CMs/SYNCHDMLC_cm.mortran
@@ -783,15 +783,6 @@ ELSEIF(ORIENT_$SYNCHDMLC=3)[
 ]
 ;
 
-"write(*,*)' in howfar, ir, med ',IR_$SYNCHDMLC,med(IRL);
-
-"IF(Z(NP)>47.62001 & Z(NP)<54.37999)["
-        "IF(med(irl) = 0)["
-         " write(*,*)' iq,x,y,z ',iq(np),x(np),y(np),z(np);"
-         " write(*,*)' ir,med ',ir(np),med(irl);"
-         " write(*,*)'  ihstry ',ihstry;"
-       " ]"
-  " ]"
 STEP_UNIT=0.0;
 COUNT = 0;
 OUTOFCMFLAG=0;
@@ -820,7 +811,6 @@ OUTOFMLCFLAG=0;
 IF(N_GAP_$SYNCHDMLC=1 & IR_$SYNCHDMLC=4 ) ["Check the air gap first"
    IF(W(NP)>0.0) ["Particle going forward"
       DIST = (ZMIN_$SYNCHDMLC - Z(NP))/W(NP); "distance to front of CM"
-      "write(*,*)' dist ',dist;
       IF(DIST <= 0.0) [
           USTEP=0.0;
           $SYNCHDMLC_FIND(IR_$SYNCHDMLC, 0.0);
@@ -852,7 +842,6 @@ ELSEIF(IR_$SYNCHDMLC=1 | IR_$SYNCHDMLC=2 | IR_$SYNCHDMLC=3) [
    ELSEIF(W(NP) < 0.0) [
      "particle going backward"
       DIST = (ZMIN_$SYNCHDMLC - Z(NP))/W(NP); "distance to back of CM "
-      "write(*,*)' dist, ustep ',dist,ustep;
       IF(DIST>USTEP) [ STEP_UNIT=USTEP;  ]
       ELSE [ STEP_UNIT=DIST; ]
       IF(DIST <= 0.0) [
@@ -870,18 +859,13 @@ ELSEIF(IR_$SYNCHDMLC=1 | IR_$SYNCHDMLC=2 | IR_$SYNCHDMLC=3) [
 
 $SYNCHDMLC_FIND(REGION_$SYNCHDMLC,0.0);
 
-"write(*,*)' leafis, region, dist ',leafis,REGION_$SYNCHDMLC,DIST;
-
 IF(REGION_$SYNCHDMLC=4) [ "Particle is in the air gap "
   REGION_$SYNCHDMLC=IRSTART_$SYNCHDMLC-1+REGION_$SYNCHDMLC;
    IF(W(NP) > 0.0)[
     IF(DIST <= USTEP ) ["particle to be moved to region boundary"
-      "write(*,*)' ustep, dist here ',ustep,dist;
       $GEO_SHIFT_1_(DIST);
-      "write(*,*)' and ustep, dist here ',ustep,dist;
       USTEP = DIST;
       $SYNCHDMLC_FIND(NEWREGION_$SYNCHDMLC,USTEP);
-      "write(*,*)' ustep, newregion ',ustep,NEWREGION_$SYNCHDMLC;
       IF( NEWREGION_$SYNCHDMLC=1) [IRNEW =IRSTART_$SYNCHDMLC; RETURN;]
       ELSE[IRNEW =IRSTART_$SYNCHDMLC+1; RETURN;]
     ]
@@ -946,11 +930,6 @@ IF(OUTOFMLCFLAG=1) [
 
 $SYNCHDMLC_MINDISTANCE(DIST1);
 
-"write(*,*)' dist1 ',dist1;
-"write(*,*)' x,y,z ',x(np),y(np),z(np);
-"write(*,*)' leafis,nx,ny,nz ',leafis,nx,ny,nz;
-"write(*,*)' dist1,step_unit ',DIST1,STEP_UNIT;
-
 "OUTPUT REGION_$SYNCHDMLC;(I3);"
 
 IF((DIST1>=STEP_UNIT) | (DIST1<0)) [
@@ -983,9 +962,6 @@ ELSEIF( (STEP_UNIT>DIST1) & (DIST1>=0)) [
         NEWREGION_$SYNCHDMLC=NEWREGION_$SYNCHDMLC+IRSTART_$SYNCHDMLC-1;
        " OUTPUT REGION_$SYNCHDMLC;(I3);"
 
-       "write(*,*)' region, newregion, dist1 ',REGION_$SYNCHDMLC,
-       "     NEWREGION_$SYNCHDMLC,DIST1;
-
         IF( OUTOFCMFLAG=1 ) [
           EXIT;
         ] ELSEIF( (NEWREGION_$SYNCHDMLC~=REGION_$SYNCHDMLC)
@@ -1018,7 +994,6 @@ ELSEIF( (STEP_UNIT>DIST1) & (DIST1>=0)) [
           RETURN;
        ] ELSE [
          $SYNCHDMLC_MINDISTANCE(DIST2);
-         "write(*,*)' dist1, dist2 ',dist1,dist2;
          DIST1=DIST1+DIST2;
        ];
 
@@ -3730,8 +3705,6 @@ IF(LEAFIS~=0) [
   ]
 
   "Check if we are in groove on leaf and move to adjacent leaf"
-"write(*,*)' before check leafis, leaftype ',leafis,LEAFTYPE_$SYNCHDMLC(LEAFIS);
-"write(*,*)' nx, ny, nz ',nx,ny,nz; 
   IF(LEAFTYPE_$SYNCHDMLC(LEAFIS)=2 | LEAFTYPE_$SYNCHDMLC(LEAFIS)=4) [
 				   "target leaf"
     IF(LEAFIS~=TOT_LEAF_$SYNCHDMLC & ((((NY=5|NY=6) & (NZ=7|NZ=8|NZ=9)) &
@@ -4145,9 +4118,6 @@ IF(LEAFIS~=0) [
 
     DIST=MIN(ABS(ZREG_$SYNCHDMLC(LEAFIS,I1)-Z(NP)),
              ABS(ZREG_$SYNCHDMLC(LEAFIS,I2)-Z(NP)));
-    "write(*,*)' x,y,z ',x(np),y(np),z(np);
-    "write(*,*)' nx,ny,nz ',nx,ny,nz;
-    "write(*,*)' leafis,i1,i2,zdist ',leafis,i1,i2,dist;
 
     "Calculates distance to closest boundary in z direction"
     "--> 3) 7 Jul 2009"
@@ -4656,7 +4626,6 @@ IF(LEAFIS~=0) [
           SQRT(1+TEMP2**2));
 
     DIST = MIN(DIST,TEMP1,TEMP2);
-    "write(*,*)' ydist ',dist;
  "NOW CHECK DISTANCE IN X DIRECTION"
     IF(ENDTYPE_$SYNCHDMLC=1) [ "focused leaf end"
       TEMP1=SURPARA2_B_$SYNCHDMLC(LEAFIS);
diff --git a/HEN_HOUSE/omega/beamnrc/CMs/SYNCHDMLC_macros.mortran b/HEN_HOUSE/omega/beamnrc/CMs/SYNCHDMLC_macros.mortran
@@ -518,8 +518,6 @@ ELSE ["inside leaf bank"
       NY=NY+1;
     ];
 
-    "write(*,*)' leafis, xyfl1, ny ',leafis,xyfl(1),ny;"
-
     IF(LEAFTYPE_$SYNCHDMLC(LEAFIS)=2 | LEAFTYPE_$SYNCHDMLC(LEAFIS)=4) [
        "Originally on Target (Quarter or Half) Leaf"
 
@@ -698,9 +696,6 @@ ELSE ["inside leaf bank"
     ELSEIF((XYFL(2)>TEMP4) | (XYFL(2)=TEMP4 & UVL(2)>=0)) [
          NX=5;]
  ];
- "write(*,*)' find: leafis, x, y, z ',leafis,xyfl(2),xyfl(1),zfl;"
- "write(*,*)' nx, ny, nz, subindex ',nx,ny,nz,"
- "            SUBINDEX_$SYNCHDMLC(LEAFIS,NX,NY,NZ);"
  {P1}=SUBINDEX_$SYNCHDMLC(LEAFIS,NX,NY,NZ);
 
  "OUTPUT LEAFIS,XYFL(2),XYFL(1),ZFL,{P1};"
@@ -1401,7 +1396,6 @@ REPLACE {$SYNCHDMLC_MINDISTANCE(#);} WITH {;
   ELSE [ YN = 1.0E20; ];
   IF((YP>=0.0) & (YN>=0.0) ) [ YDIST =MIN(YP, YN); ]
   ELSE [ YDIST = MAX(YP,YN);];
-  "write(*,*)' i,j, rhs, lhs, ydist ',i,j,rhs,lhs,ydist;"
  IF(ENDTYPE_$SYNCHDMLC = 1) ["focused leaf end"
       HOLE=HOLEPOS_$SYNCHDMLC(LEAFIS);
       TEMP3=LEAFB_$SYNCHDMLC(LEAFIS)-HOLE;
@@ -1456,7 +1450,6 @@ REPLACE {$SYNCHDMLC_MINDISTANCE(#);} WITH {;
     TEMP3 = UVL(2)*(XYFL(2)-XoP)+W(NP)*(ZFL-Zo); "B in eqn"
     TEMP4 = (XYFL(2)-XoN)**2 + (ZFL-Zo)**2; "part of C"
     TEMP5 = (XYFL(2)-XoP)**2 + (ZFL-Zo)**2; "part of C"
-    "write(*,*)' nx ',nx;"
     IF((XYFL(2)>=XoN) & (SQRT((XYFL(2)-XoN)**2+(ZFL-Zo)**2)-
       LEAFRADIUS_$SYNCHDMLC=0.0) & (UVL(2)*(XYFL(2)-XoN)+W(NP)*(ZFL-Zo)=0.0)) [
          XN=0.0;  XP=1.0E20; "on edge of negative leaf"
@@ -1572,7 +1565,6 @@ REPLACE {$SYNCHDMLC_MINDISTANCE(#);} WITH {;
  ];
 IF( (XP>=0.0) & (XN>=0.0) ) [ XDIST =MIN(XP, XN); ]
 ELSE [ XDIST = MAX(XP,XN); ];
-"write(*,*)' xdist ',xdist;"
 IF( ZDIST>=0.0) [TEMP=ZDIST;
     IF( YDIST>=0.0 ) [TEMP=MIN(TEMP,YDIST);];
     IF( XDIST>=0.0 ) [TEMP=MIN(TEMP,XDIST);];
diff --git a/HEN_HOUSE/src/egsnrc.mortran b/HEN_HOUSE/src/egsnrc.mortran
@@ -1041,117 +1041,6 @@ $start_new_particle;
 
 :NEWELECTRON:LOOP
 [
-IF(Y(NP) > -0.73059 & Y(NP) < -0.689 &
-    Z(NP) > 50.998 & ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-  write(*,*)' half air test 1';
-  write(*,*)' e,x,y,z,ir ',E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.8258 & Y(NP) < -0.689 &
-    Z(NP) > 47.716 & Z(NP) < 47.979 & ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-  write(*,*)' half air test 2';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.693 & Y(NP) < -0.689 & ABS(X(NP))<10.0 & Z(NP) > 47.716 &
-   MED(IR(NP))=2)[
-  write(*,*)' half air test 3';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9259 & Y(NP) <-0.8234 & Z(NP) >53.998 & ABS(X(NP))<10.0 &
-  MED(IR(NP))=2)[
-  write(*,*)' half air test 4';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) < -0.689 & Z(NP) >54.1161 & ABS(X(NP))<10.0 &
-  MED(IR(NP))=2)[
-  write(*,*)' half air test 5';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) < -0.926 & Z(NP) >50.996 & ABS(X(NP))<10.0 &
-  MED(IR(NP))=2)[
-  write(*,*)' half air test 6';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9259 & Y(NP) <-0.7306 & Z(NP) <53.996 & Z(NP)>47.981 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 7';
- write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9259 & Y(NP) <-0.6933 & Z(NP) <50.996 & Z(NP)>47.981 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 8';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) <-0.689 & Z(NP) <47.714 & Z(NP)>47.6201 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 9';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) <-0.826 & Z(NP) <50.994 & Z(NP)>47.6201 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 10';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.823 & Y(NP) <-0.7306 & Z(NP) <54.116 & Z(NP)>53.998 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 11';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.4566 & Z(NP) <47.884 & Z(NP)>47.621 &
-"   ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 1';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.7266 & Z(NP) <54.284 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 2';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"  stop;
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.6893 & Z(NP) <51.103 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 3';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.6106 & Z(NP) <54.284 & Z(NP)>54.021 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 4';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6142 & Y(NP) <-0.4566 & Z(NP) <48.00 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 5';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.4936 & Y(NP) <-0.4566 & Z(NP) <51.2 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 6';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.4566 & Z(NP) <54.375 & Z(NP)>54.286 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 7';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6104 & Y(NP) <-0.4566 & Z(NP) <54.375 & Z(NP)>51.206 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 8';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7265 & Y(NP) <-0.4938 & Z(NP) <54.02 & Z(NP)>51.105 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 9';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6892 & Y(NP) <-0.4938 & Z(NP) <51.103 & Z(NP)>48.005 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 10';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6892 & Y(NP) <-0.6143 & Z(NP) <48.003 & Z(NP)>47.886 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 11';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
     "Go once through this loop for each 'new' electron whose charge and
     "energy has not been checked
 
@@ -1293,7 +1182,6 @@ IF(Y(NP) > -0.823 & Y(NP) <-0.7306 & Z(NP) <54.116 & Z(NP)>53.998 &
                 $SET-TUSTEP-EM-FIELD; "optional tustep restriction in EM field
 
                 $CALL-HOWNEAR(tperp);
-                "write(*,*)'ustep, tperp ',ustep,tperp;
                 dnear(np) = tperp;
                 $RANGE-DISCARD;       "optional regional range rejection for"
                                       "particles below e_max_rr if i_do_rr set"
@@ -6592,118 +6480,6 @@ $EVALUATE GMFPR0 USING GMFP(GLE);]
 
 :PTRANS:
 LOOP["PHOTON TRANSPORT LOOP"
-IF(Y(NP) > -0.73059 & Y(NP) < -0.689 &
-    Z(NP) > 50.998 & ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-  write(*,*)' half air test 1';
-  write(*,*)' e,x,y,z,ir ',E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.8258 & Y(NP) < -0.689 &
-    Z(NP) > 47.716 & Z(NP) < 47.979 & ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-  write(*,*)' half air test 2';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.693 & Y(NP) < -0.689 & ABS(X(NP))<10.0 & Z(NP) > 47.716 &
-   MED(IR(NP))=2)[
-  write(*,*)' half air test 3';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9259 & Y(NP) <-0.8234 & Z(NP) >53.998 & ABS(X(NP))<10.0 &
-  MED(IR(NP))=2)[
-  write(*,*)' half air test 4';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) < -0.689 & Z(NP) >54.1161 & ABS(X(NP))<10.0 &
-  MED(IR(NP))=2)[
-  write(*,*)' half air test 5';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) < -0.926 & Z(NP) >50.996 & ABS(X(NP))<10.0 &
-  MED(IR(NP))=2)[
-  write(*,*)' half air test 6';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9259 & Y(NP) <-0.7306 & Z(NP) <53.996 & Z(NP)>47.981 &
- ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 7';
- write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9259 & Y(NP) <-0.6933 & Z(NP) <50.996 & Z(NP)>47.981 &
- ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 8';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) <-0.689 & Z(NP) <47.714 & Z(NP)>47.6201 &
- ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 9';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.9632 & Y(NP) <-0.826 & Z(NP) <50.994 & Z(NP)>47.6201 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 10';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-IF(Y(NP) > -0.823 & Y(NP) <-0.7306 & Z(NP) <54.116 & Z(NP)>53.998 &
-  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-  write(*,*)' half med test 11';
-  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.4566 & Z(NP) <47.884 & Z(NP)>47.621 &
-"   ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 1';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.7266 & Z(NP) <54.284 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 2';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.6893 & Z(NP) <51.103 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 3';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"  stop;
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.6106 & Z(NP) <54.284 & Z(NP)>54.021 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 4';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6142 & Y(NP) <-0.4566 & Z(NP) <48.00 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 5';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.4936 & Y(NP) <-0.4566 & Z(NP) <51.2 & Z(NP)>47.621 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=2)[
-"  write(*,*)' iso air test 6';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7274 & Y(NP) <-0.4566 & Z(NP) <54.375 & Z(NP)>54.286 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 7';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6104 & Y(NP) <-0.4566 & Z(NP) <54.375 & Z(NP)>51.206 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 8';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.7265 & Y(NP) <-0.4938 & Z(NP) <54.02 & Z(NP)>51.105 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 9';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6892 & Y(NP) <-0.4938 & Z(NP) <51.103 & Z(NP)>48.005 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 10';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-"IF(Y(NP) > -0.6892 & Y(NP) <-0.6143 & Z(NP) <48.003 & Z(NP)>47.886 &
-"  ABS(X(NP))<10.0 & MED(IR(NP))=1)[
-"  write(*,*)' iso med test 11';
-"  write(*,*)' iq,e,x,y,z,ir ',IQ(NP),E(NP),X(NP),Y(NP),Z(NP),IR(NP);
-"]
-
 IF (MEDIUM.EQ.0)[TSTEP=VACDST;]
 ELSE [$SET-RHOF;    "DENSITY RATIO SCALING TEMPLATE"
 GMFP=GMFPR0/RHOF;
