Changeset 5488 for LMDZ6/branches/contrails/libf/phylmd/lmdz_aviation.f90

Timestamp:

Jan 17, 2025, 6:12:48 PM (13 days ago)

Author:

aborella

Message:

Added dependance of rei to contrail fraction

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_aviation.f90 (modified) (1 diff)

LMDZ6/branches/contrails/libf/phylmd/lmdz_aviation.f90

@@ -678 +678 @@
 !  INCLUDE 'netcdf.inc'
 !
-!  !--------------------------------------------------------
-!  !--input variables 
-!  !--------------------------------------------------------
-!  LOGICAL, INTENT(IN) :: debut
-!  REAL, INTENT(IN)    :: pphis(klon), pplay(klon,klev), paprs(klon,klev+1), t_seri(klon,klev)
-!
-!  !--------------------------------------------------------
-!  !    ... Local variables
-!  !--------------------------------------------------------
-!
-!  CHARACTER (LEN=20) :: modname='airplane_mod'
-!  INTEGER :: i, k, kori, iret, varid, error, ncida, klona
-!  INTEGER,SAVE :: nleva, ntimea
-!!$OMP THREADPRIVATE(nleva,ntimea)
-!  REAL, ALLOCATABLE :: pkm_airpl_glo(:,:,:)    !--km/s
-!  REAL, ALLOCATABLE :: ph2o_airpl_glo(:,:,:)   !--molec H2O/cm3/s
-!  REAL, ALLOCATABLE, SAVE :: zmida(:), zinta(:)
-!  REAL, ALLOCATABLE, SAVE :: pkm_airpl(:,:,:)
-!  REAL, ALLOCATABLE, SAVE :: ph2o_airpl(:,:,:)
-!!$OMP THREADPRIVATE(pkm_airpl,ph2o_airpl,zmida,zinta)
-!  REAL :: zalt(klon,klev+1)
-!  REAL :: zrho, zdz(klon,klev), zfrac
-!
-!  !
-!  IF (debut) THEN
-!  !--------------------------------------------------------------------------------
-!  !       ... Open the file and read airplane emissions
-!  !--------------------------------------------------------------------------------
-!  !
-!  IF (is_mpi_root .AND. is_omp_root) THEN
-!      !
-!      iret = nf_open('aircraft_phy.nc', 0, ncida)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to open aircraft_phy.nc file',1)
-!      ! ... Get lengths
-!      iret = nf_inq_dimid(ncida, 'time', varid)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get time dimid in aircraft_phy.nc file',1)
-!      iret = nf_inq_dimlen(ncida, varid, ntimea)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get time dimlen aircraft_phy.nc file',1)
-!      iret = nf_inq_dimid(ncida, 'vector', varid)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get vector dimid aircraft_phy.nc file',1)
-!      iret = nf_inq_dimlen(ncida, varid, klona)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get vector dimlen aircraft_phy.nc file',1)
-!      iret = nf_inq_dimid(ncida, 'lev', varid)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get lev dimid aircraft_phy.nc file',1)
-!      iret = nf_inq_dimlen(ncida, varid, nleva)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get lev dimlen aircraft_phy.nc file',1)
-!      !
-!      IF ( klona /= klon_glo ) THEN
-!        WRITE(lunout,*) 'klona & klon_glo =', klona, klon_glo
-!        CALL abort_physic(modname,'problem klon in aircraft_phy.nc file',1)
-!      ENDIF
-!      !
-!      IF ( ntimea /= 12 ) THEN
-!        WRITE(lunout,*) 'ntimea=', ntimea
-!        CALL abort_physic(modname,'problem ntime<>12 in aircraft_phy.nc file',1)
-!      ENDIF
-!      !
-!      ALLOCATE(zmida(nleva), STAT=error)
-!      IF (error /= 0) CALL abort_physic(modname,'problem to allocate zmida',1)
-!      ALLOCATE(pkm_airpl_glo(klona,nleva,ntimea), STAT=error)
-!      IF (error /= 0) CALL abort_physic(modname,'problem to allocate pkm_airpl_glo',1)
-!      ALLOCATE(ph2o_airpl_glo(klona,nleva,ntimea), STAT=error)
-!      IF (error /= 0) CALL abort_physic(modname,'problem to allocate ph2o_airpl_glo',1)
-!      !
-!      iret = nf_inq_varid(ncida, 'lev', varid)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get lev dimid aircraft_phy.nc file',1)
-!      iret = nf_get_var_double(ncida, varid, zmida)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to read zmida file',1)
-!      !
-!      iret = nf_inq_varid(ncida, 'emi_co2_aircraft', varid)  !--CO2 as a proxy for m flown - 
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get emi_distance dimid aircraft_phy.nc file',1)
-!      iret = nf_get_var_double(ncida, varid, pkm_airpl_glo)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to read pkm_airpl file',1)
-!      !
-!      iret = nf_inq_varid(ncida, 'emi_h2o_aircraft', varid)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to get emi_h2o_aircraft dimid aircraft_phy.nc file',1)
-!      iret = nf_get_var_double(ncida, varid, ph2o_airpl_glo)
-!      IF (iret /= NF_NOERR) CALL abort_physic(modname,'problem to read ph2o_airpl file',1)
-!      ! 
-!     ENDIF    !--is_mpi_root and is_omp_root
-!     !
-!     CALL bcast(nleva)
-!     CALL bcast(ntimea)
-!     ! 
-!     IF (.NOT.ALLOCATED(zmida)) ALLOCATE(zmida(nleva), STAT=error)
-!     IF (.NOT.ALLOCATED(zinta)) ALLOCATE(zinta(nleva+1), STAT=error)
-!     !
-!     ALLOCATE(pkm_airpl(klon,nleva,ntimea))
-!     ALLOCATE(ph2o_airpl(klon,nleva,ntimea))
-!     !
-!     ALLOCATE(flight_m(klon,klev))
-!     ALLOCATE(flight_h2o(klon,klev))
-!     !
-!     CALL bcast(zmida)
-!     zinta(1)=0.0                                   !--surface
-!     DO k=2, nleva
-!       zinta(k) = (zmida(k-1)+zmida(k))/2.0*1000.0  !--conversion from km to m
-!     ENDDO
-!     zinta(nleva+1)=zinta(nleva)+(zmida(nleva)-zmida(nleva-1))*1000.0 !--extrapolation for last interface
-!     !print *,'zinta=', zinta
-!     !
-!     CALL scatter(pkm_airpl_glo,pkm_airpl)
-!     CALL scatter(ph2o_airpl_glo,ph2o_airpl)
-!     !
-!!$OMP MASTER
-!     IF (is_mpi_root .AND. is_omp_root) THEN
-!        DEALLOCATE(pkm_airpl_glo)
-!        DEALLOCATE(ph2o_airpl_glo)
-!     ENDIF   !--is_mpi_root
-!!$OMP END MASTER
-!
-!  ENDIF !--debut
-!!
-!!--compute altitude of model level interfaces
-!!
-!  DO i = 1, klon
-!    zalt(i,1)=pphis(i)/RG         !--in m
-!  ENDDO
-!!
-!  DO k=1, klev
-!    DO i = 1, klon
-!      zrho=pplay(i,k)/t_seri(i,k)/RD
-!      zdz(i,k)=(paprs(i,k)-paprs(i,k+1))/zrho/RG
-!      zalt(i,k+1)=zalt(i,k)+zdz(i,k)   !--in m
-!    ENDDO
-!  ENDDO
-!!
-!!--vertical reprojection 
-!!
-!  flight_m(:,:)=0.0 
-!  flight_h2o(:,:)=0.0 
-!!
-!  DO k=1, klev
-!    DO kori=1, nleva 
-!      DO i=1, klon
-!        !--fraction of layer kori included in layer k
-!        zfrac=max(0.0,min(zalt(i,k+1),zinta(kori+1))-max(zalt(i,k),zinta(kori)))/(zinta(kori+1)-zinta(kori))
-!        !--reproject
-!        flight_m(i,k)=flight_m(i,k) + pkm_airpl(i,kori,mth_cur)*zfrac 
-!        !--reproject
-!        flight_h2o(i,k)=flight_h2o(i,k) + ph2o_airpl(i,kori,mth_cur)*zfrac   
-!      ENDDO
-!    ENDDO
-!  ENDDO
-!!
-!  DO k=1, klev
-!    DO i=1, klon
-!      !--molec.cm-3.s-1 / (molec/mol) * kg CO2/mol * m2 * m * cm3/m3 / (kg CO2/m) => m s-1 per cell
-!      flight_m(i,k)=flight_m(i,k)/RNAVO*44.e-3*cell_area(i)*zdz(i,k)*1.e6/16.37e-3
-!      flight_m(i,k)=flight_m(i,k)*100.0  !--x100 to augment signal to noise 
-!      !--molec.cm-3.s-1 / (molec/mol) * kg H2O/mol * m2 * m * cm3/m3 => kg H2O s-1 per cell
-!      flight_h2o(i,k)=flight_h2o(i,k)/RNAVO*18.e-3*cell_area(i)*zdz(i,k)*1.e6
-!    ENDDO
-!  ENDDO
 !!
 !END SUBROUTINE airplane