Hi!

I would like to get the temperature in each cell that a particle track passes through. Is that possible?

Regards

Bo

Hi Bo,

Yes it is. The variable idc(idr) gives you the cell number a particle is in. You can print these data from the cell number. If idc.eq.0 the particle is stuck in the domain, If idc.lt.0 it has left the domain.

Michiel

This was sent to mee as an example in POSDAT.f by starcd

COMMON /C33001/ IDC

C --- DROPLET DATA

OPEN (UNIT=100,FILE='domain.out')

OPEN (UNIT=101,FILE='evap.out')

OPEN (UNIT=102,FILE='stuck.out')

REWIND (100)

REWIND (101)

REWIND (102)

MASSWO = 0.0

HEATWO = 0.0

DO IDR=1,NDR C --- DROPLET (IDR) IS LOCATED IN CELL (IC)

IC=IDC(IDR)

MASFRA = DRMAIC(IDR,1)/DRM(IDR)

IF (IC.GT.0) THEN

WRITE (100,911) IDR,IDT(IDR),DRD(IDR),DRM(IDR),DRT(IDR),MASFRA

ELSE IF (IC.EQ.0) THEN

WRITE (101,911) IDR,IDT(IDR),DRD(IDR),DRM(IDR),DRT(IDR),MASFRA

ELSE

WRITE (102,911) IDR,IDT(IDR),DRD(IDR),DRM(IDR),DRT(IDR),MASFRA

ENDIF

MASSWO = DRMAIC(IDR,1) + MASSWO

HEATWO = DRT(IDR)*

&(4181.72*DRMAIC(IDR,1)+1750*(DRM(IDR)-DRMAIC(IDR,1)))+HEATWO

ENDDO

CLOSE (100)

CLOSE (101)

CLOSE (102)