Changeset 4092 for trunk/LMDZ.VENUS/libf/phyvenus/nirdata.F90

Timestamp:

Feb 28, 2026, 6:36:58 PM (6 weeks ago)

Author:

emillour

Message:

Venus PCM:
Code cleanup: remove "hedin.h" and make variables included there be module
variables of compo_hedin_mod2.F90. Turn "nirdat.h" into module nirdata.F90
and incroporate initialization routine nir_leedat.F in it.
EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/nirdata.F90 (moved) (moved from trunk/LMDZ.VENUS/libf/phyvenus/nirdata.h) (1 diff)

trunk/LMDZ.VENUS/libf/phyvenus/nirdata.F90

@@ +1 @@
+module nirdata_mod
 
-      integer npres                ! Number of pressures in NIR correction
-      parameter (npres=42)         ! table
+implicit none
 
-      common /NIRdata/ pres1d,corgcm,oco21d,alfa,p1999
-      real    pres1d(npres)
-      real    corgcm(npres)
-      real    oco21d(npres),alfa(npres),p1999(npres)
+! Number of pressures in NIR correction table
+integer,parameter :: npres=42   
+
+real,save,protected :: pres1d(npres)
+real,save,protected :: corgcm(npres)
+real,save,protected :: oco21d(npres)
+real,save,protected :: alfa(npres)
+real,save,protected :: p1999(npres)
+!$OMP THREADPRIVATE(pres1d,corgcm,oco21d,alfa,p1999)
+
+contains
+
+  subroutine nir_leedat 
+
+!       reads parameters for NIR NLTE calculation    
+
+!       nov 2011    fgg+malv    first version                
+!***********************************************************************
+
+    use mod_phys_lmdz_para, only: is_master, bcast
+
+    implicit none
+  
+    integer :: ind
+    character (len=100) :: datafile="HIGHATM"
     
+    if (is_master) then ! only the master needs read the file
+      open(43,file=trim(datafile)//'/NIRcorrection_feb2011.dat', &
+              status='old')
+      do ind=1,9
+         read(43,*)
+      enddo
+      
+      do ind=1,npres
+         read(43,*)pres1d(ind),corgcm(ind),oco21d(ind),p1999(ind), &
+                   alfa(ind)
+         ! Tabulated pressure to Pa
+         pres1d(ind)=pres1d(ind)*100.
+      enddo
+      
+      close(43)
+    
+    endif ! of if (is_master)
+  
+    ! broadcast to all cores
+    call bcast(pres1d)
+    call bcast(corgcm)
+    call bcast(oco21d)
+    call bcast(p1999)
+    call bcast(alfa)
+
+  end subroutine nir_leedat 
+
+end module nirdata_mod