Extracting information on a particle track
Hi!
I would like to get the temperature in each cell that a particle track passes through. Is that possible? Regards Bo |
Re: Extracting information on a particle track
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) |
All times are GMT -4. The time now is 00:19. |