source: trunk/LMDZ.PLUTO.old/libf/phypluto/lymalpha.F90 @ 3436

Last change on this file since 3436 was 3175, checked in by emillour, 11 months ago

Pluto PCM:
Add the old Pluto LMDZ for reference (required prior step to making
an LMDZ.PLUTO using the same framework as the other physics packages).
TB+EM

File size: 1.6 KB
RevLine 
[3175]1      subroutine lymalpha(pls,pflux)
2      use datafile_mod   
3      implicit none
4
5!==================================================================
6!     Purpose
7!     -------
8!     Get the lyman alpha flux according to the solar longitude
9
10!     Inputs
11!     ------
12!     pls                 solar longitude
13!     
14!     Outputs
15!     -------
16!     
17!     Both
18!     ----
19!
20!     Authors
21!     -------
22!     Tanguy Bertrand     
23!==================================================================
24
25#include "comcstfi.h"
26!-----------------------------------------------------------------------
27!     Arguments
28
29      REAL,INTENT(IN) :: pls
30      REAL,INTENT(OUT) :: pflux
31!-----------------------------------------------------------------------
32!     Local variables
33      REAL :: vectls
34      REAL :: vectflux
35      LOGICAL firstcall
36      SAVE firstcall
37      DATA firstcall/.true./
38
39      !!read lyman alpha flux
40      integer Nfine
41      parameter(Nfine=13281)
42      integer ifine
43      character(len=100) :: file_path
44      real,save :: lsdat(Nfine),fluxdat(Nfine)
45
46!---------------- INPUT ------------------------------------------------
47     
48      IF (firstcall) then
49        firstcall=.false.
50        file_path=trim(datadir)//'/sol_uv_flux.txt'
51        print*,file_path
52        open(222,file=file_path,form='formatted')
53        do ifine=1,Nfine
54           read(222,*) lsdat(ifine), fluxdat(ifine)
55        enddo
56        close(222)
57      ENDIF
58
59      CALL interp_line(lsdat,fluxdat,Nfine,pls*180./pi,pflux,1)
60      !write(*,*) 'flux=',pflux
61     
62
63      end subroutine lymalpha
64
Note: See TracBrowser for help on using the repository browser.