! ! $Id $ ! SUBROUTINE dynetat0(fichnom,vcov,ucov, . teta,q,masse,ps,phis,time) USE infotrac use netcdf, only: nf90_get_var use control_mod, only : planet_type USE comconst_mod, ONLY: cpp, daysec, dtvr, g, im, jm, kappa, & lllm, omeg, rad USE serre_mod, ONLY: clon,clat,grossismx,grossismy IMPLICIT NONE c======================================================================= c c Auteur: P. Le Van / L.Fairhead c ------- c c objet: c ------ c c Lecture de l'etat initial c c======================================================================= c----------------------------------------------------------------------- c Declarations: c ------------- #include "dimensions.h" #include "paramet.h" #include "temps.h" #include "comvert.h" #include "comgeom2.h" #include "ener.h" #include "netcdf.inc" #include "description.h" #include "logic.h" #include "iniprint.h" c Arguments: c ---------- CHARACTER*(*) fichnom REAL vcov(iip1, jjm,llm),ucov(iip1, jjp1,llm),teta(iip1, jjp1,llm) REAL q(iip1,jjp1,llm,nqtot),masse(iip1, jjp1,llm) REAL ps(iip1, jjp1),phis(iip1, jjp1) REAL time c Variables c INTEGER length,iq PARAMETER (length = 100) REAL tab_cntrl(length) ! tableau des parametres du run INTEGER ierr, nid, nvarid INTEGER idecal c----------------------------------------------------------------------- c Ouverture NetCDF du fichier etat initial ierr = NF_OPEN (fichnom, NF_NOWRITE,nid) IF (ierr.NE.NF_NOERR) THEN write(lunout,*)'dynetat0: Pb d''ouverture du fichier start.nc' write(lunout,*)' ierr = ', ierr CALL ABORT_gcm("dynetat0", "", 1) ENDIF c ierr = NF_INQ_VARID (nid, "controle", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, tab_cntrl) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echoue pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF !!! AS: idecal is a hack to be able to read planeto starts... !!! .... while keeping everything OK for LMDZ EARTH if (planet_type.eq."generic") then print*,'NOTE NOTE NOTE : Planeto-like start files' idecal = 4 annee_ref = 2000 else print*,'NOTE NOTE NOTE : Earth-like start files' idecal = 5 annee_ref = tab_cntrl(5) endif im = tab_cntrl(1) jm = tab_cntrl(2) lllm = tab_cntrl(3) day_ref = tab_cntrl(4) rad = tab_cntrl(idecal+1) omeg = tab_cntrl(idecal+2) g = tab_cntrl(idecal+3) cpp = tab_cntrl(idecal+4) kappa = tab_cntrl(idecal+5) daysec = tab_cntrl(idecal+6) dtvr = tab_cntrl(idecal+7) etot0 = tab_cntrl(idecal+8) ptot0 = tab_cntrl(idecal+9) ztot0 = tab_cntrl(idecal+10) stot0 = tab_cntrl(idecal+11) ang0 = tab_cntrl(idecal+12) pa = tab_cntrl(idecal+13) preff = tab_cntrl(idecal+14) c clon = tab_cntrl(idecal+15) clat = tab_cntrl(idecal+16) grossismx = tab_cntrl(idecal+17) grossismy = tab_cntrl(idecal+18) c IF ( tab_cntrl(idecal+19).EQ.1. ) THEN fxyhypb = . TRUE . c dzoomx = tab_cntrl(25) c dzoomy = tab_cntrl(26) c taux = tab_cntrl(28) c tauy = tab_cntrl(29) ELSE fxyhypb = . FALSE . ysinus = . FALSE . IF( tab_cntrl(idecal+22).EQ.1. ) ysinus = . TRUE. ENDIF day_ini = tab_cntrl(30) itau_dyn = tab_cntrl(31) start_time = tab_cntrl(32) c ................................................................. c c write(lunout,*)'dynetat0: rad,omeg,g,cpp,kappa', & rad,omeg,g,cpp,kappa IF( im.ne.iim ) THEN PRINT 1,im,iim STOP ELSE IF( jm.ne.jjm ) THEN PRINT 2,jm,jjm STOP ELSE IF( lllm.ne.llm ) THEN PRINT 3,lllm,llm STOP ENDIF ierr = NF_INQ_VARID (nid, "rlonu", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, rlonu) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "rlatu", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, rlatu) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "rlonv", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, rlonv) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "rlatv", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, rlatv) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour rlatv" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "cu", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, cu) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "cv", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, cv) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "aire", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, aire) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "phisinit", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = nf90_get_var(nid, nvarid, phis) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Lecture echouee pour " CALL ABORT_gcm("dynetat0", "", 1) ENDIF ierr = NF_INQ_VARID (nid, "temps", nvarid) IF (ierr .NE. NF_NOERR) THEN write(lunout,*)"dynetat0: Le champ est absent" write(lunout,*)"dynetat0: J essaie