SPAM !!!!!!!!!!!!!!!! [THREAD UFFICIALE!!!]
-
Rankine ha scritto:
ma si, anch'io, però una porsche gt2 con un impianto audio da 5000W non me la toglie nessuno, una villetta fuori dai maroni e il resto me lo godo in viaggi..... e probabilmente comprerei la pizzeria dove lavoro e ci metterei 4 o 5 dipendenti. Comunque lavorerei anch'io, perchè certamente non sto studiando come un animale per nientePer le auto potrei dire anche una lamborghine o altro perchè ipoteticamente potrei permettermela, però non mi importerebbe.
Una fuori-serie ma niente di che, per i motivi citati sopra
-
Rankine ha scritto:
ma si, anch'io, però una porsche gt2 con un impianto audio da 5000W non me la toglie nessuno, una villetta fuori dai maroni e il resto me lo godo in viaggi..... e probabilmente comprerei la pizzeria dove lavoro e ci metterei 4 o 5 dipendenti. Comunque lavorerei anch'io, perchè certamente non sto studiando come un animale per nientecioè tu vincendo palate di soldi ti prenderesti una gt2?? come il ranzani di turno???
qui in europa io mi presenteri con questa:
leggi leggi qui: http://paultan.org/archives/2005/08/13/nissan-skyline-r34-gt-r-z-tune/
:fiufiu: -
beh, la gt2 è l'auto dei miei sogni, ma pernso che potrei spendere un paio di milioni di euro per farmi il parco macchine....
-
8c competizione per me
suv per mamma
casa al mare per mamma
casa in montagna per me e papà
quad per papà
un 20000 euro di hardware per benchare
basta così
il resto probabilmente andrebbe in tasse e benzina
-
io mi farò fare il pc direttamente dai fornitori di pc per la NASA,hahahahah!!!
-
ho giocato il sistema....ma non mi sono risultate delle sestine così schifose dal simulatore di oggi.... mah, sperem!!!:boh:
-
ecco...non ho vinto una mazza.... ho centrato 2 numeri...
-
Rankine ha scritto:
ecco...non ho vinto una mazza.... ho centrato 2 numeri...bisogna dire che è stato un ottimo sistema per perdere 52 euro:clapclap:
beh consolati....con 52 euro di hardware non ci usciva nulla:asd:
-
Le085 ha scritto:
bisogna dire che è stato un ottimo sistema per perdere 52 euro:clapclap:beh consolati....con 52 euro di hardware non ci usciva nulla:asd:
Questo lo dici tu
-
swattolo ha scritto:
Questo lo dici tu
a te si ma a lui no
-
Le085 ha scritto:
a te si ma a lui no
Come no, prende una versione tarocca made in China di 5° mano e ci prende anche un nehalem con 50€
:2funny: -
-
Adesso mi guardo il terzo episodio di Fringe. Finalmente una serata tranquilla.....
-
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
-
swattolo ha scritto:
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......
-
antobov ha scritto:
: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
-
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.
-
antobov ha scritto:
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: -
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!!!!!!)
-
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
Ciao! Sembra che tu sia interessato a questa conversazione, ma non hai ancora un account.
Stanco di dover scorrere gli stessi post a ogni visita? Quando registri un account, tornerai sempre esattamente dove eri rimasto e potrai scegliere di essere avvisato delle nuove risposte (tramite email o notifica push). Potrai anche salvare segnalibri e votare i post per mostrare il tuo apprezzamento agli altri membri della comunità.
Con il tuo contributo, questo post potrebbe essere ancora migliore 💗
Registrati Accedi