SPAM !!!!!!!!!!!!!!!! [THREAD UFFICIALE!!!]

[antobov]

[antobov]

Adesso mi guardo il terzo episodio di Fringe. Finalmente una serata tranquilla.....

[swattolo]

Io sono esaustoe mi fanno male le gambe, ogni tanto mi viene pure un crampo alla gamba destra, che macello :hypo:

 

Mannaggia alle belle ragazze e lo zio mio allenatore O.o

[antobov]

Io sono esaustoe mi fanno male le gambe, ogni tanto mi viene pure un crampo alla gamba destra, che macello :hypo:

 

Mannaggia alle belle ragazze e lo zio mio allenatore O.o

 

:cheazz: Hai fatto qualche posizione strana? O forse non hai già più l'eta per fare certe cose......

[swattolo]

:cheazz: Hai fatto qualche posizione strana? O forse non hai già più l'eta per fare certe cose......

 

ahahahhahahahahh magari, ora ti spiego XD

Praticamente ho cambiato scietà di atletica e con quella attuale con me c'è una ragazza molto carina.

Tempo fa ad una corsa vicino al mio allenatore (ora ex) c'era un uomo che non conoscevo, così al mio allenatore (ora ex) faccio: ammazza quella che carina, sopratutto un suo lato

L'uomo accanto risponde: ci credo è mia nipote

 

E d'altronte ora quell'uomo è il mio allenatore

 

oggi dopo aver fatto 3,5km di riscaldamento e 10*200m senza riposo, in pratica ritornavi in dietro correndo più piano e ripartivi, dopo due giri defaticanti che io ho fatto affianco a lei, chissà perchè mi fà:

se ti piace tanto stare in compagnia di mia nipote perchè non ci metti a nche 30min di fondo? (corsa lenta) :muro:

Io: vabbhè, almeno sto in buona compagnia

 

Tutti gli altri a farsi la doccia, alc alduccioe Noi due al massacro

 

Vabbhè me la sono cercata

[antobov]

 

Tutti gli altri a farsi la doccia' date=' alc alduccioe Noi due al massacro

 

Vabbhè me la sono cercata [/quote']

 

 

 

Se ne vale la pena, hai fatto bene un po' di sofferenza in questi casi non fa' male.

[swattolo]

Se ne vale la pena, hai fatto bene un po' di sofferenza in questi casi non fa' male.

 

La conosco da pochissimo, però oltre ad essere molto carina è molto sicevole

 

Apparte la mia ex che si sta comportando in modo strano con me e il mio angioletto (quella che mi piace) che misà è ancora innamorata del suo ex non posso lamentarmi :muro:

[Darkcool]

Leggete tutti la recensione della matrox g200!!!!

http://www.xtremehardware.com/forum/f89/retrorecensione_matrox_g200_8mb_matrox_millenium_g200-11580/

 

(all'interno la sfida g200 vs voodoo2!!!!!!)

[apix_1024]

C*************************************************************************

SUBROUTINE POSDAT(KEY,VOL,U,TE,ED,T,P,VIST,DEN,CP,VISM,CON,

* F,ICLMAP,ICTID,RESOR,VF,FORCB,IRN,PREFM,LEVEL)

C Post-process data

C*************************************************************************

INCLUDE 'comdb.inc'

 

DIMENSION KEY(-NBMAXU:NCTMXU),VOL(NCTMXU),U(3,-NBMAXU:NCMAXU),

* TE(-NBMAXU:NCMAXU),ED(-NBMAXU:NCMAXU),T(-NBMAXU:NCTMXU,1+NSCU),

* P(-NBMAXU:NCMAXU),VIST(-NBMAXU:NCMAXU),DEN(-NBMAXU:NCTMXU),

* CP(-NBMAXU:NCTMXU),VISM(-NBMXVU:NCMXVU),CON(-NBMXCU:NCMXCU),

* F(3,-NBMAXU:NCMAXU),ICLMAP(NCTMXU),ICTID(NCTMXU),

* RESOR(89,-100:100),VF(NCDMXU),FORCB(3,NWLMX),IRN(NWLMX)

DOUBLE PRECISION P

DIMENSION PREFM(4)

INCLUDE 'usrdat.inc'

C-------------------------------------------------------------------------------------------------------------------------

C----------------------------------------POSDAT INIEZIONE-----------------------------------------------------------------------

C---------------------------------------

c Calcola grandezze medie delle droplets sull'intera camera (posdat-tot.dat), su 3 sezioni (posdat-sez*.dat) e,

c su ciascuna delle 3 sezioni, in 3 posizioni radiali,alla distanza R dall'iniettore(posdat-sez*R*.dat).

 

c Definire i valori delle sezioni (sez1, sez2, sez3) in mm, definire i valori delle distanze radiali dall'asse iniettore

c (R1,R2,R3) in mm. Mettere 0 se non si vogliono calcolare le grandezze in queste posizioni.

 

c Per il calcolo della penetrazione, indicare se lo spray è un cone o se si tratta di in hollow-cone.

c Nel caso dell'hollow-cone, si stima il main spray escludendo le droplets il cui diametro è inferiore a 2*D10 (D1AVG).

c Occorre inoltre definire una tolleranza (pene-tol, distanza oltre la quale la posdat non prenderà in considerazione

c le droplets per stimare la penetrazione. E' un vincolo necessario per escludere il pre-spray che nella seconda metà della

c simulazione potrebbe essere molto "breakuppat".

 

 

 

 

 

 

C-----------------------------------------------------------------------

C

C *** NEW POINTERS FOR STAR 3.2 - DROPLET ***

C

C-------------- NOMENCLATURE FOR DROPLET -------------------------------

C IDC(i) .... Cell Number of the Droplet (i)

C IDT ....

C DRD .... diameter of the droplet (m)

C DRM .... droplet mass (kg)

C DENDR .... density of droplet (kg/m3)

C IDN .... number of droplets injected

C IDR .... droplet number

C IDRT .... droplet type

C RNDP .... number of droplets represented by the droplet parcel.

C ICO .... droplet component number, or 0 when droplet component number is

C irrelevant.

 

POINTER(P_IDC, IDC(NDRMAX))

POINTER(P_IDNINJ, IDNINJ(NCDMAX))

POINTER(P_IDT, IDT(NDRMAX))

POINTER(P_DRC, DRC(NDRMAX,3))

POINTER(P_DRD, DRD(NDRMAX))

POINTER(P_DRM, DRM(NDRMAX))

POINTER(P_DRT, DRT(NDRMAX))

POINTER(P_DRV, DRV(NDRMAX,3))

POINTER(P_RNDP,RNDP(NDRMAX))

POINTER(P_SPDR,SPDR(NCDMAX))

c POINTER(P_VF, VF(NCDMAX))

POINTER(P_NDRCOM,NDRCOM(1))

 

COMMON /C33001/ P_IDC

COMMON /C33002/ P_IDT

COMMON /C33003/ P_DRD

COMMON /C33004/ P_DRM

COMMON /C33005/ P_DRT

COMMON /C33006/ P_SPDR

c COMMON /C33007/ P_VF

COMMON /C33008/ P_RNDP

COMMON /C33009/ P_IDNINJ

COMMON /DRP000/ P_DRC

COMMON /DRP001/ P_DRV

COMMON /DROP1A1/P_NDRCOM

COMMON /DROP1A /NDR

 

 

dimension posglo(3),posloc(3),vecglo(3),vecloc(3)

 

 

 

c-----------------------------------------------------------------------

c--------------------DEFINIRE I PARAMETRI GENERALI----------------------

c----------------------------------------------------------------------

INTEGER TYPE_POSITION

c---------------------------------------------------------------------

c Se CONE = 1 ----->solid cone spray

c Se CONE = 0 ----->hollow-cone spray

 

parameter (cone=0)

c-------------------------------------------------------------------

c Sezioni su cui calcolare le grandezze

 

parameter (sez1=0.005)

parameter (sez2=0.010)

parameter (sez3=0.015)

c------------------------------------------------------------------

c distanze radiali dall'iniettore

c Se R* = 0 non calcola nulla

 

parameter (R1=0.00)

parameter (R2=0.00)

parameter (R3=0.00)

c------------------------------------------------------------------

c sistema di riferimento dell'iniettore

 

parameter (iniettore=11)

c------------------------------------------------------------------

c semi-angolo di cono, definisce la tolleranza radiale del singolo getto, per calcolare

c la penetrazione, verranno considerate solo le gocce nell'intervallo -z*tan(teta)z*tan(teta)

 

parameter (teta=6.5)

 

c-----------------------------------------------------------------

c Se CONE = 0 settare la tolleranza per la penetrazione

 

parameter (penetol=0.075)

c-----------------------------------------------------------------------

 

 

C***********************************************************************

C********* CALCULATION OF DROPLETS QUANTITIES * START ******************

C***********************************************************************

C ONLY GO THROUGH THE DROPLETS PART IF REQUEST BY POST PROCESSING

C ... total liquid mass in domain and start some averaging

 

IF (LEVEL.EQ.2) THEN

 

C***********************************************************************

C************** Calcolo delle droplets su tutta la camera*********************

 

C***********Inizializzazione************************************************

Droplet_sez1 = 0.

N_Parcel_sez1 = 0.

SUMDM_sez1 = 0.

SUM_ND1_sez1 = 0.

SUM_ND2_sez1 = 0.

SUM_ND3_sez1 = 0.

SUM_D2_sez1 = 0.

SUM_D3_sez1 = 0.

SUM_MD2_sez1 = 0.

SUM_MD3_sez1 = 0.

SUM_VAVG_sez1 = 0.

 

 

C************Calcolo delle grandezze*****************************************

c icsys = iniettore

DO 96 N=1,NDR

IF ( IDC(N).gt.0 )THEN

 

 

N_Parcel = N_Parcel + 1

Droplet = Droplet + RNDP(N)

SUMDM = SUMDM + DRM(N)

 

DRD22 = DRD(N)**2

DRD33 = DRD(N)*DRD22

 

RNDPDRD = RNDP(N) * DRD(N)

RNDPDRD2 = RNDP(N) * DRD22

RNDPDRD3 = RNDP(N) * DRD33

 

SUM_ND1 = SUM_ND1 + RNDPDRD

SUM_ND2 = SUM_ND2 + RNDPDRD2

SUM_ND3 = SUM_ND3 + RNDPDRD3

 

SUM_D2 = SUM_D2 + DRD22

SUM_D3 = SUM_D3 + DRD33

 

DRMDRD2 = DRM(N)*DRD22

DRMDRD3 = DRM(N)*DRD33

 

SUM_MD2 = SUM_MD2 + DRMDRD2

SUM_MD3 = SUM_MD3 + DRMDRD3

 

RNDPVMAG=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG = SUM_VAVG + RNDPVMAG

 

endif

 

96 continue

 

c-----------------------------------------------------------------------------

C***************Fine calcolo su tutta la camera*********************************

c------------------------------------------------------------------------------

 

 

 

 

C****************************************************************************

C************** Calcolo delle droplets sulla sezione 1 (sez1)*****************

 

C***********Inizializzazione**************************************************

 

Droplet_sez1 = 0.

N_Parcel_sez1 = 0.

SUMDM_sez1 = 0.

SUM_ND1_sez1 = 0.

SUM_ND2_sez1 = 0.

SUM_ND3_sez1 = 0.

SUM_D2_sez1 = 0.

SUM_D3_sez1 = 0.

SUM_MD2_sez1 = 0.

SUM_MD3_sez1 = 0.

SUM_VAVG_sez1 = 0.

 

Droplet_sez1R1 = 0.

N_Parcel_sez1R1 = 0.

SUMDM_sez1R1 = 0.

SUM_ND1_sez1R1 = 0.

SUM_ND2_sez1R1 = 0.

SUM_ND3_sez1R1 = 0.

SUM_D2_sez1R1 = 0.

SUM_D3_sez1R1 = 0.

SUM_MD2_sez1R1 = 0.

SUM_MD3_sez1R1 = 0.

SUM_VAVG_sez1R1 = 0.

 

Droplet_sez1R2 = 0.

N_Parcel_sez1R2 = 0.

SUMDM_sez1R2 = 0.

SUM_ND1_sez1R2 = 0.

SUM_ND2_sez1R2 = 0.

SUM_ND3_sez1R2 = 0.

SUM_D2_sez1R2 = 0.

SUM_D3_sez1R2 = 0.

SUM_MD2_sez1R2 = 0.

SUM_MD3_sez1R2 = 0.

SUM_VAVG_sez1R2 = 0.

 

Droplet_sez1R3 = 0.

N_Parcel_sez1R3 = 0.

SUMDM_sez1R3 = 0.

SUM_ND1_sez1R3 = 0.

SUM_ND2_sez1R3 = 0.

SUM_ND3_sez1R3 = 0.

SUM_D2_sez1R3 = 0.

SUM_D3_sez1R3 = 0.

SUM_MD2_sez1R3 = 0.

SUM_MD3_sez1R3 = 0.

SUM_VAVG_sez1R3 = 0.

 

 

C************Calcolo delle grandezze*****************************************

icsys = iniettore

DO 97 N=1,NDR

IF ( IDC(N).gt.0 )THEN

POSGLO(1)=DRC (N,1)

POSGLO(2)=DRC (N,2)

POSGLO(3)=DRC (N,3)

CALL GLOLOC(POSGLO, POSLOC, ICSYS)

 

if ( posloc(3) .gt. sez1-0.005 ) then

if ( posloc(3) .lt. sez1+0.005 ) then

 

N_Parcel_sez1 = N_Parcel_sez1 + 1

Droplet_sez1 = Droplet_sez1 + RNDP(N)

SUMDM_sez1 = SUMDM_sez1 + DRM(N)

 

DRD22_sez1 = DRD(N)**2

DRD33_sez1 = DRD(N)*DRD22_sez1

 

RNDPDRD_sez1 = RNDP(N) * DRD(N)

RNDPDRD2_sez1 = RNDP(N) * DRD22_sez1

RNDPDRD3_sez1 = RNDP(N) * DRD33_sez1

 

SUM_ND1_sez1 = SUM_ND1_sez1 + RNDPDRD_sez1

SUM_ND2_sez1 = SUM_ND2_sez1 + RNDPDRD2_sez1

SUM_ND3_sez1 = SUM_ND3_sez1 + RNDPDRD3_sez1

 

SUM_D2_sez1 = SUM_D2_sez1 + DRD22_sez1

SUM_D3_sez1 = SUM_D3_sez1 + DRD33_sez1

 

DRMDRD2_sez1 = DRM(N)*DRD22_sez1

DRMDRD3_sez1 = DRM(N)*DRD33_sez1

 

SUM_MD2_sez1 = SUM_MD2_sez1 + DRMDRD2_sez1

SUM_MD3_sez1 = SUM_MD3_sez1 + DRMDRD3_sez1

 

RNDPVMAG_sez1=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez1 = SUM_VAVG_sez1 + RNDPVMAG_sez1

 

c end if

c end if

 

 

 

 

c------------ Calcolo grandezze delle droplets sulla sez1 alla distanza R1 dall'asse iniettore

 

if (R1 .gt. 0) then

if (posloc(1).gt.R1-0.005) then

if (posloc(1).lt.R1+0.005) then

N_Parcel_sez1R1 = N_Parcel_sez1R1 + 1

Droplet_sez1R1 = Droplet_sez1R1 + RNDP(N)

SUMDM_sez1R1 = SUMDM_sez1R1 + DRM(N)

 

DRD22_sez1R1 = DRD(N)**2

DRD33_sez1R1 = DRD(N)*DRD22_sez1R1

 

RNDPDRD_sez1R1 = RNDP(N) * DRD(N)

RNDPDRD2_sez1R1 = RNDP(N) * DRD22_sez1R1

RNDPDRD3_sez1R1 = RNDP(N) * DRD33_sez1R1

 

SUM_ND1_sez1R1 = SUM_ND1_sez1R1 + RNDPDRD_sez1R1

SUM_ND2_sez1R1 = SUM_ND2_sez1R1 + RNDPDRD2_sez1R1

SUM_ND3_sez1R1 = SUM_ND3_sez1R1 + RNDPDRD3_sez1R1

 

SUM_D2_sez1R1 = SUM_D2_sez1R1 + DRD22_sez1R1

SUM_D3_sez1R1 = SUM_D3_sez1R1 + DRD33_sez1R1

 

DRMDRD2_sez1R1 = DRM(N)*DRD22_sez1R1

DRMDRD3_sez1R1 = DRM(N)*DRD33_sez1R1

 

SUM_MD2_sez1R1 = SUM_MD2_sez1R1 + DRMDRD2_sez1R1

SUM_MD3_sez1R1 = SUM_MD3_sez1R1 + DRMDRD3_sez1R1

 

RNDPVMAG_sez1R1=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez1R1 = SUM_VAVG_sez1R1 + RNDPVMAG_sez1R1

 

 

 

endif

 

endif

 

endif

 

c-----------fine sezione 1 distanza R1 dal foro-------------------------------

 

 

 

c------------ Calcolo grandezze delle droplets sulla sez1 alla distanza R2 dall'asse iniettore

 

if (R2 .NE. 0) then

if (posloc(1).gt.R2-0.005) then

if (posloc(1).lt.R2+0.005) then

N_Parcel_sez1R2 = N_Parcel_sez1R2 + 1

Droplet_sez1R2 = Droplet_sez1R2 + RNDP(N)

SUMDM_sez1R2 = SUMDM_sez1R2 + DRM(N)

 

DRD22_sez1R2 = DRD(N)**2

DRD33_sez1R2 = DRD(N)*DRD22_sez1R2

 

RNDPDRD_sez1R2 = RNDP(N) * DRD(N)

RNDPDRD2_sez1R2 = RNDP(N) * DRD22_sez1R2

RNDPDRD3_sez1R2 = RNDP(N) * DRD33_sez1R2

 

SUM_ND1_sez1R2 = SUM_ND1_sez1R2 + RNDPDRD_sez1R2

SUM_ND2_sez1R2 = SUM_ND2_sez1R2 + RNDPDRD2_sez1R2

SUM_ND3_sez1R2 = SUM_ND3_sez1R2 + RNDPDRD3_sez1R2

 

SUM_D2_sez1R2 = SUM_D2_sez1R2 + DRD22_sez1R2

SUM_D3_sez1R2 = SUM_D3_sez1R2 + DRD33_sez1R2

 

DRMDRD2_sez1R2 = DRM(N)*DRD22_sez1R2

DRMDRD3_sez1R2 = DRM(N)*DRD33_sez1R2

 

SUM_MD2_sez1R2 = SUM_MD2_sez1R2 + DRMDRD2_sez1R2

SUM_MD3_sez1R2 = SUM_MD3_sez1R2 + DRMDRD3_sez1R2

 

RNDPVMAG_sez1R2=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez1R2 = SUM_VAVG_sez1R2 + RNDPVMAG_sez1R2

 

 

 

 

endif

 

endif

 

endif

 

c-----------fine sezione 1 distanza R2 dal foro-------------------------------

 

 

 

c------------ Calcolo grandezze delle droplets sulla sez1 alla distanza R3 dall'asse iniettore

 

if (R3 .NE. 0) then

if ( posloc(1) .gt. R3-0.005) then

if (posloc(1) .lt. R3+0.005 ) then

N_Parcel_sez1R3 = N_Parcel_sez1R3 + 1

Droplet_sez1R3 = Droplet_sez1R3 + RNDP(N)

SUMDM_sez1R3 = SUMDM_sez1R3 + DRM(N)

 

DRD22_sez1R3 = DRD(N)**2

DRD33_sez1R3 = DRD(N)*DRD22_sez1R3

 

RNDPDRD_sez1R3 = RNDP(N) * DRD(N)

RNDPDRD2_sez1R3 = RNDP(N) * DRD22_sez1R3

RNDPDRD3_sez1R3 = RNDP(N) * DRD33_sez1R3

 

SUM_ND1_sez1R3 = SUM_ND1_sez1R3 + RNDPDRD_sez1R3

SUM_ND2_sez1R3 = SUM_ND2_sez1R3 + RNDPDRD2_sez1R3

SUM_ND3_sez1R3 = SUM_ND3_sez1R3 + RNDPDRD3_sez1R3

 

SUM_D2_sez1R3 = SUM_D2_sez1R3 + DRD22_sez1R3

SUM_D3_sez1R3 = SUM_D3_sez1R3 + DRD33_sez1R3

 

DRMDRD2_sez1R3 = DRM(N)*DRD22_sez1R3

DRMDRD3_sez1R3 = DRM(N)*DRD33_sez1R3

 

SUM_MD2_sez1R3 = SUM_MD2_sez1R3 + DRMDRD2_sez1R3

SUM_MD3_sez1R3 = SUM_MD3_sez1R3 + DRMDRD3_sez1R3

 

RNDPVMAG_sez1R3=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez1R3 = SUM_VAVG_sez1R3 + RNDPVMAG_sez1R3

 

 

 

endif

 

endif

 

endif

 

c-----------fine sezione 1 distanza R3 dal foro-------------------------------

 

 

 

c else

c goto 97

end if

end if

 

END IF

c 97 CONTINUE

c ENDDO

97 continue

c-----------------------------------------------------------------------------

C***************Fine sezione 1*********************************

c------------------------------------------------------------------------------

 

 

C***********************************************************************

C************** Calcolo delle droplets sulla sezione 2 (sez2)*********************

C**********************************************************************

C.............................Inizializzazione...................................................................................

 

Droplet_sez2 = 0.

N_Parcel_sez2 = 0.

SUMDM_sez2 = 0.

 

SUM_ND1_sez2 = 0.

SUM_ND2_sez2 = 0.

SUM_ND3_sez2 = 0.

SUM_D2_sez2 = 0.

SUM_D3_sez2 = 0.

SUM_MD2_sez2 = 0.

SUM_MD3_sez2 = 0.

SUM_VAVG_sez2 = 0.

 

Droplet_sez2R1 = 0.

N_Parcel_sez2R1 = 0.

SUMDM_sez2R1 = 0.

SUM_ND1_sez2R1 = 0.

SUM_ND2_sez2R1 = 0.

SUM_ND3_sez2R1 = 0.

SUM_D2_sez2R1 = 0.

SUM_D3_sez2R1 = 0.

SUM_MD2_sez2R1 = 0.

SUM_MD3_sez2R1 = 0.

SUM_VAVG_sez2R1 = 0.

 

Droplet_sez2R2 = 0.

N_Parcel_sez2R2 = 0.

SUMDM_sez2R2 = 0.

SUM_ND1_sez2R2 = 0.

SUM_ND2_sez2R2 = 0.

SUM_ND3_sez2R2 = 0.

SUM_D2_sez2R2 = 0.

SUM_D3_sez2R2 = 0.

SUM_MD2_sez2R2 = 0.

SUM_MD3_sez2R2 = 0.

SUM_VAVG_sez2R2 = 0.

 

Droplet_sez2R3 = 0.

N_Parcel_sez2R3 = 0.

SUMDM_sez2R3 = 0.

SUM_ND1_sez2R3 = 0.

SUM_ND2_sez2R3 = 0.

SUM_ND3_sez2R3 = 0.

SUM_D2_sez2R3 = 0.

SUM_D3_sez2R3 = 0.

SUM_MD2_sez2R3 = 0.

SUM_MD3_sez2R3 = 0.

SUM_VAVG_sez2R3 = 0.

 

C************Calcolo delle grandezze*****************************************

 

DO 98 N=1,NDR

IF( IDC(N).gt.0 )THEN

POSGLO(1)=DRC (N,1)

POSGLO(2)=DRC (N,2)

POSGLO(3)=DRC (N,3)

CALL GLOLOC(POSGLO, POSLOC, ICSYS)

 

if ( DRC(N,3) .gt. sez2-0.005 ) then

if ( DRC(N,3) .lt. sez2+0.005 ) then

 

N_Parcel_sez2 = N_Parcel_sez2 + 1

Droplet_sez2 = Droplet_sez2 + RNDP(N)

SUMDM_sez2 = SUMDM_sez2 + DRM(N)

 

DRD22_sez2 = DRD(N)**2

DRD33_sez2 = DRD(N)*DRD22_sez2

 

RNDPDRD_sez2 = RNDP(N) * DRD(N)

RNDPDRD2_sez2 = RNDP(N) * DRD22_sez2

RNDPDRD3_sez2 = RNDP(N) * DRD33_sez2

 

SUM_ND1_sez2 = SUM_ND1_sez2 + RNDPDRD_sez2

SUM_ND2_sez2 = SUM_ND2_sez2 + RNDPDRD2_sez2

SUM_ND3_sez2 = SUM_ND3_sez2 + RNDPDRD3_sez2

 

SUM_D2_sez2 = SUM_D2_sez2 + DRD22_sez2

SUM_D3_sez2 = SUM_D3_sez2 + DRD33_sez2

 

DRMDRD2_sez2 = DRM(N)*DRD22_sez2

DRMDRD3_sez2 = DRM(N)*DRD33_sez2

 

SUM_MD2_sez2 = SUM_MD2_sez2 + DRMDRD2_sez2

SUM_MD3_sez2 = SUM_MD3_sez2 + DRMDRD3_sez2

 

RNDPVMAG_sez2=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez2 = SUM_VAVG_sez2 + RNDPVMAG_sez2

 

 

c end if

 

c end if

c------------ Calcolo grandezze delle droplets sulla sez2 alla distanza R1 dall'asse iniettore

 

if (R1 .NE. 0) then

if ( posloc(1) .gt. R1-0.005) then

if ( posloc(1) .lt. R1+0.005 ) then

N_Parcel_sez2R1 = N_Parcel_sez2R1 + 1

Droplet_sez2R1 = Droplet_sez2R1 + RNDP(N)

SUMDM_sez2R1 = SUMDM_sez2R1 + DRM(N)

 

DRD22_sez2R1 = DRD(N)**2

DRD33_sez2R1 = DRD(N)*DRD22_sez2R1

 

RNDPDRD_sez2R1 = RNDP(N) * DRD(N)

RNDPDRD2_sez2R1 = RNDP(N) * DRD22_sez2R1

RNDPDRD3_sez2R1 = RNDP(N) * DRD33_sez2R1

 

SUM_ND1_sez2R1 = SUM_ND1_sez2R1 + RNDPDRD_sez2R1

SUM_ND2_sez2R1 = SUM_ND2_sez2R1 + RNDPDRD2_sez2R1

SUM_ND3_sez2R1 = SUM_ND3_sez2R1 + RNDPDRD3_sez2R1

 

SUM_D2_sez2R1 = SUM_D2_sez2R1 + DRD22_sez2R1

SUM_D3_sez2R1 = SUM_D3_sez2R1 + DRD33_sez2R1

 

DRMDRD2_sez2R1 = DRM(N)*DRD22_sez2R1

DRMDRD3_sez2R1 = DRM(N)*DRD33_sez2R1

 

SUM_MD2_sez2R1 = SUM_MD2_sez2R1 + DRMDRD2_sez2R1

SUM_MD3_sez2R1 = SUM_MD3_sez2R1 + DRMDRD3_sez2R1

 

RNDPVMAG_sez2R1=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez2R1 = SUM_VAVG_sez2R1 + RNDPVMAG_sez2R1

 

 

 

endif

 

endif

 

endif

 

c-----------fine sezione 2 distanza R1 dal foro-------------------------------

 

 

c------------ Calcolo grandezze delle droplets sulla sez2 alla distanza R2 dall'asse iniettore

 

if (R2 .NE. 0) then

if ( posloc(1) .gt. R2-0.005) then

if ( posloc(1) .lt. R2+0.005 ) then

N_Parcel_sez2R2 = N_Parcel_sez2R2 + 1

Droplet_sez2R2 = Droplet_sez2R2 + RNDP(N)

SUMDM_sez2R2 = SUMDM_sez2R2 + DRM(N)

 

DRD22_sez2R2 = DRD(N)**2

DRD33_sez2R2 = DRD(N)*DRD22_sez2R2

 

RNDPDRD_sez2R2 = RNDP(N) * DRD(N)

RNDPDRD2_sez2R2 = RNDP(N) * DRD22_sez2R2

RNDPDRD3_sez2R2 = RNDP(N) * DRD33_sez2R2

 

SUM_ND1_sez2R2 = SUM_ND1_sez2R2 + RNDPDRD_sez2R2

SUM_ND2_sez2R2 = SUM_ND2_sez2R2 + RNDPDRD2_sez2R2

SUM_ND3_sez2R2 = SUM_ND3_sez2R2 + RNDPDRD3_sez2R2

 

SUM_D2_sez2R2 = SUM_D2_sez2R2 + DRD22_sez2R2

SUM_D3_sez2R2 = SUM_D3_sez2R2 + DRD33_sez2R2

 

DRMDRD2_sez2R2 = DRM(N)*DRD22_sez2R2

DRMDRD3_sez2R2 = DRM(N)*DRD33_sez2R2

 

SUM_MD2_sez2R2 = SUM_MD2_sez2R2 + DRMDRD2_sez2R2

SUM_MD3_sez2R2 = SUM_MD3_sez2R2 + DRMDRD3_sez2R2

 

RNDPVMAG_sez2R2=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez2R2 = SUM_VAVG_sez2R2 + RNDPVMAG_sez2R2

 

 

 

endif

 

endif

 

endif

 

c-----------fine sezione 2 distanza R2 dal foro-------------------------------

 

 

 

c------------ Calcolo grandezze delle droplets uslla sez2 alla distanza R3 dall'asse iniettore

 

if (R3 .NE. 0) then

if ( posloc(1) .gt. R3-0.005) then

if ( posloc(1) .lt. R3+0.005 ) then

N_Parcel_sez2R3 = N_Parcel_sez2R3 + 1

Droplet_sez2R3 = Droplet_sez2R3 + RNDP(N)

SUMDM_sez2R3 = SUMDM_sez2R3 + DRM(N)

 

DRD22_sez2R3 = DRD(N)**2

DRD33_sez2R3 = DRD(N)*DRD22_sez2R3

 

RNDPDRD_sez2R3 = RNDP(N) * DRD(N)

RNDPDRD2_sez2R3 = RNDP(N) * DRD22_sez2R3

RNDPDRD3_sez2R3 = RNDP(N) * DRD33_sez2R3

 

SUM_ND1_sez2R3 = SUM_ND1_sez2R3 + RNDPDRD_sez2R3

SUM_ND2_sez2R3 = SUM_ND2_sez2R3 + RNDPDRD2_sez2R3

SUM_ND3_sez2R3 = SUM_ND3_sez2R3 + RNDPDRD3_sez2R3

 

SUM_D2_sez2R3 = SUM_D2_sez2R3 + DRD22_sez2R3

SUM_D3_sez2R3 = SUM_D3_sez2R3 + DRD33_sez2R3

 

DRMDRD2_sez2R3 = DRM(N)*DRD22_sez2R3

DRMDRD3_sez2R3 = DRM(N)*DRD33_sez2R3

 

SUM_MD2_sez2R3 = SUM_MD2_sez2R3 + DRMDRD2_sez2R3

SUM_MD3_sez2R3 = SUM_MD3_sez2R3 + DRMDRD3_sez2R3

 

RNDPVMAG_sez2R3=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez2R3 = SUM_VAVG_sez2R3 + RNDPVMAG_sez2R3

 

 

 

endif

 

endif

 

endif

 

 

c-----------fine sezione 3 distanza R3 dal foro-------------------------------

 

c else

c goto 98

end if

end if

END IF

98 CONTINUE

c END DO

c--------------------------------------------------------------------------------

C***************Fine sezione 2*********************************

c-------------------------------------------------------------------------------

 

 

 

C***********************************************************************

C************** Calcolo delle droplets sulla sezione 3 (sez3)*********************

C**********************************************************************

C.............................Inizializzazione...................................................................................

 

Droplet_sez3 = 0.

N_Parcel_sez3 = 0.

SUMDM_sez3 = 0.

 

SUM_ND1_sez3 = 0.

SUM_ND2_sez3 = 0.

SUM_ND3_sez3 = 0.

SUM_D2_sez3 = 0.

SUM_D3_sez3 = 0.

SUM_MD2_sez3 = 0.

SUM_MD3_sez3 = 0.

SUM_VAVG_sez3 = 0.

 

Droplet_sez3R1 = 0.

N_Parcel_sez3R1 = 0.

SUMDM_sez3R1 = 0.

SUM_ND1_sez3R1 = 0.

SUM_ND2_sez3R1 = 0.

SUM_ND3_sez3R1 = 0.

SUM_D2_sez3R1 = 0.

SUM_D3_sez3R1 = 0.

SUM_MD2_sez3R1 = 0.

SUM_MD3_sez3R1 = 0.

SUM_VAVG_sez3R1 = 0.

 

Droplet_sez3R2 = 0.

N_Parcel_sez3R2 = 0.

SUMDM_sez3R2 = 0.

SUM_ND1_sez3R2 = 0.

SUM_ND2_sez3R2 = 0.

SUM_ND3_sez3R2 = 0.

SUM_D2_sez3R2 = 0.

SUM_D3_sez3R2 = 0.

SUM_MD2_sez3R2 = 0.

SUM_MD3_sez3R2 = 0.

SUM_VAVG_sez3R2 = 0.

 

Droplet_sez3R3 = 0.

N_Parcel_sez3R3 = 0.

SUMDM_sez3R3 = 0.

SUM_ND1_sez3R3 = 0.

SUM_ND2_sez3R3 = 0.

SUM_ND3_sez3R3 = 0.

SUM_D2_sez3R3 = 0.

SUM_D3_sez3R3 = 0.

SUM_MD2_sez3R3 = 0.

SUM_MD3_sez3R3 = 0.

SUM_VAVG_sez3R3 = 0.

 

 

C************Calcolo delle grandezze sulla sezione sez3**********************************

 

DO 99 N=1,NDR

IF( IDC(N).gt.0 )THEN

POSGLO(1)=DRC (N,1)

POSGLO(2)=DRC (N,2)

POSGLO(3)=DRC (N,3)

CALL GLOLOC(POSGLO, POSLOC, ICSYS)

 

if ( DRC(N,3) .gt. sez3-0.005 ) then

if ( DRC(N,3) .lt. sez3+0.005 ) then

 

N_Parcel_sez3 = N_Parcel_sez3 + 1

Droplet_sez3 = Droplet_sez3 + RNDP(N)

SUMDM_sez3 = SUMDM_sez3 + DRM(N)

 

DRD22_sez3 = DRD(N)**2

DRD33_sez3 = DRD(N)*DRD22_sez3

 

RNDPDRD_sez3 = RNDP(N) * DRD(N)

RNDPDRD2_sez3 = RNDP(N) * DRD22_sez3

RNDPDRD3_sez3 = RNDP(N) * DRD33_sez3

 

SUM_ND1_sez3 = SUM_ND1_sez3 + RNDPDRD_sez3

SUM_ND2_sez3 = SUM_ND2_sez3 + RNDPDRD2_sez3

SUM_ND3_sez3 = SUM_ND3_sez3 + RNDPDRD3_sez3

 

SUM_D2_sez3 = SUM_D2_sez3 + DRD22_sez3

SUM_D3_sez3 = SUM_D3_sez3 + DRD33_sez3

 

DRMDRD2_sez3 = DRM(N)*DRD22_sez3

DRMDRD3_sez3 = DRM(N)*DRD33_sez3

 

SUM_MD2_sez3 = SUM_MD2_sez3 + DRMDRD2_sez3

SUM_MD3_sez3 = SUM_MD3_sez3 + DRMDRD3_sez3

 

RNDPVMAG_sez3=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez3 = SUM_VAVG_sez3 + RNDPVMAG_sez3

 

c------Calcolo grandezze delle droplets sulla sezione sez3 alla distanza R1 dall'asse iniettore

 

if (R1 .NE. 0) then

if ( posloc(1) .gt. R1-0.005) then

if ( posloc(1) .lt. R1+0.005 ) then

N_Parcel_sez3R1 = N_Parcel_sez3R1 + 1

Droplet_sez3R1 = Droplet_sez3R1 + RNDP(N)

SUMDM_sez3R1 = SUMDM_sez3R1 + DRM(N)

 

DRD22_sez3R1 = DRD(N)**2

DRD33_sez3R1 = DRD(N)*DRD22_sez3R1

 

RNDPDRD_sez3R1 = RNDP(N) * DRD(N)

RNDPDRD2_sez3R1 = RNDP(N) * DRD22_sez3R1

RNDPDRD3_sez3R1 = RNDP(N) * DRD33_sez3R1

 

SUM_ND1_sez3R1 = SUM_ND1_sez3R1 + RNDPDRD_sez3R1

SUM_ND2_sez3R1 = SUM_ND2_sez3R1 + RNDPDRD2_sez3R1

SUM_ND3_sez3R1 = SUM_ND3_sez3R1 + RNDPDRD3_sez3R1

 

SUM_D2_sez3R1 = SUM_D2_sez3R1 + DRD22_sez3R1

SUM_D3_sez3R1 = SUM_D3_sez3R1 + DRD33_sez3R1

 

DRMDRD2_sez3R1 = DRM(N)*DRD22_sez3R1

DRMDRD3_sez3R1 = DRM(N)*DRD33_sez3R1

 

SUM_MD2_sez3R1 = SUM_MD2_sez3R1 + DRMDRD2_sez3R1

SUM_MD3_sez3R1 = SUM_MD3_sez3R1 + DRMDRD3_sez3R1

 

RNDPVMAG_sez3R1=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez3R1 = SUM_VAVG_sez3R1 + RNDPVMAG_sez3R1

 

 

endif

endif

endif

 

c-----------fine sezione 2 distanza R1 dal foro-------------------------------

 

 

 

c------------ Calcolo grandezze delle droplets sulla sez3 alla distanza R2 dall'asse iniettore

 

if (R2 .NE. 0) then

if ( posloc(1) .gt. R2-0.005) then

if ( posloc(1) .lt. R2+0.005 ) then

N_Parcel_sez3R2 = N_Parcel_sez3R2 + 1

Droplet_sez3R2 = Droplet_sez3R2 + RNDP(N)

SUMDM_sez3R2 = SUMDM_sez3R2 + DRM(N)

 

DRD22_sez3R2 = DRD(N)**2

DRD33_sez3R2 = DRD(N)*DRD22_sez3R2

 

RNDPDRD_sez3R2 = RNDP(N) * DRD(N)

RNDPDRD2_sez3R2 = RNDP(N) * DRD22_sez3R2

RNDPDRD3_sez3R2 = RNDP(N) * DRD33_sez3R2

 

SUM_ND1_sez3R2 = SUM_ND1_sez3R2 + RNDPDRD_sez3R2

SUM_ND2_sez3R2 = SUM_ND2_sez3R2 + RNDPDRD2_sez3R2

SUM_ND3_sez3R2 = SUM_ND3_sez3R2 + RNDPDRD3_sez3R2

 

SUM_D2_sez3R2 = SUM_D2_sez3R2 + DRD22_sez3R2

SUM_D3_sez3R2 = SUM_D3_sez3R2 + DRD33_sez3R2

 

DRMDRD2_sez3R2 = DRM(N)*DRD22_sez3R2

DRMDRD3_sez3R2 = DRM(N)*DRD33_sez3R2

 

SUM_MD2_sez3R2 = SUM_MD2_sez3R2 + DRMDRD2_sez3R2

SUM_MD3_sez3R2 = SUM_MD3_sez3R2 + DRMDRD3_sez3R2

 

RNDPVMAG_sez3R2=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez3R2 = SUM_VAVG_sez3R2 + RNDPVMAG_sez3R2

 

 

 

endif

 

endif

 

endif

 

 

c-----------fine sezione 1 distanza R2 dal foro-------------------------------

 

c------------ Calcolo grandezze delle droplets sulla sez3 alla distanza R3 dall'asse iniettore

 

if (R3 .NE. 0) then

if ( posloc(1) .gt. R3-0.005) then

if (posloc(1) .lt. R3+0.005 ) then

N_Parcel_sez3R3 = N_Parcel_sez3R3 + 1

Droplet_sez3R3 = Droplet_sez3R3 + RNDP(N)

SUMDM_sez3R3 = SUMDM_sez3R3 + DRM(N)

 

DRD22_sez3R3 = DRD(N)**2

DRD33_sez3R3 = DRD(N)*DRD22_sez3R3

 

RNDPDRD_sez3R3 = RNDP(N) * DRD(N)

RNDPDRD2_sez3R3 = RNDP(N) * DRD22_sez3R3

RNDPDRD3_sez3R3 = RNDP(N) * DRD33_sez3R3

 

SUM_ND1_sez3R3 = SUM_ND1_sez3R3 + RNDPDRD_sez3R3

SUM_ND2_sez3R3 = SUM_ND2_sez3R3 + RNDPDRD2_sez3R3

SUM_ND3_sez3R3 = SUM_ND3_sez3R3 + RNDPDRD3_sez3R3

 

SUM_D2_sez3R3 = SUM_D2_sez3R3 + DRD22_sez3R3

SUM_D3_sez3R3 = SUM_D3_sez3R3 + DRD33_sez3R3

 

DRMDRD2_sez3R3 = DRM(N)*DRD22_sez3R3

DRMDRD3_sez3R3 = DRM(N)*DRD33_sez3R3

 

SUM_MD2_sez3R3 = SUM_MD2_sez3R3 + DRMDRD2_sez3R3

SUM_MD3_sez3R3 = SUM_MD3_sez3R3 + DRMDRD3_sez3R3

 

RNDPVMAG_sez3R3=RNDP(N)*SQRT(DRV(N,1)**2+

& DRV(N,2)**2+

& DRV(N,3)**2)

SUM_VAVG_sez3R3 = SUM_VAVG_sez3R3 + RNDPVMAG_sez3R3

 

 

 

endif

 

endif

 

endif

 

 

 

c-----------fine sezione 3 distanza R3 dal foro-------------------------------

 

 

 

 

 

end if

end if

END IF

99 CONTINUE

c END DO

 

C-------------------------------------------------------------------------------

C***************Fine sezione 3*********************************

C--------------------------------------------------------------------------------

 

 

 

 

C..............................Calcolo parallelo.......................................................

 

IF (NHPC.GT.1) THEN

 

SUMDM = GSUM(SUMDM)

N_Parcel = IGSUM(N_Parcel)

Droplet = GSUM(Droplet)

SUM_ND1 = GSUM(SUM_ND1)

SUM_ND2 = GSUM(SUM_ND2)

SUM_ND3 = GSUM(SUM_ND3)

SUM_D2 = GSUM(SUM_D2)

SUM_D3 = GSUM(SUM_D3)

SUM_MD2 = GSUM(SUM_MD2)

SUM_MD3 = GSUM(SUM_MD3)

SUM_VAVG = GSUM(SUM_VAVG)

 

 

 

SUMDM_sez1 = GSUM(SUMDM_sez1)

N_Parcel_sez1 = IGSUM(N_Parcel_sez1)

Droplet_sez1 = GSUM(Droplet_sez1)

SUM_ND1_sez1 = GSUM(SUM_ND1_sez1)

SUM_ND2_sez1 = GSUM(SUM_ND2_sez1)

SUM_ND3_sez1 = GSUM(SUM_ND3_sez1)

SUM_D2_sez1 = GSUM(SUM_D2_sez1)

SUM_D3_sez1 = GSUM(SUM_D3_sez1)

SUM_MD2_sez1 = GSUM(SUM_MD2_sez1)

SUM_MD3_sez1 = GSUM(SUM_MD3_sez1)

SUM_VAVG_sez1 = GSUM(SUM_VAVG_sez1)

 

SUMDM_sez2 = GSUM(SUMDM_sez2)

N_Parcel_sez2 = IGSUM(N_Parcel_sez2)

Droplet_sez2 = GSUM(Droplet_sez2)

SUM_ND1_sez2 = GSUM(SUM_ND1_sez2)

SUM_ND2_sez2 = GSUM(SUM_ND2_sez2)

SUM_ND3_sez2 = GSUM(SUM_ND3_sez2)

SUM_D2_sez2 = GSUM(SUM_D2_sez2)

SUM_D3_sez2 = GSUM(SUM_D3_sez2)

SUM_MD2_sez2 = GSUM(SUM_MD2_sez2)

SUM_MD3_sez2 = GSUM(SUM_MD3_sez2)

SUM_VAVG_sez2 = GSUM(SUM_VAVG_sez2)

 

SUMDM_sez3 = GSUM(SUMDM_sez3)

N_Parcel_sez3 = IGSUM(N_Parcel_sez3)

Droplet_sez3 = GSUM(Droplet_sez3)

SUM_ND1_sez3 = GSUM(SUM_ND1_sez3)

SUM_ND2_sez3 = GSUM(SUM_ND2_sez3)

SUM_ND3_sez3 = GSUM(SUM_ND3_sez3)

SUM_D2_sez3 = GSUM(SUM_D2_sez3)

SUM_D3_sez3 = GSUM(SUM_D3_sez3)

SUM_MD2_sez3 = GSUM(SUM_MD2_sez3)

SUM_MD3_sez3 = GSUM(SUM_MD3_sez3)

SUM_VAVG_sez3 = GSUM(SUM_VAVG_sez3)

 

 

 

 

SUMDM_sez1R1 = GSUM(SUMDM_sez1R1)

N_Parcel_sez1R1 = IGSUM(N_Parcel_sez1R1)

Droplet_sez1R1 = GSUM(Droplet_sez1R1)

SUM_ND1_sez1R1 = GSUM(SUM_ND1_sez1R1)

SUM_ND2_sez1R1 = GSUM(SUM_ND2_sez1R1)

SUM_ND3_sez1R1 = GSUM(SUM_ND3_sez1R1)

SUM_D2_sez1R1 = GSUM(SUM_D2_sez1R1)

SUM_D3_sez1R1 = GSUM(SUM_D3_sez1R1)

SUM_MD2_sez1R1 = GSUM(SUM_MD2_sez1R1)

SUM_MD3_sez1R1 = GSUM(SUM_MD3_sez1R1)

SUM_VAVG_sez1R1 = GSUM(SUM_VAVG_sez1R1)

 

SUMDM_sez1R2 = GSUM(SUMDM_sez1R2)

N_Parcel_sez1R2 = IGSUM(N_Parcel_sez1R2)

Droplet_sez1R2 = GSUM(Droplet_sez1R2)

SUM_ND1_sez1R2 = GSUM(SUM_ND1_sez1R2)

SUM_ND2_sez1R2 = GSUM(SUM_ND2_sez1R2)

SUM_ND3_sez1R2 = GSUM(SUM_ND3_sez1R2)

SUM_D2_sez1R2 = GSUM(SUM_D2_sez1R2)

SUM_D3_sez1R2 = GSUM(SUM_D3_sez1R2)

SUM_MD2_sez1R2 = GSUM(SUM_MD2_sez1R2)

SUM_MD3_sez1R2 = GSUM(SUM_MD3_sez1R2)

SUM_VAVG_sez1R2 = GSUM(SUM_VAVG_sez1R2)

 

SUMDM_sez1R3 = GSUM(SUMDM_sez1R3)

N_Parcel_sez1R3 = IGSUM(N_Parcel_sez1R3)

Droplet_sez1R3 = GSUM(Droplet_sez1R3)

SUM_ND1_sez1R3 = GSUM(SUM_ND1_sez1R3)

SUM_ND2_sez1R3 = GSUM(SUM_ND2_sez1R3)

SUM_ND3_sez1R3 = GSUM(SUM_ND3_sez1R3)

SUM_D2_sez1R3 = GSUM(SUM_D2_sez1R3)

SUM_D3_sez1R3 = GSUM(SUM_D3_sez1R3)

SUM_MD2_sez1R3 = GSUM(SUM_MD2_sez1R3)

SUM_MD3_sez1R3 = GSUM(SUM_MD3_sez1R3)

SUM_VAVG_sez1R3 = GSUM(SUM_VAVG_sez1R3)

 

 

SUMDM_sez2R1 = GSUM(SUMDM_sez2R1)

N_Parcel_sez2R1 = IGSUM(N_Parcel_sez2R1)

Droplet_sez2R1 = GSUM(Droplet_sez2R1)

SUM_ND1_sez2R1 = GSUM(SUM_ND1_sez2R1)

SUM_ND2_sez2R1 = GSUM(SUM_ND2_sez2R1)

SUM_ND3_sez2R1 = GSUM(SUM_ND3_sez2R1)

SUM_D2_sez2R1 = GSUM(SUM_D2_sez2R1)

SUM_D3_sez2R1 = GSUM(SUM_D3_sez2R1)

SUM_MD2_sez2R1 = GSUM(SUM_MD2_sez2R1)

SUM_MD3_sez2R1 = GSUM(SUM_MD3_sez2R1)

SUM_VAVG_sez2R1 = GSUM(SUM_VAVG_sez2R1)

 

SUMDM_sez2R2 = GSUM(SUMDM_sez2R2)

N_Parcel_sez2R2 = IGSUM(N_Parcel_sez2R2)

Droplet_sez2R2 = GSUM(Droplet_sez2R2)

SUM_ND1_sez2R2 = GSUM(SUM_ND1_sez2R2)

SUM_ND2_sez2R2 = GSUM(SUM_ND2_sez2R2)

SUM_ND3_sez2R2 = GSUM(SUM_ND3_sez2R2)

SUM_D2_sez2R2 = GSUM(SUM_D2_sez2R2)

SUM_D3_sez2R2 = GSUM(SUM_D3_sez2R2)

SUM_MD2_sez2R2 = GSUM(SUM_MD2_sez2R2)

SUM_MD3_sez2R2 = GSUM(SUM_MD3_sez2R2)

SUM_VAVG_sez2R2 = GSUM(SUM_VAVG_sez2R2)

 

SUMDM_sez2R3 = GSUM(SUMDM_sez2R3)

N_Parcel_sez2R3 = IGSUM(N_Parcel_sez2R3)

Droplet_sez2R3 = GSUM(Droplet_sez2R3)

SUM_ND1_sez2R3 = GSUM(SUM_ND1_sez2R3)

SUM_ND2_sez2R3 = GSUM(SUM_ND2_sez2R3)

SUM_ND3_sez2R3 = GSUM(SUM_ND3_sez2R3)

SUM_D2_sez2R3 = GSUM(SUM_D2_sez2R3)

SUM_D3_sez2R3 = GSUM(SUM_D3_sez2R3)

SUM_MD2_sez2R3 = GSUM(SUM_MD2_sez2R3)

SUM_MD3_sez2R3 = GSUM(SUM_MD3_sez2R3)

SUM_VAVG_sez2R3 = GSUM(SUM_VAVG_sez2R3)

 

 

SUMDM_sez3R1 = GSUM(SUMDM_sez3R1)

N_Parcel_sez3R1 = IGSUM(N_Parcel_sez3R1)

Droplet_sez3R1 = GSUM(Droplet_sez3R1)

SUM_ND1_sez3R1 = GSUM(SUM_ND1_sez3R1)

SUM_ND2_sez3R1 = GSUM(SUM_ND2_sez3R1)

SUM_ND3_sez3R1 = GSUM(SUM_ND3_sez3R1)

SUM_D2_sez3R1 = GSUM(SUM_D2_sez3R1)

SUM_D3_sez3R1 = GSUM(SUM_D3_sez3R1)

SUM_MD2_sez3R1 = GSUM(SUM_MD2_sez3R1)

SUM_MD3_sez3R1 = GSUM(SUM_MD3_sez3R1)

SUM_VAVG_sez3R1 = GSUM(SUM_VAVG_sez3R1)

 

SUMDM_sez3R2 = GSUM(SUMDM_sez3R2)

N_Parcel_sez3R2 = IGSUM(N_Parcel_sez3R2)

Droplet_sez3R2 = GSUM(Droplet_sez3R2)

SUM_ND1_sez3R2 = GSUM(SUM_ND1_sez3R2)

SUM_ND2_sez3R2 = GSUM(SUM_ND2_sez3R2)

SUM_ND3_sez3R2 = GSUM(SUM_ND3_sez3R2)

SUM_D2_sez3R2 = GSUM(SUM_D2_sez3R2)

SUM_D3_sez3R2 = GSUM(SUM_D3_sez3R2)

SUM_MD2_sez3R2 = GSUM(SUM_MD2_sez3R2)

SUM_MD3_sez3R2 = GSUM(SUM_MD3_sez3R2)

SUM_VAVG_sez3R2 = GSUM(SUM_VAVG_sez3R2)

 

SUMDM_sez3R3 = GSUM(SUMDM_sez3R3)

N_Parcel_sez3R3 = IGSUM(N_Parcel_sez3R3)

Droplet_sez3R3 = GSUM(Droplet_sez3R3)

SUM_ND1_sez3R3 = GSUM(SUM_ND1_sez3R3)

SUM_ND2_sez3R3 = GSUM(SUM_ND2_sez3R3)

SUM_ND3_sez3R3 = GSUM(SUM_ND3_sez3R3)

SUM_D2_sez3R3 = GSUM(SUM_D2_sez3R3)

SUM_D3_sez3R3 = GSUM(SUM_D3_sez3R3)

SUM_MD2_sez3R3 = GSUM(SUM_MD2_sez3R3)

SUM_MD3_sez3R3 = GSUM(SUM_MD3_sez3R3)

SUM_VAVG_sez3R3 = GSUM(SUM_VAVG_sez3R3)

end if

 

 

 

c............................ DROPLET AVERAGING.....................................................

IF (SUMDM.gt.0) THEN

VMAVG = SUM_VAVG / Droplet

D1AVG = SUM_ND1 / Droplet

D32N = SUM_ND3 / SUM_ND2

D32 = SUM_D3 / SUM_D2

D32M = SUM_MD3 / SUM_MD2

 

ELSE

 

VMAVG = 0.0

D1AVG = 0.0

D32N = 0.0

D32 = 0.0

D32M = 0.0

END IF

 

 

IF (SUMDM_sez1.gt.0) THEN

VMAVG_sez1 = SUM_VAVG_sez1 / Droplet_sez1

D1AVG_sez1 = SUM_ND1_sez1 / Droplet_sez1

D32N_sez1 = SUM_ND3_sez1 / SUM_ND2_sez1

D32_sez1 = SUM_D3_sez1 / SUM_D2_sez1

D32M_sez1 = SUM_MD3_sez1 / SUM_MD2_sez1

 

ELSE

 

VMAVG_sez1 = 0.0

D1AVG_sez1 = 0.0

D32N_sez1 = 0.0

D32_sez1 = 0.0

D32M_sez1 = 0.0

END IF

 

IF (SUMDM_sez2.gt.0) THEN

VMAVG_sez2 = SUM_VAVG_sez2 / Droplet_sez2

D1AVG_sez2 = SUM_ND1_sez2 / Droplet_sez2

D32N_sez2 = SUM_ND3_sez2 / SUM_ND2_sez2

D32_sez2 = SUM_D3_sez2 / SUM_D2_sez2

D32M_sez2 = SUM_MD3_sez2 / SUM_MD2_sez2

ELSE

VMAVG_sez2 = 0.0

D1AVG_sez2 = 0.0

D32N_sez2 = 0.0

D32_sez2 = 0.0

D32M_sez2 = 0.0

END IF

 

IF (SUMDM_sez3.gt.0) THEN

VMAVG_sez3 = SUM_VAVG_sez3 / Droplet_sez3

D1AVG_sez3 = SUM_ND1_sez3 / Droplet_sez3

D32N_sez3 = SUM_ND3_sez3 / SUM_ND2_sez3

D32_sez3 = SUM_D3_sez3 / SUM_D2_sez3

D32M_sez3 = SUM_MD3_sez3 / SUM_MD2_sez3

 

ELSE

 

VMAVG_sez3 = 0.0

D1AVG_sez3 = 0.0

D32N_sez3 = 0.0

D32_sez3 = 0.0

D32M_sez3 = 0.0

END IF

 

 

IF (SUMDM_sez1R1.gt.0) THEN

VMAVG_sez1R1 = SUM_VAVG_sez1R1 / Droplet_sez1R1

D1AVG_sez1R1 = SUM_ND1_sez1R1 / Droplet_sez1R1

D32N_sez1R1 = SUM_ND3_sez1R1 / SUM_ND2_sez1R1

D32_sez1R1 = SUM_D3_sez1R1 / SUM_D2_sez1R1

D32M_sez1R1 = SUM_MD3_sez1R1 / SUM_MD2_sez1R1

ELSE

VMAVG_sez1R1 = 0.0

D1AVG_sez1R1 = 0.0

D32N_sez1R1 = 0.0

D32_sez1R1 = 0.0

D32M_sez1R1 = 0.0

END IF

 

IF (SUMDM_sez1R2.gt.0) THEN

VMAVG_sez1R2 = SUM_VAVG_sez1R2 / Droplet_sez1R2

D1AVG_sez1R2 = SUM_ND1_sez1R2 / Droplet_sez1R2

D32N_sez1R2 = SUM_ND3_sez1R2 / SUM_ND2_sez1R2

D32_sez1R2 = SUM_D3_sez1R2 / SUM_D2_sez1R2

D32M_sez1R2 = SUM_MD3_sez1R2 / SUM_MD2_sez1R2

ELSE

VMAVG_sez1R2 = 0.0

D1AVG_sez1R2 = 0.0

D32N_sez1R2 = 0.0

D32_sez1R2 = 0.0

D32M_sez1R2 = 0.0

END IF

 

 

IF (SUMDM_sez1R3.gt.0) THEN

VMAVG_sez1R3 = SUM_VAVG_sez1R3 / Droplet_sez1R3

D1AVG_sez1R3 = SUM_ND1_sez1R3 / Droplet_sez1R3

D32N_sez1R3 = SUM_ND3_sez1R3 / SUM_ND2_sez1R3

D32_sez1R3 = SUM_D3_sez1R3 / SUM_D2_sez1R3

D32M_sez1R3 = SUM_MD3_sez1R3 / SUM_MD2_sez1R3

ELSE

VMAVG_sez1R3 = 0.0

D1AVG_sez1R3 = 0.0

D32N_sez1R3 = 0.0

D32_sez1R3 = 0.0

D32M_sez1R3 = 0.0

 

END IF

 

IF (SUMDM_sez2R1.gt.0) THEN

 

VMAVG_sez2R1 = SUM_VAVG_sez2R1 / Droplet_sez2R1

D1AVG_sez2R1 = SUM_ND1_sez2R1 / Droplet_sez2R1

D32N_sez2R1 = SUM_ND3_sez2R1 / SUM_ND2_sez2R1

D32_sez2R1 = SUM_D3_sez2R1 / SUM_D2_sez2R1

D32M_sez2R1 = SUM_MD3_sez2R1 / SUM_MD2_sez2R1

ELSE

VMAVG_sez2R1 = 0.0

D1AVG_sez2R1 = 0.0

D32N_sez2R1 = 0.0

D32_sez2R1 = 0.0

D32M_sez2R1 = 0.0

END IF

 

IF (SUMDM_sez2R2.gt.0) THEN

VMAVG_sez2R2 = SUM_VAVG_sez2R2 / Droplet_sez2R2

D1AVG_sez2R2 = SUM_ND1_sez2R2 / Droplet_sez2R2

D32N_sez2R2 = SUM_ND3_sez2R2 / SUM_ND2_sez2R2

D32_sez2R2 = SUM_D3_sez2R2 / SUM_D2_sez2R2

D32M_sez2R2 = SUM_MD3_sez2R2 / SUM_MD2_sez2R2

ELSE

VMAVG_sez2R2 = 0.0

D1AVG_sez2R2 = 0.0

D32N_sez2R2 = 0.0

D32_sez2R2 = 0.0

D32M_sez2R2 = 0.0

END IF

 

IF (SUMDM_sez2R3.gt.0) THEN

VMAVG_sez2R3 = SUM_VAVG_sez2R3 / Droplet_sez2R3

D1AVG_sez2R3 = SUM_ND1_sez2R3 / Droplet_sez2R3

D32N_sez2R3 = SUM_ND3_sez2R3 / SUM_ND2_sez2R3

D32_sez2R3 = SUM_D3_sez2R3 / SUM_D2_sez2R3

D32M_sez2R3 = SUM_MD3_sez2R3 / SUM_MD2_sez2R3

ELSE

VMAVG_sez2R3 = 0.0

D1AVG_sez2R3 = 0.0

D32N_sez2R3 = 0.0

D32_sez2R3 = 0.0

D32M_sez2R3 = 0.0

END IF

 

IF (SUMDM_sez3R1.gt.0) THEN

 

VMAVG_sez3R1 = SUM_VAVG_sez3R1 / Droplet_sez3R1

D1AVG_sez3R1 = SUM_ND1_sez3R1 / Droplet_sez3R1

D32N_sez3R1 = SUM_ND3_sez3R1 / SUM_ND2_sez3R1

D32_sez3R1 = SUM_D3_sez3R1 / SUM_D2_sez3R1

D32M_sez3R1 = SUM_MD3_sez3R1 / SUM_MD2_sez3R1

ELSE

VMAVG_sez3R1 = 0.0

D1AVG_sez3R1 = 0.0

D32N_sez3R1 = 0.0

D32_sez3R1 = 0.0

D32M_sez3R1 = 0.0

END IF

 

IF (SUMDM_sez3R2.gt.0) THEN

VMAVG_sez3R2 = SUM_VAVG_sez3R2 / Droplet_sez3R2

D1AVG_sez3R2 = SUM_ND1_sez3R2 / Droplet_sez3R2

D32N_sez3R2 = SUM_ND3_sez3R2 / SUM_ND2_sez3R2

D32_sez3R2 = SUM_D3_sez3R2 / SUM_D2_sez3R2

D32M_sez3R2 = SUM_MD3_sez3R2 / SUM_MD2_sez3R2

ELSE

VMAVG_sez3R2 = 0.0

D1AVG_sez3R2 = 0.0

D32N_sez3R2 = 0.0

D32_sez3R2 = 0.0

D32M_sez3R2 = 0.0

 

END IF

 

IF (SUMDM_sez3R3.gt.0) THEN

VMAVG_sez3R3 = SUM_VAVG_sez3R3 / Droplet_sez3R3

D1AVG_sez3R3 = SUM_ND1_sez3R3 / Droplet_sez3R3

D32N_sez3R3 = SUM_ND3_sez3R3 / SUM_ND2_sez3R3

D32_sez3R3 = SUM_D3_sez3R3 / SUM_D2_sez3R3

D32M_sez3R3 = SUM_MD3_sez3R3 / SUM_MD2_sez3R3

ELSE

VMAVG_sez3R3 = 0.0

D1AVG_sez3R3 = 0.0

D32N_sez3R3 = 0.0

D32_sez3R3 = 0.0

D32M_sez3R3 = 0.0

 

END IF

 

 

C***************************************************************

c********************Stima della penetrazione*******************

C***************************************************************

IF (cone .eq. 1) then

 

PENEMAX = 0.0d0

PENEMAXM=0.0d0

MAXR= 0.0

MINR= 0.0

DO 115 N=1,NDR

IF( IDC(N).gt.0 )THEN

icsys = iniettore

POSGLO(1)=DRC (N,1)

POSGLO(2)=DRC (N,2)

POSGLO(3)=DRC (N,3)

CALL GLOLOC(POSGLO, POSLOC, ICSYS)

c MAXR=posloc(3)*tan(teta)

MAXR=0.006

c MINR=-posloc(3)*tan(teta)

if (posloc(2) .lt. MAXR) then

c if (posloc(2) .gt. MINR) then

if (posloc(3) .gt. PENEMAX) then

PENEMAX = posloc(3)

c endif

endif

endif

ENDIF

115 continue

ELSE

 

PENEMAX = 0.0d0

PENEPRE = 0.0d0

DO 116 N=1,NDR

IF( IDC(N).gt.0 )THEN

icsys = iniettore

POSGLO(1)=DRC (N,1)

POSGLO(2)=DRC (N,2)

POSGLO(3)=DRC (N,3)

CALL GLOLOC(POSGLO, POSLOC, ICSYS)

if (DRD(N) .le. (2*D1AVG)) then

c POSGLO(1)=DRC (N,1)

c POSGLO(2)=DRC (N,2)

c POSGLO(3)=DRC (N,3)

c CALL GLOLOC(POSGLO, POSLOC, ICSYS)

if (posloc(3) .gt. PENEMAX .and.

& posloc(3).lt.penetol) then

PENEMAX = posloc(3)

endif

endif

 

c if (DRD(N) .gt. D1AVG) then

c POSGLO(1)=DRC (N,1)

c POSGLO(2)=DRC (N,2)

c POSGLO(3)=DRC (N,3)

c CALL GLOLOC(POSGLO, POSLOC, ICSYS)

if (posloc(3) .gt. PENEPRE) then

PENEPRE = posloc(3)

c endif

endif

endif

116 continue

 

ENDIF

 

 

IF (NHPC .GT. 1) THEN

PENEMAX = GMAX(PENEMAX)

PENEMAXM = PENEMAX * 1000

if (cone .eq. 0) then

PENEPRE = GMAX(PENEPRE)

PENEPREM = PENEPRE * 1000

endif

ENDIF

 

C***************************************************************

 

 

 

C*************** DROPLETS OUTPUT FILE*************************

 

c NAME = 'penetrazione.dat'

 

if (cone .eq. 1) then

 

CALL APPEND(88,'penetrazione.dat')

WRITE(88,3144)ITER,TIME,PENEMAX,PENEMAXM

3144 FORMAT (I6,2X,E12.5,3X,2(2x,E12.5))

CLOSE(88)

 

else

 

CALL APPEND(89,'penetrazione.dat')

WRITE(89,3145)ITER,TIME,PENEMAX,PENEMAXM,

& PENEPRE,PENEPREM

3145 FORMAT (I6,2X,E12.5,3X,4(2x,E12.5))

CLOSE(89)

 

endif

 

 

c NAME = 'posdat-tot.dat'

CALL APPEND(90,'posdat-tot.dat')

if (iter .eq. 0.0) then

WRITE(90,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (90,3146) ITER,TIME,N_Parcel,

& SUMDM,D32N,D1AVG,D32,D32M,

& VMAVG

3146 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(90)

endif

 

c NAME = 'posdat-sez1.dat'

CALL APPEND(91,'posdat-sez1.dat')

if (iter .eq. 0.0) then

WRITE(91,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (91,3147) ITER,TIME,N_Parcel_sez1,

& SUMDM_sez1,D32N_sez1,D1AVG_sez1,D32_sez1,D32M_sez1,

& VMAVG_sez1

3147 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(91)

endif

 

if (R1 .ne. 0) then

c NAME = 'posdat-sez1R1.dat'

CALL APPEND(191,'posdat-sez1R1.dat')

if (iter .eq. 0.0) then

WRITE(191,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (191,4147) ITER,TIME,N_Parcel_sez1R1,

& SUMDM_sez1R1,D32N_sez1R1,D1AVG_sez1R1,D32_sez1R1,D32M_sez1R1,

& VMAVG_sez1R1

4147 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(191)

endif

endif

 

if (R2 .ne. 0) then

c NAME = 'posdat-sez1R2.dat'

CALL APPEND(192,'posdat-sez1R2.dat')

if (iter .eq. 0.0) then

WRITE(192,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (192,4148) ITER,TIME,N_Parcel_sez1R2,

& SUMDM_sez1R2,D32N_sez1R2,D1AVG_sez1R2,D32_sez1R2,D32M_sez1R2,

& VMAVG_sez1R2

4148 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(192)

endif

endif

 

if (R3 .ne. 0) then

c NAME = 'posdat-sez1R3.dat'

CALL APPEND(193,'posdat-sez1R3.dat')

if (iter .eq. 0.0) then

WRITE(193,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (193,4149) ITER,TIME,N_Parcel_sez1R3,

& SUMDM_sez1R3,D32N_sez1R3,D1AVG_sez1R3,D32_sez1R3,D32M_sez1R3,

& VMAVG_sez1R3

4149 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(193)

endif

endif

 

 

 

c NAME = 'posdat-sez2.dat'

CALL APPEND(92,'posdat-sez2.dat')

if (iter .eq. 0) then

WRITE(92,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (92,3148) ITER,TIME,N_Parcel_sez2,

& SUMDM_sez2,D32N_sez2,D1AVG_sez2,D32_sez2,D32M_sez2,

& VMAVG_sez2

3148 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(92)

endif

 

if (R1 .ne. 0) then

c NAME = 'posdat-sez2R1.dat'

CALL APPEND(291,'posdat-sez2R1.dat')

if (iter .eq. 0.0) then

WRITE(291,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (291,5147) ITER,TIME,N_Parcel_sez2R1,

& SUMDM_sez2R1,D32N_sez2R1,D1AVG_sez2R1,D32_sez2R1,D32M_sez2R1,

& VMAVG_sez2R1

5147 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(291)

endif

endif

 

if (R2 .ne. 0) then

c NAME = 'posdat-sez2R2.dat'

CALL APPEND(292,'posdat-sez2R2.dat')

if (iter .eq. 0.0) then

WRITE(292,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (292,5148) ITER,TIME,N_Parcel_sez2R2,

& SUMDM_sez2R2,D32N_sez2R2,D1AVG_sez2R2,D32_sez2R2,D32M_sez2R2,

& VMAVG_sez2R2

5148 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(292)

endif

endif

 

if (R3 .ne. 0) then

c NAME = 'posdat-sez2R3.dat'

CALL APPEND(293,'posdat-sez2R3.dat')

if (iter .eq. 0.0) then

WRITE(293,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (293,5149) ITER,TIME,N_Parcel_sez2R3,

& SUMDM_sez2R3,D32N_sez2R3,D1AVG_sez2R3,D32_sez2R3,D32M_sez2R3,

& VMAVG_sez2R3

5149 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(293)

endif

endif

 

c NAME = 'posdat-sez3.dat'

CALL APPEND(93,'posdat-sez3.dat')

if (iter .eq. 0) then

WRITE(93,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (93,3149) ITER,TIME,N_Parcel_sez3,

& SUMDM_sez3,D32N_sez3,D1AVG_sez3,D32_sez3,D32M_sez3,

& VMAVG_sez3

3149 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(93)

endif

 

if (R1 .ne. 0) then

c NAME = 'posdat-sez3R1.dat'

CALL APPEND(391,'posdat-sez3R1.dat')

if (iter .eq. 0.0) then

WRITE(391,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (391,6147) ITER,TIME,N_Parcel_sez3R1,

& SUMDM_sez3R1,D32N_sez3R1,D1AVG_sez3R1,D32_sez3R1,D32M_sez3R1,

& VMAVG_sez3R1

6147 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(391)

endif

endif

 

if (R2 .ne. 0) then

c NAME = 'posdat-sez3R2.dat'

CALL APPEND(392,'posdat-sez3R2.dat')

if (iter .eq. 0.0) then

WRITE(392,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (392,6148) ITER,TIME,N_Parcel_sez3R2,

& SUMDM_sez3R2,D32N_sez3R2,D1AVG_sez3R2,D32_sez3R2,D32M_sez3R2,

& VMAVG_sez3R2

6148 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(392)

endif

endif

 

if (R3 .ne. 0) then

c NAME = 'posdat-sez3R3.dat'

CALL APPEND(393,'posdat-sez3R3.dat')

if (iter .eq. 0.0) then

WRITE(393,*)'!T-STEP TIME TOT-PARC-NUM TOT-MASS(kg)',

&' D32N D1AVG D32 D32M',

&' VMAVG(m/s)'

else

WRITE (393,6149) ITER,TIME,N_Parcel_sez3R3,

& SUMDM_sez3R3,D32N_sez3R3,D1AVG_sez3R3,D32_sez3R3,D32M_sez3R3,

& VMAVG_sez3R3

6149 FORMAT (I6,2X,E12.5,3X,I9,5X,6(1X,E12.5))

CLOSE(393)

endif

endif

 

 

endif

return

end

 

C********* CALCULATION OF DROPLETS QUANTITIES * STOP *******************

 

 

 

 

C**********************************************************************C

SUBROUTINE APPEND (IU,FNAME)

C**********************************************************************C

 

c Open the file or append to the file

 

C-----------------------------------------------------------------------

CHARACTER*(*) FNAME

CHARACTER CHARS*1

LOGICAL EXISTS

c

INQUIRE (FILE=FNAME,EXIST=EXISTS)

IF (EXISTS) THEN

OPEN (IU,FILE=FNAME,STATUS='OLD')

c WRITE (61,1010) FNAME

100 READ (IU,1000,END=200,ERR=900) CHARS

GO TO 100

ELSE

OPEN (IU,FILE=FNAME,STATUS='NEW')

c WRITE (61,1020) FNAME

END IF

200 CONTINUE

RETURN

900 PRINT *,'***ERROR-POSDAT*** Cannot read data file: ',FNAME

RETURN

1000 FORMAT (A1)

1010 FORMAT (1X,'***APPEND1*** Append to the data file: ',A)

1020 FORMAT (1X,'***APPEND2*** Open a new data file: ',A)

END

[apix_1024]

tutto chiaro vero?? per me no :2funny::asd::2funny:

[penguin86]

ahahahhahahahahh magari, ora ti spiego XD

Praticamente ho cambiato scietà di atletica e con quella attuale con me c'è una ragazza molto carina.

Tempo fa ad una corsa vicino al mio allenatore (ora ex) c'era un uomo che non conoscevo, così al mio allenatore (ora ex) faccio: ammazza quella che carina, sopratutto un suo lato

L'uomo accanto risponde: ci credo è mia nipote

 

E d'altronte ora quell'uomo è il mio allenatore

 

oggi dopo aver fatto 3,5km di riscaldamento e 10*200m senza riposo, in pratica ritornavi in dietro correndo più piano e ripartivi, dopo due giri defaticanti che io ho fatto affianco a lei, chissà perchè mi fà:

se ti piace tanto stare in compagnia di mia nipote perchè non ci metti a nche 30min di fondo? (corsa lenta) :muro:

Io: vabbhè, almeno sto in buona compagnia

 

Tutti gli altri a farsi la doccia, alc alduccioe Noi due al massacro

 

Vabbhè me la sono cercata

 

Uahahahah! Questo è un mitoooooo!

 

Se ne vale la pena, hai fatto bene un po' di sofferenza in questi casi non fa' male.

 

Questo è vero.

 

C*************************************************************************

SUBROUTINE

 

[...]

 

END

 

tutto chiaro vero?? per me no :2funny::asd::2funny:

Uhm... no... xke usi gli if in cascata invece di un case?

[apix_1024]

Uhm... no... xke usi gli if in cascata invece di un case?

 

 

perchè me l'ha passata così il prof che l'ha fatta ed io voglia di metterla a posto nn ne ho

[swattolo]

Uahahahah! Questo è un mitoooooo!

 

 

 

Questo è vero.

 

 

 

 

Uhm... no... xke usi gli if in cascata invece di un case?

 

Ma che mito, a me mi ammazza

 

Oggi la mia ex si è dimostrata più affettuosa del soltio :cheazz:

Non lo so nemmeno io cosa voglio...

[antobov]

Ma che mito, a me mi ammazza

 

Oggi la mia ex si è dimostrata più affettuosa del soltio :cheazz:

Non lo so nemmeno io cosa voglio...

 

Magari ha bisogno qualcosa.. :perfido:

 

mollerei tutte queste paranoie e mi concentrerei sulla Nipote..

[swattolo]

Magari ha bisogno qualcosa.. :perfido:

 

mollerei tutte queste paranoie e mi concentrerei sulla Nipote..

 

antobov è una brava ragazza

 

Prima mi dimentico dela ex

Poi esco con il "mio angioletto", la ragazza conosciuta alla festa

poi con sorpresa scopro che la nipote ed io siamo della stess squadra

Poi la mia ex oggi si è comportata in modo strano ed io gli "davo corda"

 

Che casino...:(

[Rankine]

e io cosa devo dire che in agosto mi sono preso una fulminata per una bellissima ragazza che poi si è rivelata una mia cugina di secondo grado?!?!!?!!:cheazz:

[ShadoW]

e io cosa devo dire che in agosto mi sono preso una fulminata per una bellissima ragazza che poi si è rivelata una mia cugina di secondo grado?!?!!?!!:cheazz:

 

 

ahahahah

 

 

 

ragaaaaa

 

guardati questo video :D:D:D

 

YouTube - Dancing Bear

[swattolo]

e io cosa devo dire che in agosto mi sono preso una fulminata per una bellissima ragazza che poi si è rivelata una mia cugina di secondo grado?!?!!?!!:cheazz:

 

Come è andata a finire ?

Ora mi trovo su msn un mex di lei che mi fà: oggi vedendo la mia amica con il suo raga ho pensato sempre a te, sei una persona molto speciale e non sai quanto di posso voler bene. Tua tata, 4 Ever :(

[Le085]

Come è andata a finire ?

Ora mi trovo su msn un mex di lei che mi fà: oggi vedendo la mia amica con il suo raga ho pensato sempre a te, sei una persona molto speciale e non sai quanto di posso voler bene. Tua tata, 4 Ever :(

 

la tua ex?

o quella dell'atletica?

o quella della sagra?

o quella che dici è l'angioletto tuo?

o quella del tuo avatar?

o quella dell'altra sera?

o quell'altra ancora dell'atletica?

o quella di messenger?

certo che se tu hai le idee confuse figuriamoci noi:2funny:

[swattolo]

la tua ex?

o quella dell'atletica?

o quella della sagra?

o quella che dici è l'angioletto tuo?

o quella del tuo avatar?

o quella dell'altra sera?

o quell'altra ancora dell'atletica?

o quella di messenger?

certo che se tu hai le idee confuse figuriamoci noi:2funny:

 

Ok la mia ex mi ha detto che vuole rimettersi con me :(

prima di andar a far cena (1omin fà) mi ha scritto:

Domani alle 11 troviamoci dietro la scuola, ti farò capire cosa provo veramente per te :( :(

 

Ma porca miseria, non so neppure io quello che fare :( :cry:

[Le085]

Domani alle 11 troviamoci dietro la scuola, ti farò capire cosa provo veramente per te :( :(

 

 

 

beh io ci andrei solo per curiosità di quello che ha intenzione di fare

poi se vedi che lo fa apposta perchè ti ha visto con l'altra la tratterei in malo modo:asd:

+ che altro perchè non mi sembra giusto che tu debba stare ai suoi comodi...cioè finchè tu ci stavi ancora sotto lei non ne voleva sapere ora che finalmente ti eri riuscito a distaccare un po' lei torna alla carica?

[swattolo]

 

 

 

beh io ci andrei solo per curiosità di quello che ha intenzione di fare

poi se vedi che lo fa apposta perchè ti ha visto con l'altra la tratterei in malo modo:asd:

+ che altro perchè non mi sembra giusto che tu debba stare ai suoi comodi...cioè finchè tu ci stavi ancora sotto lei non ne voleva sapere ora che finalmente ti eri riuscito a distaccare un po' lei torna alla carica?

 

Leo posso immaginare quindi... vedete di non pensrae male perchè è sempre una brava ragazza ^^

 

Sò cosa è successo e ci ho riflettuto anche io, però comunque non so cosa fare, cioè se gli direi di si io starei benissimo, certo rinuncierei alle altre ma non importa.

Spero solo che abbia imparato dal passato :(

[Le085]

Leo posso immaginare quindi... vedete di non pensrae male perchè è sempre una brava ragazza ^^

 

Sò cosa è successo e ci ho riflettuto anche io, però comunque non so cosa fare, cioè se gli direi di si io starei benissimo, certo rinuncierei alle altre ma non importa.

Spero solo che abbia imparato dal passato :(

lo spero per te, di sicuro te sai meglio di noi che ragazza è...

[swattolo]

Sisi per quelel cose garantisco al 100%, è proprio una brava ragazza, non è che se gli metti le mani li ti fà tutto quello che vuole, siamo chiari...

[Le085]

Sisi per quelel cose garantisco al 100%, è proprio una brava ragazza, non è che se gli metti le mani li ti fà tutto quello che vuole, siamo chiari...

ah però lo vuole...solo che nn lo fa:P

 

scherzo

 

scusa anzi, che se poi ti ci devi rimettere insieme mi pare anche brutto offenderla così gratuitamente...:(

 

dai speriamo che finalmente vi chiarite, almeno smetti anche di sbavare dietro ogni ragazza che passa:asd:

 

domani alle 13 e 30 tutti sintonizzati sul thread dello spam per vedere la 235° puntata di "swattolo in love":uglystupid2: