Changeset 3195 for trunk/LMDZ.PLUTO/libf/dynphy_lonlat

Timestamp:

Jan 31, 2024, 4:36:51 PM (2 years ago)

Author:

afalco

Message:

Pluto PCM:
Imported condense n2 from pluto.old.
Aerosol data from Pluto.old not yet working.
AF

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/lect_start_archive.F (modified) (44 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/lect_start_archive.F

@@ -7 +7 @@
 
       USE comsoil_h, ONLY: nsoilmx, layer, mlayer, volcapa, inertiedat
-      USE tracer_h, ONLY: igcm_n2_ice
+      USE tracer_h, ONLY: igcm_haze
       USE infotrac, ONLY: tname, nqtot
 !      USE slab_ice_h, ONLY: noceanmx
@@ -39 +39 @@
 
 
-c Variables pour les lectures des fichiers "ini" 
+c Variables pour les lectures des fichiers "ini"
 c--------------------------------------------------
 !      INTEGER sizei,
@@ -53 +53 @@
 !      character (len=50) :: tmpname
 
-c Variable histoire 
+c Variable histoire
 c------------------
       REAL,INTENT(OUT) :: vcov(iip1,jjm,llm),ucov(iip1,jjp1,llm) ! vents covariants
@@ -59 +59 @@
       REAL,INTENT(OUT) :: q(iip1,jjp1,llm,nqtot)
 
-c Physique sur grille scalaire 
+c Physique sur grille scalaire
 c----------------------------
 
@@ -112 +112 @@
 
 
-c Variable de l'ancienne grille 
+c Variable de l'ancienne grille
 c---------------------------------------------------------
 
@@ -150 +150 @@
       logical :: therminertia_3D=.true. ! flag
 ! therminertia_3D=.true. if thermal inertia is 3D and read from datafile
-c Variable intermediaires iutilise pour l'extrapolation verticale 
+c Variable intermediaires iutilise pour l'extrapolation verticale
 c----------------------------------------------------------------
-      real, dimension(:,:,:), allocatable :: var,varp1 
+      real, dimension(:,:,:), allocatable :: var,varp1
       real, dimension(:), allocatable :: oldgrid, oldval
       real, dimension(:), allocatable :: newval
@@ -158 +158 @@
       real,intent(out) :: surfith(iip1,jjp1) ! surface thermal inertia
       ! surface thermal inertia at old horizontal grid resolution
-      real, dimension(:,:), allocatable :: surfithold 
+      real, dimension(:,:), allocatable :: surfithold
 
       character(len=30) :: txt ! to store some text
@@ -172 +172 @@
 !-----------------------------------------------------------------------
 
-! 1.2 Read the various dimension lengths of data in file 
+! 1.2 Read the various dimension lengths of data in file
 
       ierr= NF_INQ_DIMID(nid,"Time",dimid)
@@ -245 +245 @@
 
 ! 1.3 Report dimensions
-      
+
       write(*,*) "Start_archive dimensions:"
       write(*,*) "longitude: ",imold
@@ -259 +259 @@
       endif
       write(*,*) "time lenght: ",timelen
-      write(*,*) 
+      write(*,*)
 
 !-----------------------------------------------------------------------
@@ -311 +311 @@
 
 C-----------------------------------------------------------------------
-c 3.1. Lecture du tableau des parametres du run 
+c 3.1. Lecture du tableau des parametres du run
 c     (pour  la lecture ulterieure de "ptotalold" et "n2icetotalold")
 c-----------------------------------------------------------------------
@@ -355 +355 @@
          PRINT*, "lect_start_archive: Field <rlatu> not found"
          CALL abort
-      ENDIF 
+      ENDIF
 #ifdef NC_DOUBLE
       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlatuold)
@@ -440 +440 @@
 c 3.4 Read Soil layers depths
 c-----------------------------------------------------------------------
-     
+
       ierr=NF_INQ_VARID(nid,"soildepth",nvarid)
       if (ierr.ne.NF_NOERR) then
@@ -540 +540 @@
       ptotalold = tab_cntrl(tab0+49)
       n2icetotalold = tab_cntrl(tab0+50)
- 
+
 c-----------------------------------------------------------------------
 c 4. Lecture du temps et choix
 c-----------------------------------------------------------------------
- 
+
 c  lecture du temps
 c
@@ -559 +559 @@
       IF (ierr .NE. NF_NOERR) THEN
          ierr = NF_INQ_VARID (nid, "temps", nvarid)
-      endif 
+      endif
 #ifdef NC_DOUBLE
       ierr = NF_GET_VAR_DOUBLE(nid, nvarid, timelist)
@@ -584 +584 @@
 
    6  FORMAT(i7,i7,f9.3)
- 
+
       write(*,*)
       write(*,*) 'Choice for the date'
@@ -595 +595 @@
         endif
       end do
- 
+
       if (memo.eq.0) then
         write(*,*)
@@ -618 +618 @@
 c 5.1 Lecture des champs 2D (n2ice, emis,ps,tsurf,Tg[10], qsurf)
 c-----------------------------------------------------------------------
- 
+
       start=(/1,1,memo,0/)
       count=(/imold+1,jmold+1,1,0/)
-       
+
       ierr = NF_INQ_VARID (nid, "emis", nvarid)
       IF (ierr .NE. NF_NOERR) THEN
@@ -756 +756 @@
 !          PRINT*, "lect_start_archive: Failed loading <sea_ice>"
 !       ENDIF
-  
+
 !       ENDIF! ok_slab_ocean
 c
@@ -763 +763 @@
       write(*,*)
 
-! Surface tracers:      
+! Surface tracers:
       ! initialize all surface tracers to zero
       qsurfold(1:imold+1,1:jmold+1,1:nqtot)=0
@@ -777 +777 @@
           endif
 
-        
+
         write(*,*) "lect_start_archive: loading tracer ",trim(txt)
         ierr = NF_INQ_VARID (nid,txt,nvarid)
@@ -838 +838 @@
 c
        enddo ! of do isoil=1,nsoilold
-      
+
       ! reset 'start' and 'count' to "3D" behaviour
       start=(/1,1,1,memo/)
       count=(/imold+1,jmold+1,nsoilold,1/)
-      
+
       else
        write(*,*) "lect_start_archive: loading tsoil "
@@ -856 +856 @@
 #endif
        endif ! of if (ierr.ne.NF_NOERR)
-       
+
       endif ! of if (olddepthdef)
 
@@ -928 +928 @@
 c
 
-! Tracers:      
+! Tracers:
       qold(1:imold+1,1:jmold+1,1:lmold,1:nqtot)=0.
 
@@ -1040 +1040 @@
 c   INTERPOLATION DANS LA NOUVELLE GRILLE et initialisation des variables
 c=======================================================================
-c  Interpolation horizontale puis passage dans la grille physique pour 
-c  les variables physique 
+c  Interpolation horizontale puis passage dans la grille physique pour
+c  les variables physique
 c  Interpolation verticale puis horizontale pour chaque variable 3D
 c=======================================================================
@@ -1048 +1048 @@
 c 6.1   Variable 2d :
 c-----------------------------------------------------------------------
-c Relief 
+c Relief
       call interp_horiz (phisold,phisold_newgrid,imold,jmold,iim,jjm,1,
      &                   rlonuold,rlatvold,rlonu,rlatv)
 
-! N2 ice is now in qsurf(igcm_n2_ice)
+! N2 ice is now in qsurf(igcm_haze)
 !      call interp_horiz (n2iceold,n2ices,imold,jmold,iim,jjm,1,
 !     &                   rlonuold,rlatvold,rlonu,rlatv)
@@ -1088 +1088 @@
 c       On assure la conservation de la masse de l'atmosphere + calottes
 c-----------------------------------------------------------------------
+      !AF: TODO: mass conservation: check this. haze?
 
       ptotal =  0.
@@ -1096 +1097 @@
       END DO
       n2icetotal = 0.
-      if (igcm_n2_ice.ne.0) then
+      if (igcm_haze.ne.0) then
         ! recast surface N2 ice on new grid
-        call interp_horiz(qsurfold(1,1,igcm_n2_ice),
-     &                  qsurfs(1,1,igcm_n2_ice),
+        call interp_horiz(qsurfold(1,1,igcm_haze),
+     &                  qsurfs(1,1,igcm_haze),
      &                  imold,jmold,iim,jjm,1,
      &                  rlonuold,rlatvold,rlonu,rlatv)
@@ -1105 +1106 @@
          DO i=1,iim
            !n2icetotal = n2icetotal + n2iceS(i,j)*aire(i,j)
-           n2icetotal=n2icetotal+qsurfS(i,j,igcm_n2_ice)*aire(i,j)
+           n2icetotal=n2icetotal+qsurfS(i,j,igcm_haze)*aire(i,j)
          END DO
        END DO
@@ -1115 +1116 @@
       write(*,*)'Old grid: atmospheric mass :',ptotalold
       write(*,*)'New grid: atmospheric mass :',ptotal
-      write (*,*) 'Ratio new atm./ old atm =', ptotal/ptotalold 
+      write (*,*) 'Ratio new atm./ old atm =', ptotal/ptotalold
       write(*,*)
       write(*,*)'Old grid: mass of N2 ice:',n2icetotalold
@@ -1131 +1132 @@
       END DO
 
-      if ( n2icetotalold.gt.0.) then 
+      if ( n2icetotalold.gt.0.) then
 !         DO j=1,jjp1
 !            DO i=1,iip1
@@ -1157 +1158 @@
      &f soil thermal inertia; might be wiser to reset it.'
         write(*,*)
-       
+
         do i=1,imold+1
          do j=1,jmold+1
@@ -1177 +1178 @@
         ! We have inertiedatold
        if((imold.ne.iim).or.(jmold.ne.jjm)) then
-       write(*,*)'lect_start_archive: WARNING: horizontal interpolation 
+       write(*,*)'lect_start_archive: WARNING: horizontal interpolation
      &of thermal inertia; might be better to reset it.'
        write(*,*)
        endif
-       
+
         ! Do horizontal interpolation
         if (depthinterpol) then
@@ -1210 +1211 @@
       call gr_dyn_fi (nsoilmx,iim+1,jjm+1,ngrid,
      &                  inertiedatS,inertiedat)
-      
+
 c-----------------------------------------------------------------------
 c 6.2.2 Soil temperature
@@ -1282 +1283 @@
         deallocate(oldval)
         deallocate(newval)
-       
+
        else
         tsoiloldnew(:,:,:)=tsoilold(:,:,:)
@@ -1303 +1304 @@
 c 6.4 Variable 3d :
 c-----------------------------------------------------------------------
-      
+
 c temperatures atmospheriques
       write (*,*) 'lect_start_archive: told ', told (1,jmold+1,1)  ! INFO
@@ -1321 +1322 @@
      &                   rlonuold,rlatvold,rlonu,rlatv)
       call gr_dyn_fi(llm,iim+1,jjm+1,ngrid,du_nonoro_gwdS,du_nonoro_gwd)
-      
+
       call interp_vert
      &    (dv_nonoro_gwdold,var,lmold,llm,apsold,bpsold,aps,bps,
@@ -1328 +1329 @@
      &                   rlonuold,rlatvold,rlonu,rlatv)
       call gr_dyn_fi(llm,iim+1,jjm+1,ngrid,dv_nonoro_gwdS,dv_nonoro_gwd)
-      
+
       call interp_vert
      &    (east_gwstressold,var,lmold,llm,apsold,bpsold,aps,bps,
@@ -1335 +1336 @@
      &                   rlonuold,rlatvold,rlonu,rlatv)
       call gr_dyn_fi(llm,iim+1,jjm+1,ngrid,east_gwstressS,east_gwstress)
-      
+
       call interp_vert
      &    (west_gwstressold,var,lmold,llm,apsold,bpsold,aps,bps,
@@ -1379 +1380 @@
           end do
         end do
-      end do 
+      end do
       write (*,*) 'lect_start_archive: ucov ', ucov (1,2,1)  ! INFO
 c     write(48,*) 'ucov',ucov
@@ -1410 +1411 @@
      &                  rlonuold,rlatvold,rlonu,rlatv)
       enddo
-cccccccccccccccccccccccccccccc      
-c  make sure that sum of q = 1      
-c dominent species is = 1 - sum(all other species)      
-cccccccccccccccccccccccccccccc      
+cccccccccccccccccccccccccccccc
+c  make sure that sum of q = 1
+c dominent species is = 1 - sum(all other species)
+cccccccccccccccccccccccccccccc
 c      iqmax=1
-c      
+c
 c      if (nqold.gt.10) then
 c       do l=1,llm
@@ -1428 +1429 @@
 c           qtot(i,j,l)=0
 c           do iq=1,nqold
-c            if (iq.ne.iqmax) then        
-c              q(i,j,l,iqmax)=q(i,j,l,iqmax)-q(i,j,l,iq)        
+c            if (iq.ne.iqmax) then
+c              q(i,j,l,iqmax)=q(i,j,l,iqmax)-q(i,j,l,iq)
 c            endif
 c           enddo !iq
@@ -1437 +1438 @@
 c     $    qtot(i,j,l)
 c           enddo !iq
-c          enddo !i   
-c         enddo !j   
-c       enddo !l  
+c          enddo !i
+c         enddo !j
+c       enddo !l
 c      endif
 ccccccccccccccccccccccccccccccc
@@ -1451 +1452 @@
          end do
       enddo
-      
+
 !      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,n2ices,n2ice)
-! no need to transfer "n2ice" any more; it is in qsurf(igcm_n2_ice)
+! no need to transfer "n2ice" any more; it is in qsurf(igcm_haze)
 
       endif !! if nqtot .ne. 0