[3006] | 1 | module nirdata_mod |
---|
[414] | 2 | |
---|
[3006] | 3 | implicit none |
---|
[414] | 4 | |
---|
[3006] | 5 | ! Number of pressures in NIR correction table |
---|
| 6 | integer,parameter :: npres=42 |
---|
| 7 | |
---|
| 8 | real,save,protected :: pres1d(npres) |
---|
| 9 | real,save,protected :: corgcm(npres) |
---|
| 10 | real,save,protected :: oco21d(npres) |
---|
| 11 | real,save,protected :: alfa(npres) |
---|
| 12 | real,save,protected :: p1999(npres) |
---|
| 13 | !$OMP THREADPRIVATE(pres1d,corgcm,oco21d,alfa,p1999) |
---|
| 14 | |
---|
| 15 | contains |
---|
| 16 | |
---|
| 17 | subroutine NIR_leedat |
---|
| 18 | |
---|
| 19 | ! reads parameters for NIR NLTE calculation |
---|
| 20 | |
---|
| 21 | ! nov 2011 fgg+malv first version |
---|
| 22 | !*********************************************************************** |
---|
| 23 | |
---|
| 24 | use datafile_mod, only: datadir |
---|
| 25 | USE mod_phys_lmdz_para, ONLY: is_master |
---|
| 26 | USE mod_phys_lmdz_transfert_para, ONLY: bcast |
---|
| 27 | |
---|
| 28 | implicit none |
---|
| 29 | |
---|
| 30 | ! local variables |
---|
| 31 | |
---|
| 32 | integer :: ind |
---|
| 33 | |
---|
| 34 | if (is_master) then |
---|
| 35 | open(43,file=trim(datadir)//'/NIRcorrection_feb2011.dat', & |
---|
| 36 | status='old') |
---|
| 37 | do ind=1,9 |
---|
| 38 | read(43,*) |
---|
| 39 | enddo |
---|
| 40 | |
---|
| 41 | do ind=1,npres |
---|
| 42 | read(43,*)pres1d(ind),corgcm(ind),oco21d(ind),p1999(ind), & |
---|
| 43 | alfa(ind) |
---|
| 44 | !Tabulated pression to Pa |
---|
| 45 | pres1d(ind)=pres1d(ind)*100. |
---|
| 46 | enddo |
---|
| 47 | close(43) |
---|
| 48 | |
---|
| 49 | endif ! if(is_master) then |
---|
| 50 | |
---|
| 51 | call bcast(pres1d) |
---|
| 52 | call bcast(corgcm) |
---|
| 53 | call bcast(oco21d) |
---|
| 54 | call bcast(p1999) |
---|
| 55 | call bcast(alfa) |
---|
| 56 | |
---|
| 57 | end subroutine NIR_leedat |
---|
| 58 | |
---|
| 59 | end module nirdata_mod |
---|